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ABSTRACT 


The anticipated application of advanced turboprop 
propulsion systems and use of composite materials in primary 
structure is expected to increase the interior noise of 
future aircraft to unacceptably high levels. The absence of 
technically and economically feasible noise source-path 
diagnostic tools has been a prime obstacle in the develop- 
ment of efficient noise control treatments for propeller 
driven aircraft. A new diagnostic method which permits the 
separation and prediction of the fully coherent airborne and 
structureborne components of the sound radiated by plates or 
thin shells has been developed. Analytical and experimental 
studies of the proposed method were performed on plates con- 
structed of both conventional and composite materials. The 
results of the study indicate that the proposed method can 
be applied to a variety of aircraft materials, could be used 
in flight, and has fewer encumbrances than the other diag- 
nostic tools currently available. The study has also re- 
vealed that the noise radiation of vibrating plates in the 
low frequency regime due to combined airborne and structure- 
borne inputs possesses a strong synergistic nature. The 
large influence of the interaction between the airborne and 
structureborne terms has been hitherto ignored by research- 
ers of aircraft interior noise problems. 
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Chapter I 
INTRODUCTION 

A_j_ INTERIOR NOISE IN PROPELLER DRIVEN AIRCRAFT 

Interior noise levels of propeller driven aircraft have 
been measured between 84 and 104 dB on the A-weighted scale. 
(See references 1 and 2.) Limited exposure to these noise 
levels can cause a temporary shift in the hearing threshold 
of the passengers. Prolonged exposure can result in perm- 
anent hearing damage. (Pilot and crew are at risk.) The 
noise can also interfere with the operational safety of the 
aircraft and the efficiency of the pilot and crew. For 
these reasons, high interior noise levels in propeller dri- 
ven aircraft have historically been a cause for concern in 
both the commercial and the military sectors of the aircraft 
industry. 

Problems with interior noise levels in future aircraft 
are expected to intensify due to the advanced turboprop pro- 
pulsion systems now being incorporated into the design of 
transport aircraft and because of the anticipated widespread 
use of composite materials in the primary structure of air- 
craft. (See Appendix I for a discussion of how the future 
trends in aircraft design are expected to effect interior 
noise levels in aircraft. ) (Also see references 1-27 for 
background information relating to interior noise in pro- 
eller driven aircraft. ) 
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The noise entering the cabin of a propeller driven air- 
craft is generally divided into two major categories, viz. 
airborne noise and structureborne noise. Airborne noise is 
generated aerodynamically by the propellers and propagates 
along an acoustic path (the fluid medium being air in this 
case). This noise then strikes the sidewalls of the air- 
craft and is transmitted into the interior of the aircraft. 
Structureborne noise has its source in the vibration of the 
wings or other structural members of the aircraft. This 
vibrational energy propagates along structural paths into 
the cabin of the aircraft where it causes vibration of the 
sidewalls and other surfaces inside the aircraft. This 
structural vibration then radiates noise inside the aircraft 
cabin . 

B. REDUCTION OF INTERIOR NOISE IN PROPELLER DRIVEN AIRCRAFT 
The distinction between airborne noise and structure- 
borne noise is important because the methods typically used 
in reducing airborne and structureborne noise are quite 
different. For example, if the predominant source of the 
noise in the aircraft is structureborne, then the problem 
might be solved through the use of vibration isolators, or 
by the application of damping materials. If the predominant 
source of noise is airborne, then the problem might be solv- 
ed by aft mounting the propellers or by adding massive mat- 
erials to the sidewalls. 
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Since much of the success or failure of a particular 
aircraft design depends on the performance of the aircraft 
in terms of the payload capacity, range, and other weight 
sensitive parameters, an aircraft designer/manufacturer can 
add mass to an aircraft for noise control purposes only on 
the areas of the fuselage where it is absolutely necessary. 
Thus, to optimize the noise control design and treatment of 
an aircraft from a cost/performance standpoint, a confident 
knowledge of the relative importance of the airborne and 
structureborne noise transmission is essential. This know- 
ledge can be obtained only if the aircraft designer and 
manufacturer has reliable and effective noise source/path 
diagnostic tools. Unfortunately, the diagnostic measurement 
and prediction methods currently available to the aircraft 
industry for determining the relative contributions of air- 
borne and structureborne noise have proven themselves either 
inadequate or technically and/or economically unfeasible. 
(See Appendix I for a brief review of the most recent 
diagnostic methods applied to aircraft.) (Also see refer- 
ences 1-27 for background information.) 

THE PURPOSE OF THIS RESEARCH 

Recently, several new noise source/path identification 
tools have come into widespread use. Among the most promis- 
ing of these new tools are several methods for measuring the 
acoustic intensity vector. (See references 28-34.) In the 
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past 5 years, researchers have begun apply these methods to 
noise transmission problems in aircraft. (See references 
35-41.) The purpose of this research is to develop a new 
measurement method, based on the two microphone cross spec- 
tral acoustic intensity measurement method, for separating 
and predicting the airborne and structureborne components of 
the noise radiated by aircraft- type panels. 
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Chapter II 
PROBLEM APPROACH 

OVERVIEW OF THE DIAGNOSTIC APPROACH 

The measurement method proposed in this study utilizes 

a two microphone acoustic intensity probe to measure the 

time averaged sound power radiated, <TT> t , and an array 

of miniature accelerometers to measure the space-time aver- 

2 

aged mean square surface velocity, <v > . , of a thin 

r , t 

shell structure. The method employs these two types of 
measurements (exclusively) along with several equations 
developed by the author (based on theoretical consider- 
ations) to separate and predict the relative amounts of the 
total (combined) sound power that can be attributed to the 
airborne and structureborne components (respectively) due to 
some unknown combination of acoustic and vibrational inputs. 
The scheme of the proposed diagnostic method is to: 

(1) determine the radiation efficiency of the structure when 
it is radiating purely airborne noise (a ) , 

(2) determine the radiation efficiency of the structure when 
it is radiating purely structureborne noise (o s ), 

(3) use the information obtained in steps (1) and (2) along 

with measurements of <7T>. and <v > . (while the 

L IT j L* 

structure is being driven by the combined noise sources of 
interest) to predict the airborne and structureborne compo- 
nents of the combined sound power. 
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B. PROOF OF CONCEPT 


In order to validate the proposed diagnostic method it 
was decided that both analytical and experimental studies 
should be performed. The purpose of the analytical study 
was to develop a theoretical foundation for the proposed 
diagnostic method (based on the theory of sound radiation of 
thin shell structures) and to simulate the performance of 
the diagnostic method under various operating conditions. 

The purpose of the experimental study was then to verify the 
expected behavior of the diagnostic method by showing that 
the trends predicted by the analytical simulations are the 
same trends that are observed under actual measurement cond- 
itions . 

It was decided that simple rectangular plates would 
serve as the test vehicle for both analytical simulations 
and experimental verification of the proposed diagnostic 
method. Plates were chosen because they possess most of the 
vibrational and sound radiative properties that are exhi- 
bited by actual aircraft sidewalls (due to their thin shell 
construction). Since the accuracy of the proposed measure- 
ment method does not depend on the complexity of structures 
used, if the method works well on simple flat plates, it 
should also work on more complicated structures. Futher- 
more, the theory of vibration and sound radiation of plates 
is well developed enough so that the author could extend it 
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to include the case of combined airborne and structureborne 
inputs. Last, but not least, the plates could easily be 
constructed and mounted with the experimental facilities 
that were available. 

The physical dimensions of the plates were chosen to be 
.4064 m x .2413 m (16 in x 9.5 in) since an extensive study 
of the transmission loss properties of plates of this size 
had already been completed by NASA researchers. Both con- 
ventional aluminum and composite materials were included in 
the study since it is expected that composites will be used 
in the primary structure of aircraft built in the 1990's and 
beyond. 

G. BASIS FOR THE ANALYTICAL AND EXPERIMENTAL STUDIES 

An analytical model of the vibration and sound radia- 
tion of plates of symmetric angle-ply laminate construction 
due to combined airborne and structureborne inputs was dev- 
eloped. Chapter III, which follows, will outline the devel- 
opment of the analytical model and present the theoretical 
basis for the proposed diagnostic method. (See Appendices 
II-VI for mathematical details of the theory used in dev- 
eloping the analytical model.) 

The analytical model assumes that the plates are rect- 
angular and simply supported in an infinite., rigid baffle. 
The plates are assumed to be symmetric laminates with multi- 
ple generally orthotropic layers (e.g. graphite, fiberglass, 
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or Kevlar regular symmetric angle-ply tape panels) or con- 
structed of simple orthotropic materials (e.g. isotropic A A 
2024 aluminum panels). The sound radiated by the plates is 
assumed to be generated by the flexural (bending) vibration- 
al response of the plate and is radiated to a free-field 
acoustic space (anechoic). Futhermore, the dynamic response 
of the plate (which generates the sound) is assumed to be 
dependent only on the incident airborne and structureborne 
forcing functions on the plate and independent of (uncoupled 
from) the sound pressure radiated by the plate. (This last 
assumption is sometimes referred to as the blocked pressure 
assumption in acoustic transmission theory.) 

A normally incident, spatially uniform pressure field 
was used to model the airborne input. A point load was used 
to model the structureborne input. The analysis was confin- 
ed to the 0-1000 Hz frequency range since this range encom- 
passes the most troublesome noise region for propeller dri- 
ven aircraft. 

The experimental study of the proposed diagnostic 
method was performed using the NASA Langley Research Cen- 
ter's acoustic transmission loss apparatus. Using this 
facility, the test panels could be subjected to the desired 
acoustic and vibrational inputs in the source room while the 
surface velocity and sound power radiated by the panels 
could be measured in the receiving room. Every effort was 
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made to emulate the simplifying assumptions used in the 
analytical modeling. Chapter IV, which follows, describes 
the special equipment that was constructed for mounting and 
driving the plates and the modifications that were made to 
the transmission loss apparatus. 

The results of the analytical simulations of the appli- 
cation of the proposed diagnostic method under various oper- 
ating conditions and the experimental verification of the 
predicted trends are presented and discussed in chapter V. 
The variation of several parameters and their effects on the 
proposed diagnostic method are presented for both the analy- 
tical and experimental investigations. Parameters studied 
include the effects on the diagnostic method due to changes 
in the relative magnitude and phase of the airborne and 
structureborne inputs, changes in the level of structural 
damping, changes in the location of the structureborne 
input, and changes in the thickness, density, and fiber 
orientation for the composites. 
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Chapter III 
ANALYSIS 

A_;_ GOVERNING DIFFERENTIAL EQUATION OF THE PLATE DYNAMICS 

The governing equation for the transverse undamped free 
vibration of a flat rectangular orthotropic plate is given 
by 

D ll3 4 u/3aJ + 2(D 12 +2D g6 )9 4 u/3a^3a2 + D 22 3 4 u/3a 4 

+ phu = 0 , (3.1) 

where and a 2 are the in-plane coordinates, u is 
the transverse displacement of the plate and the D^ are 
the plate bending rigidity constants. (See Appendix II for 
a more complete discussion of the parameters in this equa- 
tion.) This equation can be solved exactly for the case of 
simply supported boundary conditions by using the method of 
separation of variables. From Appendix II, the solution is 
given by 

u( a ]_ ’ a 2 ’ t ) = 

“ “ j(<D t+e ) 

I I Ti mn sin(m7ra 1 /d 1 )si.n(mra 2 /d 2 )e mn mn (3.2) 

m = i n= i 

where the w are the radian natural frequencies of 
mn 

the plate given by 
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(3.3) 


“mn ' ,4/ » h 


D^m/dj.)' 


+ 2(D 12 +2D 66 )(m/d 1 ) 2 (n/d 2 ) 2 + D 22 (n/d 2 ) 4 


The governing equation for the transverse undamped free 
vibration of a flat rectangular plate of symmetric angle-ply 
laminate construction contains two additional terms due to 
the coupling between shear stress and extensional deforma- 
tion and coupling between normal stresses and shear deforma- 
tion and is given by (see reference 42, pg. 273) 


D n9 4 u/aaJ 


4D 16 3 4 u/9a^3a 2 + 


2(D 12 +2D 66 )3 4 u/3a 2 3a 2 


4D 26 3 4 u/3a 1 3a 2 


D 2 2 aV3a 4 


+ phu = 0. 


(3.4) 


(See Appendix VII for a complete discussion on the calcu- 
lation of the bending rigidity constants for sym- 

metric laminates with multiple generally orthotropic 
layers . ) 

Equation (3.4) can not be solved exactly using the 
method of separation of variables because of the additional 
fourth order differentials. In fact, no exact solution of 
any kind has been found for equation (3.4). It can be shown 
by dimensional analysis (see reference 43), however, that an 
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equation analogous to the equation (3.3) for the radian 
natural frequencies of the plate is given by 


2 4 / , 

“.n = ' /ph 


D^Cm/d^ 


- C(D 16 )(m/d 1 ) 3 (n/d 2 )+ 2(D 12 +2D 66 )(m/d 1 ) 2 (n/d 2 ) 2 


C(D 26 )(m/d 1 )(n/d 2 ) 3 + D 22 (n/d 2 ) 4 


(3.5) 


where C is some unknown constant. Bert (reference 43) 
utilized results obtained from the Raleigh-Ritz analysis 
performed by Ashton and Waddoups (reference 44) and results 
obtained from Green's classical Fourier analysis performed 
by Whitney (reference 45) to estimate an upper and lower 
bound on the value of the constant C. With these two anal- 
yses it was estimated that a value of C=2 was a good approx- 
imation for the unknown constant of equation (3.5). From 
this analysis performed by Bert (reference 43), it was 
decided that an estimate of the natural frequencies for a 
regular symmetric angle-ply laminate is given by 


2 4/ . 

“mn = 1 /ph 


D^fm/d^V 


2 (D 16 ) (m/d x ) 3 ( n/d 2 ) + 2 (D 12 +2D 66 ) (m/d L ) 2 ( n/d 2 ) 2 


- 2(D 26 ) (m/dj^) (n/d 2 ) 3 + D 22 (n/d 2 ) 4 


(3.6) 
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Thus, the solution of the transverse undamped f ree-vibration 
problem governed by equation (3.4) is approximated by equa- 
tion (3.2) with the natural frequencies, w , calcul- 

mn 

ated from equation (3.6). Hence, the analytical model 
chosen for this work treats the regular symmetric angle-ply 
laminates in exactly the same fashion as the simpler ortho- 
tropic plate with the only difference being that the natural 
frequencies of the laminates are calculated using equation 
(3.6). It should be noted that this model approximates both 
the eigenvalues and the eigenfunctions of equation (3.4) 
since the sine waves of equation (3.2) satisfy the geometric 
boundary conditions but do not satisfy the natural boundary 
conditions or the differential equation. These two approxi- 
mations will cause no error in the analysis of the noise 
radiation of symmetric laminates with multiple specially 
orthotropic layers (e.g. a 0/90 lay-up) since, for that 
case, and D are zero and the problem reduces 

to the orthotropic case. Since it is known that the true 
eigenfunctions are not exactly sine waves for a symmetric 
laminate with multiple generally orthotropic layers (e.g. a 
+45/-45 lay-up), the use of equation (3.2) can be expected 
to cause some inaccuracies in the calculated dynamic res- 
ponse, and subsequently the noise radiation in those cases 
(see reference 46, pp. 101-110). Inasmuch as the behavior 
of the symmetric angle-ply plates converges to the 
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orthotropic solution as the number of plies become large, 
the sine wave approximation is considered acceptable for 
plates with 8 or 16 plies. 

THE ADDITION OF DAMPING TO THE DYNAMIC MODEL 

The governing equation for the transverse damped free 
vibration of a flat, rectangular orthotropic plate is given 
by (see Appendix II) 

D n9 4 u/3aJ + 2(D 12 + 2D 66 )ci 4 u/3c^3a2 + D 22 3 4 u/3a 4 

+Xu + phu = 0 (3.7) 

where X is an equivalent viscous damping coefficient. 

We now make the assumption that the vibrational modes 
of the plate are uncoupled. Technically, the flexural 
modeshapes of a rectangular plate are uncoupled only when 
the damping is related to the mass and stiffness properties 
of the plate in a special fashion (see reference 47, pp. 
390-394). Thus, by making this assumption, our mathematical 
relations for the orthotropic plate become approximations. 
The assumption that the nodes are uncoupled can be used 
without serious error when the damping is a second order 
effect. Furthermore, the assumption is aided by the fact 
that the modal density for flexural modes of a rectangular 
plate is nearly constant at low frequency. This means that 
there are no frequency regions of high modal density and 
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therefore little opportunity for the exchange of energy 
among modes. With this assumption it is shown in Appendix 
II that the solution to equation (3.7) is given by 

u(a 1 ,a 2 ,t) = 

CD CD 

5 1 h mn sin(mTTa 1 /d 1 )sin(mra 2 /d 2 ) 

m= i n= i 

Wmrt j (^ m t + 0 ) 

mn ran J mn mn , n 

• e e . (3.8) 


And the damped radian natural frequencies Q mn are 
given by 


=03 

mn mn 


1 - 5 . 


1/2 


mn 


(3.9) 


where the modal damping coefficient is defined by 


E = X/(2pho) ) . 
mn mn 


(3.10) 


Once it is assumed that the modal damping coefficients 
are independent of one another, a value must be chosen for 
each. Reasonable damping coefficients for 2024 aluminum are 
in the range of E = 0.01. A more moderately damped struc- 
ture might have damping coefficients in the range E = 

0.04. Because of the large energy loss due to sound radia- 
tion at low frequency, however, experimentally determined 
damping coefficients for the first few resonance frequencies 
will appear to be much, much larger than this. (On the 
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order of E, = 0.1.) These inflated damping coefficients 
are a direct result of the blocked pressure assumption dis- 
cussed earlier. Thus, the large damping coefficients at low 
resonance frequencies do not owe their existence to any 
erroneous assumptions about the nature of the damping in the 
structure, but rather are the result of neglecting the sig- 
nificant fluid loading effect in the model of the forcing 
function. For these reasons, two different damping models 


were 

investigated 

in 

this study given by 

the equations 

^mn 

• 09 ( “ll A Vn ) 

+ 

.01 , 

(3.11) 

^mn 

•°6C“ll/“ in n ) 

4 * 

.04 , 

(3.12) 


where the hyberbolic term in each equation corrects for the 
fluid loading effect at the first few resonance frequencies 
and the constant term in each equation represents the inher- 
ent damping of the structure. Equation (3.11) shall be sub- 
sequently referred to as the small damping model while equa- 
tion (3.12) shall be referred to as the moderate damping 
model. A plot of these two damping models for a plate con- 
structed of AA 2024 aluminum is shown in figure (III-1). 
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C. DEFINITION OF THE FORCING FUNCTION FOR THE DYNAMIC MODEL 


The governing differential equation for the forced 
dynamic response of a simply supported, rectangular ortho- 
tropic plate [from Appendix III) is given by 

D Ha 4 u/3a^ + 2(D 12 + 2D g6 ) 8 4 u/3a^aa 2 + D 22 3 4 u/3a 4 

+ Xu + phu = q(a^,a 2 ,t) , (3.13) 


where, for the purposes of this study, the forcing function 
q(a^,a 2 ,t) will be defined by the equation 

q(a 1 ,a 2 ,t) = 




<5 ( a i _c i ) S ( a 2 ~c 2 ) + f. 




jut 


(3.14) 


(See Appendix III for the definitions of the various terms 
in this equation.) Thus, the acoustic (airborne) input is 
chosen to be a normally incident, spatially uniform, simple 
harmonic forcing function and the vibrational (structure- 
borne) input is chosen to be a simple harmonic, point vibra- 
tional load. 

The reasons for choosing a normally incident airborne 
input are threefold. As discussed in Appendix I, the pro- 
peller noise which impinges on the sidewall of an aircraft 
is an oblique incidence problem with a slowly varying 
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spatial pressure distribution caused by the extremely long 
wavelengths of the low frequency sound. Secondly, the case 
of normal incidence produces an airborne radiation effi- 
ciency that is numerically smaller than the corresponding 
oblique or random incidence cases. (This can be verified 
theoretically and will be recognized immediately if the 
reader is familiar with acoustic transmission loss theory. ) 
Thus, the normal incidence case produces a worst case 
scenerio in which the airborne radiation efficiency is as 
small as possible. (The reasons why this is desirable will 
become evident later.) Thirdly, it is comparatively much 
easier to construct an experimental apparatus that approx- 
imates the normal incidence condition. 

The point vibrational input was chosen to model the 
structureborne input primarily because it is easily mimicked 
experimentally using a shaker. It might be argued that a 
line load is a more realistic structureborne model since 
aircraft panels have ring frames and stringers attached to 
them. Because of the wide variation in design of composite 
airframes and in wing attachment, however, it is not clear 
that a line load model would be any more realistic. Fur- 
thermore, for the purposes of verifying the measurement 
method proposed in Chapter II, the point load model will, in 
principle, work just as well as a more complicated model. 
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D. SOLUTION TO THE FORCED RESPONSE DYNAMIC PROBLEM 


From Appendix III, the steady state solution for the 
dynamic response of an orthotropic plate due to combined 
simple harmonic airborne and structureborne inputs (the 

solution to equation (3.13)) is given by 
00 00 

u(a 1 ,a 2 ,t) = J J 
m = i n= i 


-JT e 


s e 
mn 


+ a 


-JT- 


mn 


•sin(mirct^/d^)sin(niTci 2 /d 2 ) 


j 

e 


y 


(3.15) 


where s m and a mrs are the structureborne and air- 
mn mn 

borne influence coefficents, respectively, of the vibration- 
al modes of the structure and are defined by the equations 


s mn = 4f s /( d i d 2 pha) mn^ sin(mirc 1 /d 1 )sin(mrc 2 /d 2 ) 


/ 


(1-mW) 2 + (2L wA) ) 2 

mn mn mn 


1 1/2 


(3.16) 


and 

a mn 

/ 

(See 


2 

= 4f a /(mnir phu mn ) 

(l-V.i.j* ♦ (25 

Appendix III for 


( 1-cos (mn) ) ( 1-cos (ntr) ) 


a )/<i) ) 

mn mn 


1/2 


details. ) 


(3.17) 
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Simple differentiation with respect to time of the 
equation for displacement yields the equation for the 
surface velocity of the structure (see Appendix III): 

v(a 1 ,a 2 ,t) = 



-^s 

-^a 1 

jw 

s e 
mn 

+ a mn e 


•sin (mira^/d^ ) sin (rnro^A^ ) 




(3.18) 


Equation (3.18) can then be used to calculate the space- time 

2 

averaged mean square surface velocity <v > . . From the 

r , t. 

results of the calculations given in Appendix IV: 

<v 2 > . = uj 2 /8 • 

r , u 


00 CO 


1 l 

m= i n= i 


mn 


+ a“_ + 2s a cos(<b -<1) ) 
— mn mn VT s T a 


mn 


(3.19) 


Redefining the terms in equation (3.19) as discussed in 
Appendix IV, the equation becomes 


2 2 2 2 

<v > . = <v > . + <v > . + <v > . , 

r,t sr,t ar,t sxa r,t ’ 


(3.20) 
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where 

2 

<v > r ^ = the space-time averaged mean square surface 
velocity, 

2 

<v > , = the structureborne component of the 

s r j l 

space- time averaged mean square surface velocity, 

2 

<v_> . = the airborne component of the space- 

9. L j L 

time averaged mean square surface velocity, 

2 

<v „ v >„ . = the cross term component of the space- 

S Xa XT j L. 

time averaged mean square surface velocity. 

See Appendices III and IV for the derivations of the 
components of equation (3.20) in terms of series expansions 
in the modal influence coefficients of the structure. The 
last term in equation (3.20) accounts for the interaction 
between the structureborne and airborne inputs in terms of 
the panel dynamics. Thus, equation (3.20) shows that the 
dynamic response of the panel to the combined airborne and 
structureborne inputs is, in general, not equal to the sum 
of the responses to the airborne and structureborne inputs 
individually. Similar cross term components arise in the 
derivations of the power input to the plate and the power 
dissipated by the plate. (See Appendix IV. ) One of the 
important features of equation (3.20) is that although it 
contains cross terms due to the combined inputs, it contains 
no cross terms between different modes of the structure. 

(See Appendix IV for more details.) (This absence of cross 
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terms between modes is a direct consequence of the ortho- 
gonality principle between modes of the structure.) Equa- 
tion (A4.12) indicates that when the structureborne and the 
airborne inputs are precisely 90 degrees out of phase, the 
inputs are uncorrelated and the cross term is zero. Thus, 
when the two inputs are uncorrelated, the space- time aver- 
aged mean square surface velocity of the panel due to the 
combined inputs is exactly equal to the sum of the space- 
time averaged mean square surface velocities due to the 
structureborne and airborne inputs acting individually. 

THEORY OF THE SOUND GENERATION OF PLATES 

Once the surface velocity distribution of the plate has 
been found (equation (3.18)), the sound generated by the 
plate can be calculated. The classical approach for cal- 
culating the sound radiation from a vibrating plate, which 
utilizes Huygens principle of superposition of simple 
sources, was used in this study to calculate the sound power 
radiated by the plate. This theory models each incremental 
area of the vibrating plate as a point monopole source near 
an infinitely rigid reflecting surface. (A review of the 
theory of point monopole sound radiation is presented in 
Appendix V. ) The mathematical details of this classical 
theory were derived for the particular case of combined 
structureborne and airborne inputs. The results of the 
derivations are presented in Appendix VI. The scheme of the 
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derivation is as follows: An equation is derived for the 
incremental acoustic pressure at a point in space due to the 
surface velocity of an incremental area of the plate. (See 
Appendix VI, eqn. (A6.3).) This equation is then simplified 
using a far-field approximation which assumes that the point 
of observation is at a distance from the plate which is 
much, much greater than the largest of the two plate dimen- 
sions. (See Appendix VI eqn. (A6.12).) This equation for 
the incremental pressure, as a function of the surface velo- 
city of an incremental area on the plate, is then integrated 
over the entire surface of the plate. (See Appendix VI, 
eqn. (A6.13).) The expression for the surface velocity 
(equation (A3. 16)) is then substituted into this equation 
and, after several complicated integrals are evaluated, an 
equation for the pressure at a point in the far acoustic 
field, as a function of polar and azimuthal angles, is ob- 
tained. (See Appendix VI eqn. (A6.36).) Since pressure and 
particle velocity are in phase in the far acoustic field, 
this equation for acoustic pressure can be used to calculate 
the time averaged acoustic intensity. (See Appendix VI, 
eqn. (A6.40).) The sound power radiated by the plate can 
then be calculated by integrating the equation for intensity 
over the polar and azimuthal angles. (See Appendix VI, eqn. 
(A6.43).) The final step in the analysis is to evaluate the 
expression for the special case of combined airborne and 
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structureborne inputs to the plate. (See Appendix VI, eqn. 
(A6.47).) This results in the equation on the following 
page for the time averaged sound power radiated by the 
plate : 

<TT(t)>. = 


2 it tt/2 


0 0 


(p o w 2 /(2tt) ) 2 sin(e)/(2p Q c o ) 


00 00 


A A L 

k= i l=i L 


•Si l\i 2 'V 2 


OO OO 

+ Re{ Y ]_ 




m= i 


00 

00 
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kl mn k 1 m n 
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d0d<J> 


(3.21) 
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(See Appendix VI for more information on the various terms 
in this equation.) Redefining the various terms in equation 
(3.21) as outlined in Appendix VI, the equation becomes 


<TT>. = <TT >. + <TT >. + <TT >. 

t s t at sxa t 


(3.22) 


where 

<TT> t = the total time averaged sound power radiated, 

<TT > = the time averaged sound power radiated due 

solely to the airborne input, 

<T!" s >£ = the time averaged sound power radiated due 

solely to the structureborne input, 

<TT :>. = the time averaged sound power radiated due 

to cross terms between the inputs. 


See Appendices V and VI for the derivations of the com- 
ponents of equation (3.22) in terms of series expansions in 
the modal influence coefficients of the structure. The last 
term in equation (3.22) accounts for the interaction between 
the structureborne and airborne inputs in terms of the sound 
power radiated. Thus, equation (3.22) similarly shows that 
the sound power generated by the panel due to the combined 
airborne and structureborne inputs is, in general, not equal 
to the sum of the individual sound powers radiated by the 
panel due to the airborne and structureborne inputs acting 
individually. Unlike equation (3.20) for the panel dyna- 
mics, however, equation (3.22) contains cross terms between 
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the different modes of the structure. (See Appendix VI for 
more details.) The additional cross terms between different 
flexural modeshapes arise because the spatial average of the 
acoustic intensity used in calculating the sound power takes 
place over a hemi-spher ical surface in the far acoustic 
field, and not over the surface of the plate (the area of 
spatial averaging for the analysis of the dynamic response 
of the plate). Thus, the orthogonality principle does not 
apply when calculating the sound power radiated by the 
plate. Furthermore, equation (A6.51) indicates that even 
when the structureborne and airborne inputs are precisely 90 
degrees out of phase (uncorrelated), there are still many 
non- zero cross term components between the two inputs which 
contribute to the total sound power radiated due to the 
existence of the cross terms between different flexural 
modes of the structure. Therefore, in general, the total 
sound power radiated due to the simultaneously combined 
structureborne and airborne inputs is never exactly equal to 
the sum of the sound powers radiated due to the structure- 
borne and airborne inputs acting individually. 

F. THEORY OF THE ACOUSTIC RADIATION EFFICIENCY OF PLATES 
The acoustic radiation efficiency, a, of a vibrating 
thin shell structure is defined as (see reference 48) 

o = <TT> t /Cp o c o <v 2 > rit S) , (3.23) 
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where 


<7T>^ = the time averaged sound power radiated by 

the structure, 

P q c o = the characteristic acoustic impedance of 
the fluid medium, 

2 

<v > . = the space- time averaged mean square surface 

r , t. 

velocity of the structure, 

S = the surface area of the structure. 

(See Appendix I, sections P - T for more information.) 

The equations for the acoustic radiation efficiencies 
for airborne and structureborne inputs to simple thin shell 
structures can be found in the literature. (See references 
48 and 49.) Theoretical discussions in the text of refer- 
ence 48 point out that in the frequency region below the 
critical frequency of a flat plate, the radiation efficien- 
cies of airborne noise and structureborne noise are differ- 
ent. (See reference 48 or Appendix I, section U for the 
definition of the critical frequency. ) Recently, Forssen 
and Crocker (reference 41) and McGary (references 39 and 40) 
have experimentally demonstrated this difference in the air- 
borne and structureborne radiation efficiencies on panels in 
the frequency range below the critical frequency using a two 
microphone cross spectral acoustic intensity probe. 

A rigorous explanation of the difference in the air- 
borne and structureborne radiation efficiencies for thin 
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shell structures can be made in terms of the vibrational 
modes of the structure. It is well known in classical 
vibration theory that the vibrational response of a thin 
shell structure to any type of input (airborne or structure- 
borne) can always be found from the superposition of the 
individual reponses of the flexural mode shapes of the 
structure. Mathematically, the vibrational response is 
expressed as a series expansion in the mode shapes of the 
structure, with each mode in the series being multiplied by 
a modal participation factor that is uniquely associated 
with that particular mode. The modal participation factors 
consist of both the magnitude of the response, called the 
influence coefficient, and the relative phase of the res- 
ponse for each of the individual modes. The influence 
coefficients determine the relative amount of control or 
influence that each individual mode retains over the total 
response of the structure. (See Appendix II for the eigen- 
value analysis and free vibrational response of a simply 
supported, rectangular orthotropic plate in terms of its 
natural modes. See Appendix III for the complete mathe- 
matical derivation of the forced vibrational response of a 
simply supported, rectangular orthotropic plate to simul- 
taneously combined airborne and structureborne inputs in 
terms of its natural modes.) 


28 



If the total vibrational response is controlled pri- 
marily by modes whose resonance frequencies are within the 
frequency band of excitation, the vibrational response is 
said to be resonance controlled. If the total vibrational 
response is controlled by modes whose resonance frequencies 
are below the frequency band of excitation, the vibrational 
response is said to be mass controlled. In the frequency 
bands at and above the fundamental resonance frequency of 
the structure, the vibrational response of the structure to 
most types of structureborne inputs (e.g. a point or line 
load) is primarily resonance controlled, and the vibrational 
response of the structure to acoustic inputs (which are for 
the most part distributed loads) can be either resonance or 
mass controlled. Below the critical frequency of a flat 
plate, a resonance controlled panel response is a very in- 
efficient noise generator (sometimes called acoustically 
slow) , whereas a non-resonant or mass controlled panel res- 
ponse is a very efficient noise generator (sometimes called 
acoustically fast). 

The difference in acoustic radiation efficiencies (in 
the frequency range below the critical frequency) for the 
resonance and non-resonant controlled cases stems from the 
difference in the effective flexural wave speeds for the two 
cases. For the case of a resonance controlled response, the 
product of the characteristic length (the wavelengths of the 
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resonant modes) and the forcing frequency is equal to the 
free vibrational or traveling flexural wavespeed in the 
structure. Below the critical frequency of the structure 
this free vibrational wavespeed is subsonic. Hence the 
resonance controlled response is an inefficient noise gener- 
ating mechanism (acoustically slow). For the case of a non- 
resonant or mass controlled response, the product of the 
characteristic length scale (the wavelengths of the non- 
resonant modes) and the forcing frequency are much larger 
than the free vibrational or traveling flexural wavespeed in 
the structure and can exceed the speed of sound in air 
(supersonic). Hence the non-resonant or mass controlled 
response can be a very efficient noise generating mechanism. 

At and above the critical frequency of the structure, 
both the non-resonant and resonant responses are efficient 
noise generators (acoustically fast). This is because 
structures are dispersive mediums for free vibrational 
(traveling) flexural waves. (The flexural wavespeed is 
proportional to the 1/2 power of the forcing frequency. ) 
Thus, at and above the critical frequency, the free vibra- 
tional wavespeed is equal to or exceeds the speed of sound. 
Hence, there is very little difference in the acoustic rad- 
iation efficiencies for structureborne and airborne inputs 
at these high frequencies. 
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G. DIAGNOSTIC APPROACH AND ANALYSIS 


In Section D, it was shown that the total space- time 
averaged mean square surface velocity of a simply supported, 
rectangular, orthotropic plate due to simultaneously com- 
bined airborne and structureborne inputs is given by (see 
Appendix IV, equation (A4.13)) 


2 2 2 2 

<v> =<v> . + <v > . + <v > , . 

r,t sr,t ar,t sxa r,t 


(3. 20) 


Similarly, in Section E, it was shown that the total 
time averaged sound power radiated by a simply supported, 
rectangular, orthotropic plate due to simultaneously 
combined structureborne and airborne inputs is given by (see 
Appendix VI, equation (A6.52)) 


<TT> . = <TT > , + <1 T > . + <TT > , 
t s t at sxa t 


(3.22) 


Multiplying equation (3.20) by the constant p o c Q 
S, and rearranging it algebraically one obtains the relation 


2 2 
pc <v> . S - p c <v" > . S 
M o o r,t M o o sxa r,t 


2 2 
=pc <v> . S + p c <v> .3 

oo sr,t oo ar,t 


(3.24) 


Utilizing the definition of radiation efficiency 
(equation (3.23)), and substituting into equation (3.24) the 
result is 
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(3.25) 


<TT> , /a - <TT 


> , /a 

sxa t' sxa 


= <TT > , /cr + 
s t' s 


<7T >. /a 

a t' a 


where 

a = the radiation efficiency for the simultaneously 

combined inputs, 

cr = the radiation efficiency of the cross terms 

S XcL 

between the two inputs, 

a s = the radiation efficiency of the structureborne 

component, 

a a = the radiation efficiency of the airborne 

component. 

Now rearranging equations (3.22) algebraically, the 
following system of equations is obtained: 


<TT>./cr - <TT >./a = <7T >, /a + 

t / sxa t' sxa s t s 


<TT >./a 
a t' a 


(3.25) 


<1T> . 


<TT >. = <TT >. + <TT >. 

sxa t st at 


(3.26) 


Solving this system of two equations and two unknowns for 
structureborne and airborne components of the sound power, 
one obtains 


<T Vt = <Tr> t (a a /a) (a - a s )/(a a " a s :i 


<T '"sxa > t ^°a^°sxa^ ^ a sxa " a s^ // ^ a a - a s^ ’ (3.27) 
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(3.28) 


<T Vt = <TT> t (G s /a) (a a " a)/[a a “ a s ) 


<TT >. (a /a ) (a -a )/(a - a ) . 

sxa t s sxa a sxa a s 


Equations (3.27) and (3.28) are exact expressions which 
give the airborne and structureborne sound power components 
respectively in terms of the acoustic radiation efficien- 
cies, the total sound power, and the cross term sound power. 
As discussed earlier, the cross term sound power depends on 
the relative phase between the two inputs. This dependence 
on phase makes the cross term sound power, and cross term 
radiation efficiency difficult to measure in practice. In 
contrast to the cross terms, the structureborne and airborne 
radiation efficiencies do not depend on the magnitudes or 
the relative phase of the inputs, and (as will be discussed 
in a later section) can easily be measured. Similarly, the 
total sound power and radiation efficiency for the combined 
inputs can be measured with relative ease. For these rea- 
sons, the following definitions of estimates of the struc- 
tureborne and airborne components of the sound power are 
made : 

An estimate of the airborne sound power is 


<T r a > t = <Tr> t ^ a a^ a ^ " °s- )/ ^ a a “ °s-* ’ (3.29) 

so that 
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(3.30) 


<T V = <T ^"a > t 


+ <TT >, (a /a ) (a - a )/(a - a ) . 

sxa t a' sxa sxa s a s 


An estimate of the structureborne sound power is 


< TT s>t = <TT> t (a g /a) (a a - a)/(a a - a g ) , 


(3.31) 


so that 


<TT >, = <TT >. 
st s t 


+ <TT >. (a /a ) (a - a )/(a - a ) 

sxa t s 7 sxa a sxa 7/ v a s 


(3.32) 


Note that these estimates of the structureborne and 
airborne components of the sound power involve only the 
airborne and structureborne radiation efficiencies, and the 
total sound power and radiation efficiency due to the 
combined inputs. Another important feature of these 
estimates is that their sum is power preserving, i.e. 

<TT>. = <T T >. + <TT >, . (3.33) 

Cr s t a. t 

Furthermore, equations (3.29) and (3.31) suggest that 
measurable differences in radiation efficiencies for air- 
borne and structureborne noise permits the separation, and 
prediction of an estimate of the airborne and structureborne 
components of the total sound power radiated by some unknown 
combination of structureborne and airborne inputs to a 
panel . 
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Hj_ CONTRIBUTION OF THE CROSS TERM TO THE DIAGNOSTICS 

Recall equations (3.30) and (3.32) from the previous 
section of this chapter: 


<T V = <T Vt 


<"TT >. (a /a ) (a - c )/(a - a ) 

sxa t v a sxa sxa s ^ a s 


(3.30) 


<T Vt * <T Vt 


+ <TT >. (a /a ) (a -a )/(a - a ) 

sxa t v s / sxa' v a sxa' ^ a s' 


(3.32) 


An analysis was performed on these two equations to 
determine how much of the cross term sound power, 

<TT >. , is attributed to the estimated airborne and 
structureborne sound power components. The results of the 
analysis is given in figures (III-2) and (III-3). Figure 
(III-2) shows the fractional amount of the cross term sound 
power (0 to 1) which is attributed to the estimate of the 
structureborne component plotted against the fractional rad- 
iation efficiency of the cross term. Figure (III-3) shows 
the fractional amount of the cross term sound power (0 to 1) 
which is attributed to the estimate of the airborne compo- 
nent plotted against the fractional radiation efficiency of 
the cross term. The fractional radiation efficiency of the 


cross term is defined by 

a c = (a - a )/(a - a), 

f x v sxa s J / a s J 


(3.34) 
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Thus, <j£ x is equal to 0 when the cross term radiation 

efficiency, a , is equal to the structureborne radiation 
sxa. 

efficiency, a , and is equal to 1 when the cross term 
s 

radiation efficiency, a , is equal to the airborne rad- 

SX3. 

iation efficiency, a . The various curves in figures 
(III-2) and (III-3) show how the results of the analysis 
depend on the difference between the structureborne and air- 
borne radiation efficiencies. If the difference between the 

airborne and structureborne efficiencies, a - c , 

a s 

is small (e.g. 5 dB or less), the relationship between the 
fraction of the cross term attributed to the estimates and 
the fractional cross term radiation efficiency is nearly 
linear. If the difference between the structureborne and 
airborne efficiencies, a - a , is large (e.g. 20 or 
25 dB) the relationship between the fraction of the cross 
term attributed to the estimates and the fractional cross 
term efficiency is seen to be hyperbolic with most of the 
cross term sound power being attributed to the estimate of 
the airborne component irregardless of the radiation effi- 
ciency of the cross term. Thus, if the measurements show a 
large difference in the radiation efficiencies of the indi- 
vidual components, one can be reasonably certain that the 
effects of the cross terms will manifest themselves pri- 
marily in the estimates of the airborne component of the 
total sound power. 
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I. ERROR ANALYSIS OF THE PROPOSED DIAGNOSTIC METHOD 


Recall 

<TT >, = <TT> 

a t 

<TT >, = <TT> 
s t 


equations (3.29) and 
t " a )/(a 

t (a g /a) (ff & - a)/(a 


a 


a 


(3.31): 





(3.29) 

(3.31) 


An error analysis was performed to determine the 
severity of the errors in these two estimates due to errors 
in the measurements. The results of this error analysis are 
given in Appendix VIII. The analysis includes the effects 
on the estimates due to measurement errors in the structure- 
borne radiation efficiency, a g , the airborne radiation 
efficiency, a . and the radiation efficiency of the 

a. 

combined noise, a. The severity of error is found to be 

sensitive to the actual values of a and (a -a ) in 

v a s 

some cases. The main conclusion of the study is that 
measurement errors are much more likely to cause inaccur- 
acies in <TT >. than <TT >, . 

s t at 

J^_ IMPLEMENTATION OF THE ANALYTICAL MODEL 

Two computer programs were written in ANSI Standard 
FORTRAN 77 to implement the analysis outlined in the pre- 
vious sections. The first program, entitled MODAL, cal- 
culates the natural frequencies of the plate under invest- 
igation. The second program, entitled SOUND, uses the 
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output of the first program and calculates the dynamic res- 
ponse, sound power radiated, . and radiation efficiency as a 
function of forcing frequency for a specified combination of 
the airborne and structureborne inputs discussed earlier. 
Source code listings of these two programs along with de- 
tailed discussions of their calculation schemes are given in 
Appendix IX. 

Compilation and execution of these two programs were 
performed on the Analysis and Computation Division's Cyber 
175 computers (built by Control Data Corporation) under the 
Network Operating System (NOS). The ASCII data files con- 
taining the results of these computations (entitled DATA 
files) were then down loaded to a Tandy 2000 personal com- 
puter using a public domain terminal emulation program 
(COMSH ) , loaded into a spreadsheet-graphics program (LOTUS 
1-2-3), and plotted on peripheral devices such as a Radio 
Shack CGP-220 color ink jet printer or a Hewlett-Packard 
7470A pen plotter. 

Several preliminary runs of programs MODAL and SOUND 
were performed in order to ensure that the computations were 
free from serious error. The results of several simple 
tests designed to validate the computer codes are presented 
in Appendix X. 
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K. SIMULATION OF THE PROPOSED DIAGNOSTIC METHOD 

A third computer program, entitled PREDICT, was dev- 
eloped on the Tandy 2000 system to be used in conjunction 
with three of the output data files (DATA) generated by the 
program SOUND. The purpose of program PREDICT is to simu- 
late the proposed diagnostic method. PREDICT uses the rad- 
iation efficiencies calculated by program SOUND for the 
cases of purely airborne and purely structureborne inputs 
(contained in two different DATA files) to separate and 
predict the airborne and structureborne sound power compo- 
nents of some combined inputs case of interest (contained in 
a third DATA file). Thus, program PREDICT is the implement- 
ation of equations (3.29) and (3.31). A listing of program 
PREDICT is given in Appendix IX. 
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Figure (III-2) 


Fraction of the cross term sound power attributed to the estimate of the 
structurebor ne sound power. 


CONTRIBUTION TO AIRBORNE ESTIMATE 



Figure (III-3) — Fraction of the cross term sound power attributed to the estimate of the 

airborne sound power. 


Chapter IV 


THEORY OF MEASUREMENT AND EXPERIMENTAL APPARATUS 
THEORY OF MEASUREMENT 

In chapter III a diagnostic analysis was suggested 
whereby the structureborne and airborne components of the 
total sound power radiated by a plate or thin shell struc- 
ture could be predicted. To implement this diagnostic 
analysis, the airborne and structureborne radiation effi- 
ciencies must be known and the sound power and radiation 
efficiency of the structure for the combined inputs case of 
interest must be measured. The present section deals with 
how the sound power and radiation efficiency of a structure 
can be measured in theory. 

Recall that the time averaged sound power radiated by 
the surface of a structure is given by the equation (see 
Appendix I, equations (Al.50) and (Al.55)) 


<TT> 


t 


<I*n>^ dS 


<X ’ n> r , t S » 


(4.1) 


where I is the acoustic intensity vector measured at the 
surface, n is the unit normal vector to the surface, and S 
is the surface area. The acoustic intensity vector normal 
to the surface can be measured using a two microphone, cross 
spectral, acoustic intensity probe. The theoretical basis 
for this type of probe is summarized in Appendix XI. The 
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acoustic intensity is calculated from the imaginary part of 
the one-sided cross spectral density between the two signals 
produced by two closely spaced microphones. The equation 
for the time averaged acoustic intensity is given by (see 
Appendix XI, equation (All. 22)) 

<I> t = -Q 12 Af/(p o wAx) , (4.2) 

where 

Q ^2 = minus the imaginary part of the one-sided cross 

spectral density between the two microphone signals, 

Af = the frequency resolution (bandwidth) in Hertz, 

P Q = the density of the acoustic fluid medium, 

to = the radian frequency of the acoustic disturbance, 

Ax = the spacing between the two microphones. 

The quadrature spectral density, can be 

easily measured with a dual channel or a multichannel Fast 

Fourier Transform (FFT) analyzer. The remaining terms in 

equation (4.2) are either constant or are parameters (such 

as bandwidth) which are set by the FFT analyzer. The space- 

time averaged acoustic intensity, <I> . , is found in prac- 

r , t. 

tice by slowly sweeping the two microphone acoustic inten- 
sity probe near the surface of the intended measurement area 
as the FFT analyzer calculates the time averaged cross spec- 
tral density between the microphone signals. Thus, the 
space averaging and the time averaging of the cross spectral 
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density is performed simultaneously. (The consequences of 
the use of this method of space- time average for obtaining 
the space-time averaged acoustic intensity is discussed 
briefly in Appendix XI.) Once the space-time averaged cross 
spectral density between microphone signals is measured, the 
space-time averaged acoustic intensity is calculated as a 
function of frequency using equation (4.2) by computer, or 
by the FFT analyzer (depending on the analyzer's level of 
sophistication). The total sound power radiated is then 
calculated by multiplying the space-time averaged intensity 
by the measurement area. (See equation (4.1).) The 
measurement area may be the surface area of the structure, 
providing that the intensity probe is swept close enough to 
the surface, and if there are no stiffeners, etc., attached 
to the intended measurement surface. 

Recall from Appendix I, equation (Al.53) that the 
space-time averaged mean square surface velocity of a 
vibrating surface is given by the relation 


<v^> . = 1/S 

r, t 


R vv (0) dS . < 


G (f) df> , 
w v J r ’ 


(4.3) 


where G vv is the one-sided auto spectral density of the 
surface velocity. Substituting an estimate of the auto 
spectral density from the equations given in Appendix I (see 
equation (Al.38)) 
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oo 


(4.4) 


<v 2 > 


r , t 


= < 


2/T V (f)V(f) df > 


Substituting from the relationships between velocity and 
acceleration given in Appendix I (see equation (Al.9)) 


<v 2 > 


r, t 


2/T 


A*(f)A(f)/w 2 


df > 


(4.5) 


And from the definition of the estimated value of the one- 
sided auto spectral density of acceleration (see equation 
(Al. 38)) 


<v 2 > 


r, t 


= < 


G aa ( f )/w' 


df > 


(4.6) 


where G is the one-sided auto spectral density of the 

3 . 9 . 

surface acceleration. 

The auto spectral density of the surface acceleration 
can be easily measured using an accelerometer together with 
an FFT analyzer. The time averaged mean square surface 
velocity can then be calculated from equation (4.6) by using 
a computer, or by integrating directly on the FFT analyzer 
(if that function is available on the analyzer). The space- 
time averaged mean square surface velocity is obtained in 
practice by measuring and calculating the time averaged mean 
square surface velocity from each of several accelerometers 
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which are attached to the measurement surface, and then 
simply averaging the results together as follows: 


<v 2 > 


r, t 


N 

1/N } (<v 2 > 

k= i 


t-’k 


(4.7) 


This space averaged estimate will be normally distributed 
about its true value as the number of samples (acceler- 
ometers), N, increases to infinity, irregardless of the 
particular spatial distribution of the surface velocity. 
(This is a direct consequence of the central limit theorem. ) 
One must be careful, however, not to place so many acceler- 
ometers on the surface that the added mass produces a large 
bias error in the measurements. Furthermore, the acceler- 
ometers should be placed on the surface in a random pattern 
(irregular as opposed to a regular pattern) so that the 
space average is not biased toward any particular mode shape 
of the structure. 

Once the time averaged sound power and the space- time 
averaged mean square surface velocity are calculated, the 
acoustic radiation efficiency of the structure may be cal- 
culated using equation (3.23): 

a = <Tr> t /(P 0 c 0 <v2> r t s 3 • (3.23) 
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B. MEASUREMENT PROCEDURE 


In the last section, the measurement theory used in 
determining the sound power and radiation efficiency of a 
structure was discussed. The present section outlines the 
implementation of this measurement theory in terms of a 
procedure for the measurement of sound power and radiation 
efficiency. 

The following procedure is proposed for the measurement 
of the sound power and radiation efficiency of an exposed 
surface in the interior of a propeller driven aircraft: 

Collect the Necessary Equipment - The following items 
are necessary for the measurements: 

1. A two microphone acoustic intensity probe. Homemade 
probes will do just as well as the commercially available 
varieties. For measurements over the 0-1000 Hz frequency 
range, it is recommended that the probe consist of two 1/2 
inch high gain microphones (to lower the noise floor and 
increase the signal level). 

2. A minimum of two miniature (2 gram) piezoelectric 
accelerometers with the appropriate charge amplifiers (if 
charge amplifiers are necessary). 

3. A dual channel or multichannel FFT analyzer with its own 
storage medium. The analyzer must be capable of being 
interfaced with a minicomputer or microcomputer via IEEE- 
488, RS-232C, or some other interface communication system. 
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4. Acoustic absorptive material to place in the enclosed 
space in which the measurements are to take place. 
Polyurethane foam or fiberglass wedges are recommended. 

5. A barometer and thermometer. 

6. A pistonphone and vibration calibration device. 

Calibrate the Equipment - The following steps should be 
performed in calibrating the equipment: 

1. The microphones used for the acoustic intensity probe 
should be calibrated with the pistonphone to obtain the 
appropriate gain factors. If the acoustic intensity probe 
is home made, it may also be necessary to phase calibrate 
the microphones. The need for phase calibration is depend- 
ent on the frequency range over which the measurements are 
intended. (See Appendix XI for a discussion of the possible 
measurement errors associated with the two microphone method 
of acoustic intensity measurement. ) 

2. Select the microphone spacing for the acoustic intensity 
probe. The exact spacing is a function of the frequency 
range of interest of the intended measurement. (See 
Appendix XI for a brief discussion of the selection of an 
appropriate microphone spacing.) 

3. Calibrate the accelerometers for the appropriate gain 
factors using the vibration calibration device. 
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4. Record all of the gain factors, the ambient temperature, 
and the barometric pressure for later computer analysis and 
data reduction. (The computer software uses the atmospheric 
conditions to calculate the characteristic acoustic imped- 
ance of the fluid medium, p c . which is used in the 

o o 

acoustic intensity calculations.) 

Prepare the Structure for the Measurements - Select the. 
sidewall area for the intended measurements as follows: 

1. The area selected may be any area of the structure which 
is basically of a thin plate or shell type construction. 

The boundary conditions on the measurement area are unim- 
portant. Any type of boundary condition is acceptable. The 
chosen area may include stiffeners (e.g. ring frames or 
stringers), however it is best to choose the area small 
enough so that the majority of the stiffeners are along the 
boundaries. It is important that the trim material is re- 
moved from the measurement area before the measurements are 
performed. 

2. Place the acoustic absorptive material inside the re- 
ceiving space so as to minimize reflected sound. It may also 
be necessary to erect an acoustic barrier in the receiving 
space to prevent errors in the measurements due to flanking 
if there is a strong noise source near the selected measure- 
ment area. 
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Procedure - Use the following procedure: 

1. Attach the array of accelerometers to the measurement 
surface in a random pattern. (A regular pattern may bias 
the measurements toward one particular mode shape or a 
family of mode shapes.) 

2. Turn on the noise sources of interest and adjust the 
gain factors so that the resulting vibrational and acoustic 
signal levels are above the instrumentation thresholds over 
the frequency range of interest. 

3. Measure the time averaged auto spectral density of 
acceleration of the measurement surface using the array 
accelerometers and the multichannel FFT analyzer. Store 
this information on the analyzer. 

4. Use the two microphone acoustic intensity probe along 
with the FFT analyzer to measure the space-time averaged 
cross spectral density between the microphones by slowly 
sweeping the handheld probe over the surface. Store this 
information on the analyzer. 

5. Turn off the noise sources of interest. 

6. Use the information gathered in steps 3 and 4 to cal- 
culate the acoustic radiation efficiency of the measurement 
surface as outlined in section A. 
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C. PROPOSED DIAGNOSTIC PROCEDURE 


The following diagnostic procedure is proposed for the 
separation and prediction of the structureborne and airborne 
components of the total sound power radiated by an exposed 
surface in the interior of a propeller driven aircraft: 

Collect the Necessary Equipment - In addition to the 
equipment discussed in section B, the following additional 
items are necessary for the measurements: 

1. A mechanical shaker with its associated power amplifier. 

2. A loud speaker with its associated power amplifier. 

3. A white noise (broad band frequency) signal generator. 

Determine the Structureborne Radiation Efficiency - 

1. Attach the mechanical shaker to the same attachment 
point as the source of the structureborne noise. (In the 
case of an aircraft, this would be the points at which the 
engines or wings attach to the aircraft.) If the structure- 
borne noise source has more than one attachment point, then 
attach the shaker in line with the center of mass of the 
source. Keep in mind that it is not important that the 
shaker reproduce the level or the character of the actual 
structureborne noise source, but only that the vibrational 
energy takes the same path as the vibrational energy pro- 
duced by the actual structureborne noise source. 

2. Connect the white noise generator to the mechanical 
shaker's power amplifier. 
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3. Turn on the noise generator and adjust the level so that 
the resulting vibrational and acoustic signal levels are 
above the instrumentation thresholds over the frequency 
range of interest. 

4. Perform the measurment procedure outlined in section B. 

5. Turn off the noise generator, and disconnect the mech- 
anical shaker. Leave the acoustic absorptive material in 
place . 

6. Use the information gathered in step 4 to calculate the 
acoustic radiation efficiency of the measurement surface for 
the case of a purely structureborne input. 

Determine the Airborne Radiation Efficiency - 

1. Position the loud speaker so that it is approximately 
the same distance from the measurement area as the source of 
the airborne noise. (in the case of an aircraft, this may 
be in or aft of the propeller tip path plane. ) Orient the 
loud speaker so that it has approximately the same direct- 
ivity with respect to the measurement area as the actual 
airborne noise source. Keep in mind that it is not import- 
ant that the loudspeaker reproduce the level or character of 
the actual airborne noise source, but only that the path of 
the acoustic propagation is the same. 

2. Connect the white noise generator to the loudspeaker to 
the loudspeaker's amplifier. 
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3. Turn on the noise generator and adjust the level so that 
the resulting vibrational and acoustic signal levels are 
above the instrumentation thresholds over the frequency 
range of interest. 

4. Perform the measurment procedure outlined in section B. 

5. Turn off the noise generator, and remove the loud- 
speaker. Leave the acoustic absorptive material in place. 

6. Use the information gathered in step 4 to calculate the 
acoustic radiation efficiency of the measurement surface for 
the case of a purely airborne input. 

Measurements on the Noise Sources of Interest - 

1. Turn on the actual noise source of interest. (in the 
case of an aircraft, bring the engines to the desired power 
setting and feather the propellers to the desired thrust 
setting. Altitude and airspeed may also be variables of 
interest. ) 

2. Perform the measurment procedure outlined in section B. 

3. Use the information gathered in step 2 to calculate the 
sound power and the acoustic radiation efficiency of the 
measurement surface for the combined inputs case of 
interest . 

Calculations - Now use equations (3.29) and (3.31) to 
separate and predict the structureborne and airborne com- 
ponents of the total sound power radiated by the measurement 
surface for the combined inputs case of interest. 
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D. EXPERIMENTAL APPARATUS 


An experimental study of the proposed diagnostic method 
was performed using the NASA Langley Research Center's 
acoustic transmission loss apparatus. This facility is a 
hard walled, two room facility designed for acoustic trans- 
mission loss measurements using the classical room acoustics 
method. The two rooms have an adjoining wall which is de- 
signed so that simple or built-up aircraft panels can be 
mounted between the two rooms. With this arrangement, test 
panels could be subjected to the desired acoustic and vibra- 
tional inputs in the source room while the surface velocity 
and sound power radiated by the panels could be measured in 
the receiving room. Since the analytical model discussed in 
chapter III assumes that the test panels radiate sound to an 
acoustic free field condition, and since the accuracy of 
intensity measurements are in question under reverberent 
conditions, the receiving room of the transmission loss 
apparatus was modified to semi-anechoic conditions. This 
was accomplished by covering the back wall of the receiving 
room with .91 m deep acoustic wedges and covering the floor 
of the room with .46 m deep acoustic wedges. The acoustic 
wedges, constructed of polyurethane foam, are shown in the 
photograph of figure (IV-1). No further modifications of 
the transmission loss apparatus were required to perform the 
measurements. (Additional information regarding the 
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acoustic properties of the NASA acoustic transmission loss 
apparatus is available in references 50 and 51.) 

A special apparatus for mounting the test panels was 
constructed so that the experimental conditions would 
emulate the conditions assumed for the analytical modeling. 
The apparatus consisted of a speaker box which completely 
enclosed the incident side of the test panels. Two small 
holes were drilled through the enclosure so that a steel rod 
and shaker could be attached to the test panels in one of 
two locations. The test panels were then clamped in the 
mounting brackets seen in figure -(IV-2) in front of six 4 
inch diameter loudspeakers . The mounting brackets, shown 
in close-up in figure (IV-3), were constructed with a rubber 
O-ring type material so that the test panels would have some 
rotational degree of freedom, thus approximating the simply 
supported conditions assumed in the analytical modeling. 

The array of six loudspeakers shown in figure (IV-2) 
were used to produce a normally incident, spatially uniform 
acoustic (airborne) input to the test panels. The loud- 
speakers were positioned 5.7 cm from the surface of the test 
panels, thus insuring that the acoustic resonances in the 
cavity between the panels and the speakers have natural 
frequencies much greater than 1000 Hz. The elimination of 
any significant influence due to the cavity modes helped to 
produce an acoustic input that was nearly uniform over the 
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0-1000 Hz frequency range. The small distance between the 
speakers and the panels also insured that the direct sound 
field from the speakers would overwhelm the effects of any 
cross modes in the cavity, thereby approximating the spat- 
ially uniform conditions. A preliminary set of measurements 
were performed on the loudspeakers to insure that they were 
in phase and produced the same level of sound over the 0- 
1000 Hz range. A plot of the measured space- time averaged 
acoustic intensity radiated by each of the six loudspeakers 
to the free field over the 0-1000 Hz frequency range is 
given in figure (IV-4). 

A 44.5 N (force) vibration shaker was used in conjunc- 
tion with a .635 cm diameter steel rod to simulate the point 
vibrational (structureborne) input. The shaker was mounted 
outside of the speaker box by suspending it freely with 
bungee cord that was attached to a scaffold- type apparatus. 
This arrangement ensured that the shaker-rod- panel system 
had a low natural frequency and reduced any d.c. component 
of the point forcing function to a minimum. The threaded 
rod was attached to the panel, in each case, by drilling a 
hole in the panel, slipping the rod through the hole, and 
tightening a hex nut down on each side of the panel. 

A typical example of the forcing function produced by 
this apparatus over the 0-1000 Hz frequency range is given 
in the plot of figure (IV-5). 
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E. TEST PANELS 


The test structures chosen for use in this study con- 
sisted of rectangular, simply supported plates constructed 
of isotropic or composite materials. Six different types of 
test panels were chosen for study by the author. The NASA 
Langley Research Center's Materials Division then fabricated 
the test panels to the desired specifications. A summary of 
the physical characteristics and stiffness properties of the 
plates chosen for study are given in Tables IV-1, IV-2, and 
IV- 3. 

E\_ INSTRUMENTATION AND DATA ACQUISITION 

A block diagram of the instrumentation used for the 
measurements is shown in figure (IV-6). The specifications 
for the instruments are given in Appendix XII. 

The white noise generator, shown in the figure, provid- 
ed a broadband random signal (0-5000 Hz) that was used to 
simultaneously drive both the loudspeakers and the shaker 
system. This single source ensured that the airborne and 
structureborne inputs were fully coherent. The signal was 
filtered using both a high pass and a low pass filter so 
that the sound radiated by the panels would be concentrated 
in the 100-1000 Hz frequency range. A signal attenuator was 
used to adjust the level of the airborne input so that the 
relative amounts of airborne and structureborne noise radi- 
ated by the panels were roughly equal. 
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An array of four miniature (2 gram) piezoelectric 
accelerometers were used to measure the space-time averaged 
surface acceleration of the panels. The accelerometers were 
attached to the surface of the panels using double sticky 
back tape in a random pattern. Each accelerometer was used 
in four different locations for a given measurement of the 
surface acceleration. Thus, a total of sixteen different 
accelerometer locations were used to obtain the space aver- 
age for each measurement of the surface acceleration. A 
total of 100 ensemble averages at each accelerometer loca- 
tion (obtained by the FFT analyzer) were used to obtain the 
time average in each case. 

The two microphone acoustic intensity probe, used to 
measure the sound power radiated by the panels, is shown in 
figure (IV-7). The probe consisted of two 1.27 cm diameter 
high gain microphones in a face-to-face configuration. The 
solid nylon cylindrical spacer between the microphones pro- 
vided a constant separation distance of 50 mm. This separa- 
tion distance between microphones ensures that the sound 
power measurements are accurate over the 100-1000 Hz freq- 
uency range. Below 100 Hz the sound power measurements are 
suspect due to phase mismatch errors. Above 1000 Hz the 
sound power measurements are inaccurate due to finite dif- 
ference error. (See Appendix XI for details. ) Since the 
microphone interchange technique was used for the 
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measurements, the intensity probe was used twice for any 
given measurement of space-time averaged acoustic intensity. 
(This method reduces the phase mismatch error. See Appendix 
XI for details.) The FFT analyzer obtained 200 ensemble 
averages for each of the two passes of the intensity probe. 
Thus, a total of 400 ensemble averages were used to obtain 
the time averaged intensity. The space average was obtained 
by slowly sweeping the intensity probe near the surface of 
the panel (as outlined in section A). 

Calibration of the accelerometers and the microphones 
were performed using a GENRAD 1557A vibration exciter and a 
GENRAD 1986 sound level calibrator respectively. The trans- 
ducers were calibrated prior to the measurement on each day 
that a measurement was to take place. The voltages produced 
by the transducers due to a known input were read on the 
digital voltmeter shown in figure (IV-6). Calibration was 
then implemented by adjusting the gain on the accelerometer 
signal conditioners and the microphone amplifiers so that 
the gain factors were always .01 volt/(m/sec ) for the 
accelerometers and .1 volts/Pascal for the microphones. The 
oscilloscope shown in figure (IV-6) was used to insure that 
the signals received during calibration were free of distor- 
tion, thus making certain that the transducers were in 
acceptable operating condition. The scanner shown in figure 
(IV-6) was utilized during both calibration and measurement 
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so that any of the 8 data channels could be monitored with- 
out disconnecting or reconnecting any wire leads. 

The data acquisition was performed by the GENRAD 2515 
8-channel FFT analyzer. This data acquisition system is 
particularly well suited for work on aircraft because of its 
small size and portability. Two of the 2515 's 8 channels 
were used for the microphones, while four of the channels 
were used for the accelerometers. The two remaining 
channels were used to monitor the white noise input signal 
and the signal provided by the force gauge (see figure (IV- 
6). In addition to the data obtained by the analyzer, the 
atmospheric conditions including temperature and barometric 
pressure were recorded each day and entered into the data 
files. The atmospheric data were used later in the computa- 
tions to calculate the density, p , and characteristic 

o 

acoustic impedance, p q c o , of the fluid medium. (These 
quantities are used in calculating the sound power radiated 
by the plate.) 

G. DATA REDUCTION AND ANALYSIS 

Once the data were obtained, they were uploaded to the 
Acoustics Division's VAX 11/780 via the KERMIT file transfer 
protocol. A computer program, written in FORTRAN 77, was 
then used to reduce the acquired data and calculate the 
sound power radiated and radiation efficiency of the panels 
as a function of frequency. The program was written so that 
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the ASCII data files produced by the experiments had pre- 
cisely the same file format as the DATA files produced by 
the analytical program SOUND. These data files produced by 
the experiments were then down loaded to a Tandy 2000 per- 
sonal computer via the KERMIT file transfer protocol. The 
data files were then loaded into a spreadsheet-graphics 
program (LOTUS 1-2-3) and plotted on either a Radio Shack 
CGP-220 ink jet printer or a Hewlett-Packard pen plotter. A 
block diagram of the equipment used in this research for 
data acquisition, reduction, and analysis is shown in figure 
(IV-8 ) . 

H_;_ VERIFICATION OF THE PROPOSED DIAGNOSTIC METHOD 

Program PREDICT, discussed in section K of chapter III, 
was used in conjunction with the data files produced by the 
experiments to test the diagnostic method proposed in 
chapter III. PREDICT uses the measured radiation efficien- 
cies for the cases of a purely airborne and a purely struc- 
tureborne inputs (contained in two different data files) to 
separate and predict the airborne and structureborne sound 
power components of some combined inputs case of interest 
(contained in a third data file). 
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Table II-l -- Description of the plate construction. 


Plate 

No. 

Material 

Number 
of plies 

Fiber 

lay-up 

Density 

(kg/m 2 ) 

Thickness 

(mm) 

1 

Aluminum 

1 

isotropic 

2.22 

0.79 

2 

t 

i . . .. _ . 

* , . 
Graphite 



8 

0/90 

/0/90 

1.59 

1.02 

! 

3 

* , . 
Graphite 

8 

+45/-45 

/+45/-45 

1.59 

1.02 

i 4 

i 

i 

* . 1 
Graphite j 

. . . _ _ ... 

1 16 

! 

+45/-45 
/+45/-45 
/+4S/-45 
/+45/- 45 

3.05 

! 

1.85 

5 

* 

Kevlar 

8 


1.37 

1.02 

6 

1 * 

Fiberglass 

i 

8 

+45/-45 

/+45/-45 

2.21 

1.02 


Table II-2 -- Material properties of the plate plies. 


Plate 

No. 

Material 

E 11 

(Pa) 

E 22 

(Pa) 

G 

(Pa) 

u 12 

1 

1 

Aluminum 

. 731ell 

. 731ell 

. 276ell 

. 330 

2-4 | 

Graphite 

. 137el2 

. 965el0 

. 480el0 

.300 


* 

Kevlar 

. 760ell 

. 550el0 

. 210el0 

.340 

6 

★ 

Fiberglass 

. 390ell 

. 900el0 

. 240el0 

mmm 


See note on page 64. 
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Table II-3 -- Calculated bending rigidities of the plates. 


Plate 

* — 1 
i — 1 

D 

°22 

°12 

°66 

°16 

°26 


(N*m) 



(N*m) 

(N*m) 

(N*m) 

(N*m) 

(N*m) 

| 


3.42 

1.13 

1.15 

0.00 

0.00 

2 


4.40 

0.26 

0.42 

0.00 

0.00 

3 

3.82 

3.82 

2.97 

3.13 

|~ ~ 1 

1.06 

1.06 

4 

22.77 

22.77 

17.70 

18.70 

3.17 

3.17 

5 

2.09 

2.09 

1.72 

1.73 

0.59 

0.59 

6 

1.42 

1.42 

0.99 

0.96 

i 

0.25 

0.25 I 

1 

1 i 


* 

NOTE - All materials used for construction of the composite 
panels were procured as pre-preg tape from: 

FIBERITE corporation 
General Headquarters 
501 West Third Street 
Winona, Minnesota 55987 


The material designations were as follows: 
Graphite tape - Celion 6000/Epoxy resin 934 
Fiberglass tape - E glass/Epoxy resin 934 
Kevlar tape - Kevlar 49/Epoxy resin 934 


The panels were layed-up and cured by NASA Langley's 
Materials Division. 
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Figure (IV-4) -- Acoustic intensity radiated by each of the six loudspeakers 
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Figure (IV-5) 


Forcing function produced by the shaker- rod-panel system 
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Figure (IV-6) 


Block diagram of the instrumentation used for the measurements 


















Figure (IV- 7) -- Two microphone acoustic intensity probe used for the measurements. 
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Data acquisition, reduction, and analysis systems 







Chapter V 

RESULTS AND DISCUSSION 


A^ INTRODUCTION 

This chapter contains selected results of the studies. 
The analytical results shown here are intended to simulate 
the performance of the proposed diagnostic method under 
various operating conditions. The experimental results are 
then presented alongside the analytical results and are 
intended to verify the predicted trends. Since the elec- 
tronic equipment necessary for controlling the relative mag- 
nitude and phase of the inputs as a function of frequency 
was unavailable, the corresponding analytical and experi- 
mental cases presented are not exactly the same. Therefore, 
the direct comparison of the analytical and experimental 
results (overlaying the plots) would be ludicrous. The 
reader should keep in mind that the purpose of the analyti- 
cal and experimental comparisons is not to compare quanta- 
tive results, but rather is intended to show that predicted 
and observed trends are similar in a qualitative sense. 

(The electronic boards necessary for controlling the phase 
and amplitudes of the inputs as a function of frequency are 
being constructed at this writing. Future tests are planned 
for an aircraft fuselage using these gain/phase boards.) 

The analytical results presented in this chapter were 
produced using the theoretical developments and computer 
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programs (SOUND and PREDICT) discussed in chapter III. Un- 
less explicitly stated otherwise, the results are for the 
case of a normally incident, spatially uniform, 1 Pa peak 
acoustic load and a . 01 N peak point vibrational load 
located at the coordinates of ot^=. 06033 m and <* 2 = 

.13547 m, and are for the case of the small damping model 
represented by equation (3.11). In all cases the forcing 
functions were applied uniformly over the frequency range at 
2 Hz intervals. 

The experimental results presented in this chapter were 
obtained using the apparatus and computer programs discussed 
in chapter IV. The results that demonstrate the effects 
produced by changes in the relative phase of the inputs were 
obtained by reversing the polarity of the shaker system. 

This was accomplished simply by switching the wire leads 
connecting the shaker to the input signal. This reversal of 
the polarity of the shaker is tantamount to changing the 
relative phase between the acoustic and vibrational inputs 
by 180 degrees from the existing phase difference at all 
frequencies . 

Several conditions were consistently maintained for all 
of the experimental results shown. The same input signal 
was applied to both the array of speakers and the shaker so 
that the two inputs were fully coherent. The input signal 
was filtered so that the input was approximately uniform 
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over the 0-1000 Hz frequency range. The overall gain fac- 
tors of the acoustic and vibrational inputs were adjusted so 
that the airborne sound power component would be dominant in 
some frequency ranges while the structureborne sound power 
component was dominant in other ranges. Also, in order to 
match the resolution of the analytical results, the FFT 
analyzer was set with a constant bandwidth of 2 Hz. Fin- 
ally, unless explicitly stated otherwise, the experimental 
results presented here were measured data on an aluminum 
plate (plate no. 1) for the case of an approximately uniform 
normal acoustic load and a point vibrational load located at 
the coordinates of a^= . 06033 m and o<2=. 13547 m. 

The results that follow are divided into five sections 
which examine the effects of several parameters on (1) the 
airborne and structureborne noise radiative characteristics 
of plates and (2) the proposed diagnostic method. In the 
first section (section B) , the effects of the relative mag- 
nitudes and phase of the inputs is examined for the baseline 
case of an aluminum plate. In section C, the effects pro- 
duced by changing the shaker location (altering the path of 
the structureborne input) are examined. Section D examines 
the effects of adding damping treatment to the aluminum 
plate. Finally, sections E and F examine the effects 
produced by altering the plate material to various types of 
composite materials. 
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B. EFFECTS OF RELATIVE MAGNITUDE AND PHASE 


Analytical Results for a Dominant Airborne Source 

Figures (V-la) and (V-lb) show the radiation efficien- 
cies and sound power levels obtained from the analytical 
program SOUND for the aluminum plate (plate no. 1) produced 
by individual airborne and structureborne inputs. Figure 
(V-la) clearly demonstrates that the airborne noise is much 
more efficient than the structureborne noise as expected. 
Figure (V-lb) shows that, in this case, the airborne source 
is dominant. 

The radiation efficiencies and sound powers of plate 
no. 1 produced by combining the airborne and structureborne 
inputs at 0 degrees in phase, 180 degrees out of phase, and 
90 degrees are shown in figures (V-lc) and (v-ld) respec- 
tively. These plots show that the sum of the results of the 
individual inputs is roughly equivalent to results obtained 
by combining the inputs at 90 degrees over most of the freq- 
uency range. The plots show that large deviations occur, 
however, when the inputs are combined at 0 or 180 degrees. 
This result underscores the importance of the cross term 
<TT > contained in equation (3.22) by showing that its 
contribution to the overall sound power radiated can be as 
large as the dominant term. Furthermore, since the sign of 
<TT >, can be positive or negative, the overall noise 


76 




generated can vary over a tremendous range depending on the 
relative phase of the inputs. 

Figures (V-le), (V-lf), and (V-lg) shows the results of 

applying program PREDICT (the proposed noise path separation 
method) to the cases of combined inputs at 90, 0, and 180 
degrees respectively. PREDICT uses a and a shown 
in figure (V-la) along with equations (3.29) and (3.31) to 
estimate the airborne and structureborne components of the 
total sound power in figure (V-ld). In each of figures (V- 
le) through (V-lg) the actual airborne and structureborne 
components are plotted against the components estimated by 
the proposed diagnostic method. In all cases the proposed 
separation method correctly identifies the airborne input as 
the dominant noise source over most of the frequency range. 

Figure (V-le) (uncorrelated inputs) shows good agree- 
ment between the curves over nearly the entire 0-1000 Hz 
frequency range. Discrepancies are seen to occur in this 
figure only in those ranges where a - a is very 
small. Good agreement is expected in this case since the 
cross terms are smallest when the inputs are uncorrelated. 

Figure (V-lf) indicates that some anomalies occur in 

the estimated components when the inputs are in phase. The 

larger differences in the curves are caused by the strong 

influence of <TT >. . For example, in figure (V-lf) 

s xa. c 

<TT >, exceeds the combined noise in the 800-900 Hz 

3 t 
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frequency range. Since it is known that the estimates are 
power preserving, this indicates that the <TT s > t has 
a negative value in this range. From equation (3.32) it is 
known that <TT >, can have negative values only if 
< TT o „_ > . is negative or if a_„_ is greater than o_. 

SXu u SXu ci 

From the figure (V-ld), it can be inferred that <7T >. 

sxa. t 

is negative in this case and larger than <TT >. . 

S3 L. 

Figure (V-lg) shows a similar irregularity in the 
800-900 Hz frequency range except that in this case 
<TT >. is negative. Since figure (V-ld) indicates 

Cl L. 

that <7T >. is positive, it can be inferred from equa- 

S Xa u 

tion (3.30) that or _ is less than a in this case. This 

S Xa 

implies that over 100 percent of <TT > is attributed 

s xa u 

to <TT > . 
s t 

In spite of the anomalies seen in figures (V-lf) and 
(V-lg), the proposed prediction method correctly identifies 
the dominant input over most of the frequency range and pro- 
vides valuable information about the interaction between 
inputs in those ranges where the cross terms are signifi- 
cant. It is worth noting that the 800-900 Hz frequency 
range where the cross term has such a significant effect is 
a region where there is strong cross coupling between dif- 
ferent modes, a - a is small, and <TT >. is 
as at 

nearly equal to <TT s >^. (See figures (V-la) through 
(V-ld). ) 
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Analytical Results for a Dominant Structureborne Source 

Figures (V-2a) and (V-2b) show the computed radiation 
efficiencies and sound power levels, respectively, produced 
by plate no. 1 due to a .2 Pa peak acoustic load and an in- 
dependent .05 N peak vibrational load. Figure (V-2b) shows 
that, in this case, the structureborne source is dominant. 
Comparing figure (V-la) and figure (V-2a) it is clearly seen 
that the radiation efficiencies are independent of the mag- 
nitude of the inputs. 

The radiation efficiencies and sound powers of plate 
no. 1 produced by combining these two inputs at 0 degrees in 
phase, 180 degrees out of phase, and 90 degrees are shown in 
figures (V-2c) and (V-2d) respectively. These figures again 
show that large deviations occur when the inputs are combin- 
ed at 0 or 180 degrees. Particularly good agreement between 
the sum of the individual inputs and the case of uncorrelat- 
ed, combined inputs is seen to occur in these figures. 

Figures (V-2e), (V-2f), and (V-2g) show the results of 

applying program PREDICT to the cases of combined inputs at 
90, 0, and 180 degrees respectively. In all three cases the 
proposed separation method correctly identifies the struc- 
tureborne input as the dominant noise source over the entire 
frequency range. Figure (V-2e) (uncorrelated inputs) shows 
good agreement between the component curves and the diag- 
nostic curves except in those ranges where a - a 

cl S 
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is small. Figure (V-2f) (inputs in phase) and figure (V-2g) 

(inputs out of phase) show much larger differences in the 

curves due to the extreme influence of <‘TT' >.. Note 

sxa t 

that <7T >. has a negative value over almost the 
at 

entire frequency range when the inputs are combined out of 
phase (see figure (V-2g)). The much more radical influence 
of <TT >. in figures (V-2f) and (V-2g) (compared to 
figures (V-lf) and (V-lg)) can be explained by the analysis 
in section H of chapter III where it was shown that the ef- 
fects of <TT >. are more likely to manifest themselves 

S XcL L. 

in <T T >. than in <TT >. . In spite of these 
at s t 

large differences, the prediction method correctly identi- 
fies the dominant input over the entire frequency range and 
once again provides valuable information about the inter- 
action between inputs in those ranges where the cross terms 
dominate. 

Experimental Verification of the Results 

Figures (V-3a) and (V-3b) show the measured radiation 
efficiencies and sound power levels, respectively, produced 
by plate no. 1 due to the acoustic and vibrational inputs 
acting independently (alone). Figure (V-3a) shows that the 
airborne sound power generation is more efficient than the 
structureborne sound power generation as predicted by the 
analytical study. Figure (V-3b) shows that, in this case, 
the structureborne component of the sound power is dominant 


80 



in several small discrete frequency regions while the 
airborne component is dominant over most of the frequency 
range. Note that the radiation efficiency and sound power 
curves appear to be somewhat more jagged than the smooth 
curves obtained in the analytical study. This change in the 
quality of the smoothness of the curves is quite under- 
standable in light of the many approximations employed by 
the analytical model and particularly since the inputs to 
the plate are only approximately uniform over the 0-1000 Hz 
frequency range. 

The radiation efficiencies and sound powers of plate 
no. 1 produced by combining the airborne and structureborne 
components are shown in figures (V-3c) and (V-3d) respec- 
tively. In each of these two figures, the three curves show 
the results of (1) summing the individual airborne and 
structureborne components, (2) combining the airborne and 
structureborne inputs with positive polarity on the shaker, 
and (3) combining the inputs with negative polarity on the 
shaker. These figures show that the sum of the results of 
the individual inputs is roughly equivalent to the results 
obtained by combining the inputs over much of the frequency 
range. Large differences in the curves are seen to occur, 
however, predominantly in frequency regions where the radia- 
tion efficiency curve reaches a local peak. These differ- 
ences are most certainly a reflection of the large infuence 
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of the cross term <TT > in those localized frequency 

sxd L. 

ranges . 

Figures (V-3e) and (V-3f) show the results of applying 
program PREDICT (the proposed noise path separation method) 
to the two cases of combined inputs shown in figure (V-3d). 
PREDICT uses a and a shown in figure (V-3a) along 
with equations (3.29) and (3.31) to predict the airborne and 
structureborne components of the total sound power in figure 
(V-3d). Figures (V-3e) and (V-3f) show that the proposed 
separation method correctly identifies the dominant noise 
source over most of the frequency range. It is evident, 
however, that because of the close proximity of a and 
a , the separation method was unable to make a reliable 
prediction in a number of ranges. Close comparison of the 
curves in figures (V-3e) and (V-3f) suggests that the cross 
term <TT > manifests itself primarily in the airborne 
estimate <TT a >£ as predicted by the analytical study. 

As expected, figure (V-3e) indicates that <TT >, 

3 L. 

overestimates the airborne component when <7T >, is 

large and positive while figure (V-3f) indicates that 
<TT >. underestimates the airborne component when 
<TT >. is large and negative. Note that the struc- 

S X3 l. 

tureborne estimate, <TT" >^, is quite accurate in both 
figures (V-3e) and (V-3f). 
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C;_ EFFECTS OF ALTERING THE INPUT PATHS ( SHAKER LOCATION) 
Analytical Results 

Figures (V-4a) through (V-4d) show the results that 
were obtained from program SOUND for the case of a normally 
incident uniform acoustic load and a point vibrational load 
located at the new coordinates of a^=. 12065 m and 
« 2 = • 20320 m. These coordinates place the vibrational 
point load at the center of the plate, thus driving exactly 
the same modes as the acoustic input (i.e. the odd modes of 
the plate). The small damping model, represented by equa- 
tion (3.11), was used again for this simulation. 

Figures (V-4a) and (V-4b) show the radiation efficien- 
cies and sound power levels, respectively, produced by plate 
no. 1 due to a 1 Pa peak acoustic load and an independent 
.01 N peak vibrational load. Figure (V-4a) shows once again 
that the airborne noise is more efficient than the struc- 
tureborne noise over most of the frequency range. The 
regions where a - a is small or where the two 
curves intersect are seen to cover a larger portion of the 
curves than for previous analytical cases, however. Figure 
(V-4b) shows that the airborne source is dominant or at 
least of equal influence over the entire frequency range for 
these inputs. 

The radiation efficiencies and sound powers produced by 
combining these two inputs are shown in figures (V-4c) and 
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(V-4d). The figures show that the overall noise radiation 
due to combined inputs is extremely sensitive to changes in 
the relative phases between the inputs. This result indi- 
cates that the cross term <7T >. has an even larger 

influence on the overall noise radiation if and when the 
airborne and structureborne inputs drive the same modes of 
the structure. 

Figures (V-4e) through (V-4g) show the results of ap- 
plying program PREDICT to the three cases of combined inputs 
at 90 degrees, 0 degrees, and 180 degrees respectively. In 
all three cases the proposed separation method correctly 
identifies the dominant noise source except in those ranges 
where a - a is small or zero. Note that figures 
(V-4f) and (V-4g) show very large differences between the 
estimated components and actual airborne and structureborne 
components. The differences in the curves are not so large, 
however, that the information contained in the estimates is 
deemed unusable. 

Experimental Verification of the Results 

Figures (V-5a) through (V-5d) show the measured data 
that were obtained from plate no. 1 for the case of an 
approximately uniform acoustic load and a point vibrational 
load located at the new coordinates of a^=. 12065 m and 
oi2 = . 20320 m. These coordinates place the vibrational 
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point load at the center of the plate, thus driving exactly 
the same modes as the acoustic input. 

Figures (V-5a) and (V-5b) show the radiation efficien- 
cies and sound power levels, respectively, produced by plate 
no. 1 due to independent acoustic and vibrational loads. 
Figure (V-5a) shows once again that the airborne noise is 
more efficient than the structureborne noise over most of 
the frequency range. Figure (V-5b) shows that the struc- 
tureborne component of the sound power is dominant in sever- 
al small discrete frequency regions while the airborne com- 
ponent is dominant over most of the frequency range. 

Smoother structureborne curves are expected in this case, 
since the shaker drives only the odd modes of the panel. 

Note that the structureborne radiation efficiency and sound 
power curves are indeed smoother than the curves obtained 
when the panel was driven near the corner. 

The measured radiation efficiencies and sound powers of 
plate no. 1 produced by combining these two inputs are shown 
in figures (V-5c) and (V-5d) respectively. These figures 
again show that the sum of the results of the individual in- 
puts is roughly equivalent to the results obtained by com- 
bining the inputs. Significant differences in the curves 
are seen to occur, however, in many frequency regions. 

These differences can again be explained by the influence of 
the cross term <TT >,. 
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Figures (V-5e) and (V-5f) show the results of applying 
program PREDICT to the combined inputs cases of figure (V- 
5d). These figures show that the proposed separation method 
correctly identifies the dominant noise source over nearly 
the entire frequency range. The separation method was again 
unable to make a reliable prediction in some ranges due to 
the small value of a - a . The comparisons are 
seen to be quite good with significant discrepancies occur- 
ing only in regions where a - a is small. 

5L S 

D. EFFECTS OF ADDED DAMPING 
Analytical Results 

Figures (V-6a) through (V-6d) show the results that 
were obtained from program SOUND for baseline conditions 
(shaker located at a^= .06033 m and a 2 =. 13547 m) on 
the aluminum plate with the moderate damping model repre- 
sented by equation (3.12). 

Figures (V-6a) and (V-6b) show the radiation efficien- 
cies and sound power levels, respectively, produced by plate 
no. 1 due to a 1 Pa peak acoustic load and an independent 
.01 N peak vibrational load. Figure (V-6a) shows a marked 
increase in the airborne radiation efficiency and consider- 
able smoothing of the structureborne radiation efficiency 
curve. The regions where a - a is small are seen 
to cover a much smaller band of frequencies than for 
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previous cases. Figure (V-6b) shows that the airborne 
source is highly dominant over the entire frequency range. 

The radiation efficiencies and sound powers produced by 
combining these two inputs are shown in figures (V-6c) and 
(V-6d). The figures show that overall noise radiation due 
to combined inputs is still sensitive to changes in the 
relative phases between the inputs. 

Figures (V-6e) through (V-6g) show the results of 
applying program PREDICT to the cases of combined inputs. 

In all three cases the proposed separation method correctly 
identifies the dominant noise source except in one small 
range in figure (V-6g). Figure (V-6e) (uncorrelated inputs) 
shows extremely good agreement between the curves over 
nearly the entire frequency range. Figure (V-6f) (inputs in 
phase) and figure (V-6g) (inputs out of phase) show moderate 
differences in the curves with the largest deviations occur- 
ing when the inputs are in phase. 

Experimental Verification of the Results 

Figures (V-7a) through (V-7d) show the measured results 
that were obtained from plate no. 1 with damping tape added 
to the plate. The self-adhesive damping tape added to the 
panel consisted of single layer of polystyrene type foam 
material with an outer layer of aluminum foil. The tape 
added approximately 1.44 kg/m to the surface density of the 
panel. This type of damping tape is commercially available 
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from several manufacturers and is routinely used to dampen 
vibration of the sidewalls of general aviation aircraft. 

Figures (V-7a) and (V-7b) show the radiation efficien- 
cies and sound power levels, respectively, produced by the 
damped plate due to independent acoustic and vibrational 
loads. Figure (V-7a) shows a marked increase in the air- 
borne radiation efficiency and a drastic smoothing of the 
structureborne radiation efficiency curve. The regions 

where a - a is small are seen to be confined to 
a s 

frequencies below 200 Hz with crossovers confined to the 
region below 100 Hz. Figure (V-7b) shows that the airborne 
source is dominant over nearly the entire frequency range. 

The radiation efficiencies and sound powers produced by 

combining these two inputs are shown in figures (V-7c) and 

(V-7d). Figure (V-7c) shows that the radiation efficiency 

of the combined noise can vary over an extremely large range 

of values depending on the phase relationship between the 

inputs. Likewise, figure (V-7d) shows that the combined 

sound power is largely influenced by the cross term 

<TT >. . 
sxa t 

Figures (V-7e) and (V-7f) show the results of applying 
program PREDICT to the two cases of combined inputs. In 
figure (V-7e) the proposed diagnostic method correctly 
identifies the airborne component as the dominant noise 
source. In figure (V-7f), however, the airborne prediction 
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is seen to be overwhelmed by the large negative cross term 

<7T >. . The structureborne estimate <7 T >. is 

sxa t st 

seen to exhibit close agreement with the independent struc- 
tureborne sound power < TI’ s > t in both cases, however. 

E. EFFECTS OF COMPOSITE MATERIAL PLATE CONSTRUCTION 
Analytical Results 

Figures (V-8a) through (V-8d) show the results that 
were obtained from program SOUND for baseline conditions 
(shaker located at a^= . 06033 m and 012 = . 13547 m with 
the small damping model) for the case of a symmetric tape 
ply graphite epoxy plate with a 0/90 lay-up (plate no. 2 ). 
(See Tables IV-1 through IV-3 for additional information on 
the characteristics of this plate.) 

Figures (V-8a) and (V-8b) show the radiation efficien- 
cies and sound power levels, respectively, produced by plate 
no. 2 due to a 1 Pa peak acoustic load and an independent 
. 01 N peak vibrational load. Figure (V-8a) shows consider- 
able smoothing of the radiation efficiency curves when com- 
pared with the baseline aluminum case. The explanation for 
this phenomenon can be found in the eigenvalue analysis 
which indicates that plate no. 2 has fewer modes with 
natural frequencies in the 0-2000 Hz range. Thus, the 
curves are smoother for the composite plate simply because 
there are fewer modes to cause the peaks and valleys. The 
regions where a - a is small or zero are seen to 

3. S 
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cover a larger band of frequencies than for the aluminum 
plate, however. Figure (V-8b) shows that, in this case, the 
airborne source is highly dominant over the entire frequency 
range. 

The radiation efficiencies and sound powers produced by 
combining these two inputs are shown in figures (V-8c) and 
(V-8d). The figures show that overall noise radiation due 
to combined inputs is still sensitive to changes in the 
relative phases between the inputs. 

Figures (V-8e) through (V-8g) show the results of 
applying program PREDICT to the cases of combined inputs. 

In all three cases the proposed separation method correctly 
identifies the dominant noise source except in those ranges 
where a - a is small. Figure (V-8e) (uncorrelated 
inputs) shows good agreement between the curves over nearly 
the entire frequency range. Figure (V-8f) (inputs in phase) 
and figure (V-8g) (inputs out of phase) tend to indicate 
much larger differences in the curves with the largest de- 
viations occuring when the inputs are out of phase. The 
differences in these figures between the estimated compo- 
nents and the actual components are not so large as to 
render the diagnostic information unusable however. 

Experimental Verification of the Results 

Figures (V-9a) through (V-9d) show the measured results 
that were obtained from plate no. 2. 
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Figures (V-9a) and (V-9b) show the radiation efficien- 
cies and sound power levels, respectively, produced by plate 
no. 2 due to independent acoustic and vibrational loads. 
Figure (V-9a) shows radiation efficiencies that are quali- 
tatively similar to the efficiencies of figure (V-3a) (plate 
no. 1). Futher scrutiny of figure (V-9a) tends to confirm 
that the composite plate possesses fewer resonance frequen- 
cies in the 0-2000 Hz range than the aluminum plate. The 
regions of figure (V-9a) where a - a is small or 
zero are seen to occur quite frequently. Figure (V-9b) 
shows that, for this case, the structureborne source is 
dominant in many discrete frequency regions. 

The sound power produced by combining the inputs with 
positive polarity on the shaker is shown in figure (V-9c). 
This figure shows that the cross term <Tf' >. can have a 
large influence on the combined sound power. 

Figure (V-9d) shows the results of applying program 
PREDICT to the combined inputs case shown in figure (V-9c). 
Figure (V-9d) shows that the proposed diagnostic method 
correctly identifies the dominant noise source except in 
those frequency regions where a - a is small or 
zero . 
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E\ RADIATION EFFICIENCIES OF COMPOSITE TAPE PLY PANELS 
Analytical Results 

The results of an analytical study of the radiative 
characteristics of several symmetric angle-ply composite 
tape panels (plates 3-6) are given in figure (V-10). (See 
Tables IV-1 through IV-3 for additional information on the 
physical characteristics of these plates.) The data shown 
were predicted by program SOUND for the case of the baseline 
conditions outlined in section A. 

Figure (V-lOa) shows the calculated radiation effi- 
ciencies of a symmetric 8 ply graphite tape panels with a 
+45/-45 lay-up. The radiation efficiency curves are seen to 
be similar in character to the curves in figure (V-8a) for 
the 0/90 lay-up. The curves in both figures are somewhat 
smoother than the radiation efficiency curves for the alumi- 
num panel. This is a strong indication that the graphite 
panels are stiffer than the baseline aluminum panel. The 
eigenvalue analysis performed indicates that both of these 
panels have 48 modes with natural frequencies in the 0-2000 
Hz frequency range whereas the aluminum panel has 68 modes 
with natural frequencies in this range. 

Figure (V-lOb) shows the calculated radiation efficien- 
cies of a 16 ply graphite tape panel with a +45/-45 lay-up. 
The radiation efficiencies curves in this case are smoother 
and show a smaller difference in the airborne and 
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struc tureborne curves. This is a direct consequence of the 
increased stiffness of the panel caused by the larger number 
of plies. The eigenvalue analysis performed indicated that 
this panel has only 26 modes with natural frequencies in the 
0-2000 Hz frequency range. 

Figure (V-lOc) shows the calculated radiation efficien- 
cies of an 8 ply Kevlar tape panel with a +45/-45 lay-up. 

The radiation efficiencies curves are seen to be remarkably 
similar to the curves obtained in figure (V-lOa) for the 8 
ply graphite tape panel with the +45/-45 lay-up. Close 
comparison of these two figures indicates that the Kevlar 
panel has lower resonance frequencies than the graphite 
panel. The greater compliance of the Kevlar panel was 
confirmed by the eigenvalue analysis which predicted 63 
modes with natural frequencies in the 0-2000 Hz frequency 
range compared to 48 for the graphite panel. 

Figure (V-lOd) shows the calculated radiation efficien- 
cies of an 8 ply fiberglass tape panel with a +45/-45 lay- 
up. The 2.21 kg/m surface density of the fiberglass panel 
is very close to that of the baseline aluminum panel and 
therefore invites comparison. The extreme dissimilarity of 
figures (V-la) and (V-lOd) indicates that the fiberglass 
panel has many more resonance frequencies in the 0-2000 Hz 
range than the aluminum panel. The greater stiffness of the 
aluminum panel was confirmed by the eigenvalue analysis 
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which predicted 103 modes with natural frequencies in the 0- 
2000 Hz frequency range for the fiberglass panel compared to 
68 for the aluminum panel. 

Experimental Verification of the Results 

The measured radiation efficiencies of several sym- 
metric angle-ply composite tape panels (plates 3-6) are 
given in figure (V-ll). First the results for the 8 ply 
graphite panel with a +45/-45 lay-up (plate no. 3) are shown 
in figure (V-lla). Secondly, the results for the 16 ply 
graphite panel with a +45/-45 lay-up (plate no. 4) are shown 
in figure (V-llb). Next, the results for the 8 ply Kevlar 
panel with a +45/-45 lay-up (plate no. 5) are shown in fig- 
ure (V-llc). Finally, the results for the 8 ply fiberglass 
panel with a +45/-45 lay-up (plate no. 6) are shown in fig- 
ure (V-lld). 

These measurements produced results that were qualita- 
tively similar to the results obtained in the analytical 
study. For example, figure (V-llb) indicates that the 16 
ply graphite panel was again the stiffest panel showing the 
fewest natural frequencies in the 0-1000 Hz range. Figure 
(V-lld) indicates that the 8 ply fiberglass panel was again 
the most compliant panel showing the greatest number of 
natural frequencies in the 0-1000 Hz range. All five com- 
posite panels exhibited marked differences in their respect- 
ive airborne and structureborne radiation efficiencies as 


94 



predicted by the analytical model with plate no. 4 showing 

the smallest difference (as predicted). Thus, both the 

analytical and experimental studies of the composite panels 

suggest that the coincidence phenomenon, which would result 

in little or no difference in a and ct , occurs at 

a s 

some frequency outside the 0-1000 Hz range. Therefore, the 
application of the proposed noise path separation method 
would not differ significantly if used on an aircraft fuse- 
lage constructed of composite materials. 

Gj_ DISCUSSION OF THE RESULTS 

To the author's knowledge, this is the first compre- 
hensive study of the problem of the sound radiation of air- 
craft type materials due to fully coherent combined airborne 
and structureborne inputs in the low frequency regime. The 
single most significant finding in the results of this study 
was the discovery of the relative importance of the cross 

term sound power component <TT >. . The importance of 

sxa. c 

this finding can not be over emphasized in light of the fact 
that all of the previous analytical and experimental studies 
performed in the last 7 or 8 years have neglected the inter- 
action between the airborne and structureborne inputs in the 
low frequency regime. By NASA's count (see reference 3), 
some 93 publications appeared in the open literature between 
1978 and 1984 which devoted their attention strictly to the 
airborne noise transmission through the sidewalls of 
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aircraft. A smaller number of studies (see references 18- 
22), addressing the problem of structureborne noise paths in 
aircraft, similarly restricted the scope of their analysis. 
The conclusions of some of these studies now seem dubious at 
best . 

Evidence to support this claim can be found in the 
results such as those found in figures (V-6) and (V-7). 
Computations on the analytical results in these figures (for 
the damped aluminum panel) show that while <TT >. = 82 
dB overall and <TT s > t = 68.6 dB overall, the levels for 
the combined sound power ranged anywhere from 80.7 dB to 
83.3 dB overall. Thus, the overall level of the combined 
sound power varied over a 2.6 dB range depending on the 
phase relationship between the airborne and structureborne 
inputs. Simple calculations show that the sum of individual 
components is 82.2 dB or, in other words, one would expect a 
0.2 dB increase in the overall sound power level due to the 
addition of the structureborne noise to the dominant air- 
borne noise rf the sound powers were additive . The implica- 
tion here is that the structureborne component does not have 
to radiate a significant amount of noise on its own in order 
to significantly change the level of the combined sound. In 
fact, the structureborne component needs only to change the 
dynamics of the problem in order to also significantly in- 
fluence the overall noise radiation. 
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The most significant finding with respect to the pro- 
posed noise path separation method was the extremely large 
influence of the cross term <7T >. on the airborne 

S Xa 

estimate <TT >. . This effect was observed in earlier 
a t 

experimental studies performed by the author on panels and 
on an aircraft fuselage (see references 39 and 40) but, at 
that time, was not fully understood. The analysis performed 
in section H of chapter III has shown that the effects of 
the cross terms are most likely to manifest themselves in 
the airborne estimate. The results of both the analytical 
study and experimental study clearly exhibit this behavioral 
aspect of equations (3.29) and (3.31). The results obtained 
from the measurements on the damped aluminum panel (see fig- 
ures V-6 and V-7) illustrate the effect in a particularly 
dramatic fashion by showing that the cross term <TT >. 
can be larger than the airborne component <Tr a >^ even 
when it is the dominant component of the noise. This result 
should serve as a warning to a potential user of this separ- 
ation method that one must be especially cautious in inter- 
preting the airborne estimate <TT >, . The fact that 

the airborne estimate <TT >, is often negative when 

at 

the inputs are out of phase should not be seen as a weakness 
in the diagnostic approach, however. By indicating a 
negative value, the estimates tell the measurer that a de- 
crease in one component of the sound power might actually 
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increase the total or combined sound power radiated. This 
is obviously a very useful piece of information to have when 
modifying the structure for noise control purposes. 

The results also indicate that the proposed separation 
method produces meaningful results irregardless of which 
source (airborne or structureborne) is dominant. Admit- 
tedly, however, the results may be more difficult to inter- 
pret when the structureborne noise is dominant since, when 
the inputs are correlated, the cross term tends to overwhelm 

the airborne estimate <TT >. . 

a t 

The results obtained when the structureborne and air- 
borne inputs drove the same modes (shaker at the center of 
the plate) indicate that there will be less difference in 
the airborne and structureborne noise radiative character- 
istics for this case and more frequency regions where a 

3 . 

- is small. This causes added difficultly in the 
implementation of the proposed separation method. It also 
appears that the effects of the cross terms are more severe 
in this case. These difficulties are primarily of academic 
interest, however, since the actual physical occurence of 
this case is unlikely. The most important aspect of these 
results is that they demonstrate the sensitivity of the 
structureborne radiation efficiency to changes in the point 
of application of the load. Comparison of figures (V-l) and 
(V-4) and comparison of figures (V-3) and (V-5) suggest that 


98 



a s is very path dependent. Thus, if the proposed separ- 
ation method is to be successful, the measurer must be cer- 
tain that path of the structureborne inputs during calibra- 
tion and measurement are the same. 

The results obtained on the damped aluminum panel 
(figures (V-6) and (V-7)) show that the principle effects of 
the added damping are to smooth the cr s curve, reduce the 
<TT s > t component, and to dramatically increase the air- 
borne radiation efficiency a . The results also suggest 
that the cross term <TT >. has a large influence in the 
damped case as well as the undamped case. The increase in 
a due to damping is caused by the reduction of the res- 
onance controlled component of the vibration (which is very 
inefficient). The importance of this result lies in its 
potential to aid the measurer in cases where the difference 

a - a is too small to utilize the proposed separ- 
a. s 

ation method. Comparison of figures (V-la) and (V-6a) dem- 
onstrates the large change in radiation efficiencies that 
can be realized through small changes in the damping of the 
structure . 

Finally, the analytical and experimental results ob- 
tained on the composite panels indicate that the noise 
radiative properties of these materials do not differ 
radically from those of aluminum in the 0-1000 Hz range. 

The results do suggest, however, that significant tayloring 
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of the noise radiative properties in isolated frequency 
regions can be achieved by adjusting parameters such as 
panel lay-up and the number of plies. This fact might be 
used to the noise control engineer's advantage if the 
interior noise of the aircraft in question has a discrete 
frequency character (such as in propeller driven aircraft). 
History has shown, however, that structural modifications 
such as this are more attuned to aircraft performance and 
load bearing considerations than to secondary considerations 
such as noise control. 

H. OBSERVED DIFFERENCES BETWEEN THEORY AND EXPERIMENT 

Analytical and experimental results for corresponding 
cases show considerable differences in the smoothness, 
location of resonances, and other characteristics of the 
radiation efficiency curves, etc... In this section, 
several possible sources of error in the analytical 
assumptions and in experimental data taking are discussed. 

Possible sources of error in the analytical modeling 
include (1) the eigenvalue and eigenfunction analysis, 

(2) the assumed boundary conditions, (3) the damping model, 
(4) the forcing function, and (5) the radiation conditions. 

The simple eigenvalue analysis used in this study is 
considered exact for the cases involving the aluminum plate 
and the symmetric ply laminates with specially orthotropic 
layers. No measure of the accuracy of equation (3.6) for 
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the eigenvalue analysis of the angle-ply laminates is avail- 
able. The author concedes that eigenvalues obtained using 
this equation are not as accurate as those which could be 
obtained by a finite element model or even a Rayleigh-Ritz 
analysis. The necessity for using an approximative equation 
such as (3.6) arises, however, as the number of natural fre- 
quencies which contribute to the noise generation in the 0- 
1000 Hz range becomes large. And as seen from the results, 
some of the plates have as many as 100 natural frequencies 
in the 0-2000 Hz range. This large number of natural freq- 
uencies makes a numerical analysis such as the finite ele- 
ment method impracticable. The use of the sine-wave eigen- 
functions could also cause some errors in the angle-ply lam- 
inate cases. These errors should be small, however, as long 
as the number of plies is large. 

The analytical model also assumes simply supported 
boundary conditions. This assumption is likely to cause 
prediction errors in the first few natural frequencies of 
the plates. Because the actual boundary conditions are 
somewhere in between clamped and simply supported, the true 
natural frequencies of the plates are expected to be higher 
than those predicted. The severity of this error decreases 
with increasing mode numbers. 

Errors introduced by the damping model are expected to 
have their largest impact on the airborne radiation 
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efficiency curves since it was shown that a a is very 
sensitive to damping. The damping model may also introduce 
some spurious cross coupling between modes which does not 
exist in the actual case under study. 

Differences in the actual and analytically modeled 
forcing functions could be the source of discrepencies in 
the results. During the experimental phase of the study, 
the shaker- rod-panel system was attached throughout the mea- 
surement process. This means that the actual forcing func- 
tions experience an additional impedance due to the attached 
mass of the shaker at the rod attachment point. The anal- 
ytical model tends to overestimate natural frequencies and 
overestimate the forced response because of this added 
driving impedance. A second source of error is the crude 
approximation of the spatial distribution of the forcing 
function. Small differences in the placement of the shaker, 
and differences in the approximately uniform airborne input 
are apt to cause large differences in the modal response of 
the structure. Finally, the fluid loading effects at the 
first few natural frequencies of the plate have been neg- 
lected in the forcing function. An effort was made to cor- 
rect for this error through the damping model, however. 

A final source of error in the analytical model is the 
assumption of free-field acoustic radiation conditions. 
Efforts were made to approximate these conditions in the 
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experimental portion of the study by the addition of the 
acoustic wedges to the receiving room of the TL apparatus. 
Realistically, however, the free-field conditions were 
probably achieved only above 200 Hz. 

Sources of error in the data taking include (1) the 
acoustic intensity measurement method, and (2) random and 
bias errors in the space averaging and time averaging. 

The most worrisome source of error (in this study) due 
to the acoustic intensity measurement method was near field 
effects. This type of error is most severe when the mea- 
surements are performed in a highly reactive acoustic field. 
This is precisely the type of noise field that is produced 
by a structureborne noise source. The existence of large 
time-averaged pressure gradients in these fields can cause 
large errors in the intensity measurements. A brief discus- 
sion this and several other types of error associated with 
the acoustic intensity method is given in Appendix XI. 

Random errors in the the space and time averaging can 
be controlled by an appropriate choice of the number of 
accelerometer locations and the size of the area chosen for 
analysis. The measurer should also be aware of certain 
types of bias error that enter into the data taking due to 
effects such as leakage and aliasing which are inherent in 
the discrete Fourier tranform process. Appropriate measures 
can also be taken to control the influence of these errors. 
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Figure (V-l) -- Analytical results for plate no. 1 with a dominant airborne source, the 

shaker located near the corner, and the small damping model. 
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Figure (V-l) 


Analytical results for plate no. 1 with a dominant airborne source, the 
shaker located near the corner, and the small damping model (cont'd). 
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Figure (V-2) 


Analytical results 
the shaker located 


for plate no. 1 with a dominant structureborne source 
near the corner, and the small damping model. 
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Analytical results 
the shaker located 


for plate no. 1 with a dominant structureborne source 
near the corner, and the small damping model (cont'd) 
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Figure (V-3) -- Experimental results for plate 

corner (cont'd). 
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Figure (V-4) -- Analytical results for plate no. 1 with a dominant airborne source, the 

shaker located at the center, and the small damping model (cont'd). 


RADIATION EFFICIENCY, dB RADIATION EFFICIENCY. dB 


(a) EFFC ENCIES OF THE COMPONENTS 




I i i 

'i'i 

\l i 


0.4 0.6 0.8 

(Thou*and») 

FREQUENCY. Hz 

s.a. a. a. 


(c) EFFICIENCIES OF THE COMBINED NOISE 



Figure (V-5) -- Experimental results for plate 

center . 
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Experimental results Cor plate 
(cont' d) . 
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Figure (V-6) 


Analytical results for plate no. 1 with a dominant airborne source, the 
shaker located near the corner, and the moderate damping model (cont'd) 
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Figure (V-7) -- Experimental results for plate no. 1 with the shaker located near the 

corner and with damping tape added to the plate. 


(e) D,ACNOS~CS - fiPJ'S COMBINED (,-) 


(f) DIAGNOSTICS - .'MPUTS COMBINED (-) 




Figure (V-7) 


Experimental results for plate no. 1 with the shaker located near the 
corner and with damping tape added to the plate (cont'd). 
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Figure (V-8) -- Analytical results for plate no. 2 with a dominant airborne source, the 

shaker located near the corner, and the small damping model. 
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Figure (V-8) -- Analytical results for plate no. 2 with a dominant airborne source, the 

shaker located near the corner, and the small damping model (cont'd). 
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Figure (V-10) -- Analytical airborne and structureborne radiation efficiency curves for 

the symmetric angle ply composite panels (plates 3-6). 
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Figure (V-ll) -- Measured airborne and structureborne radiation efficiency curves for the 

symmetric angle ply composite panels (plates 3-6). 


Chapter VI 

SUMMARY AND CONCLUSIONS 

A_;_ REVIEW OF THE OBJECTIVE 

The purpose of this study was to develop a new diag- 
nostic measurement method, based on the acoustic intensity 
measurement method and on several equations developed by the 
author, for separating and predicting the airborne and 
structureborne components of the sound power radiated by 
aircraft- type panels due to some unknown combination of 
acoustic and vibrational inputs. In order to validate the 
proposed method, both analytical and experimental studies 
were undertaken. The purpose of the analytical study was to 
provide the theoretical basis for the method and simulate 
its behavior under various operating conditions. The pur- 
pose of the experimental study was then to verify the ex- 
pected behavior of the proposed diagnostic method. 

_B;_ MAJOR FINDINGS OF THE STUDY 

The most significant finding of this study was the 
extremely large influence that the relative phase between 
the inputs has on the combined noise radiation of the 
plates. It has been shown that phase dependent effects 
manifest themselves as cross terms in both the dynamic and 
acoustic portions of the analysis. Both the analytical and 
experimental studies show that these cross terms can 
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radically alter both the combined sound power and combined 
radiation efficiencies of plates constructed of aircraft- 
type materials. 

The diagnostic method proposed in this study has been 
used successfully to predict the relative contributions of 
the airborne and structureborne sound power components of 
fully coherent combinations of simultaneous acoustic and 
vibrational inputs to simple panels. The use of radiation 
efficiency measurements for separating the airborne and 
structureborne components of the total sound power radiated 
in aircraft appears to be a viable option to the other 
diagnostic tools currently available. The results of the 
study indicate that the method accurately predicts the domi- 
nant noise source, alerts the measurer to any strong inter- 
active components, and can be applied to plate or thin shell 
structures constructed of a variety of materials. 

Finally, both analytical and experimental studies on 
the composite panels show that, although the noise radiative 
characteristics of these materials are not radically dif- 
ferent from those of the aluminum panel, the character of 
the radiation efficiency curves can be surprisingly sensi- 
tive to panel lay-up. This is an important finding since it 
suggests that the designer might be able to selectively mod- 
ify the noise transmissive properties of the sidewalls of an 
aircraft with little or no weight penalty. 
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C. APPLICABILITY OF THE DIAGNOSTIC METHOD TO AIRCRAFT 

Although the proposed noise path separation method was 
not applied to an actual aircraft structure in this study, a 
preliminary study of the feasibility of measuring radiation 
efficiencies in aircraft was performed by the author in an 
earlier experimental study. (See reference 40.) This earl- 
ier study demonstrated that the method has considerable pro- 
mise for application in actual aircraft structures. The 
results of that earlier study coupled with the rigorous 
analysis performed in this study warrant the further invest- 
igation of this method as a noise source-path identification 
tool in aircraft. 

The major advantage of the proposed method over the 
other diagnostic approaches currently available is that this 
new method does not necessarily require any modifications to 
the aircraft and could be used during flight. Also, the 
method can be used for fully coherent acoustic and vibra- 
tional inputs, whereas the partial/multiple coherence func- 
tion methods can not, and the method works best in the low 
frequency regime where lead-wrapping fails because of the 
low transmission loss of lead. 

Limitations of the proposed method include its low 
frequency restriction and its inapplicability when a 
cr g is small. The method can only be used at frequencies 
below the critical frequency of a thin shell structure. At 
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frequencies above the critical frequency, there will be 

little or no difference in a and a . Below the 

a s 

critical frequency, the difficulties of a small value of 

a - a can be easily overcome, however, with the 
a s 

application of additional damping to the structure. 

D;_ AREAS FOR FUTURE RESEARCH 

Future research on the proposed diagnostic method and 
its applicability to aircraft should be directed at the 
study of the certainty of radiation efficiency measurements 
in aircraft. Studies should be performed to determine the 
sensitivity of radiation efficiency measurements to such 
factors as (1) changes in the acoustic source directivity, 

(2) changes in the path of the vibrational energy flow, and 

(3) changes in the number of acoustic and vibrational 
inputs . 

The applicability of the proposed separation method to 
aircraft will remain in question until these problems are 
addressed. For example, a complex structure such as an 
aircraft may have several important structural transmission 
paths. By measuring a s for each structural transmission 
path, equations (3.29) and (3.31) become a system of two 
equations with many more than two unknowns. Thus, measuring 
structureborne radiation efficiencies for many structural 
transmission paths would not necessarily solve the problem 
since the relative importance of the different structural 
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paths may be unknown. As another example, suppose a loud- 
speaker is used as an acoustic source to determine the air- 
borne radiation efficiency of some small area on a propeller 
driven aircraft fuselage. The acoustic directivity of the 
loudspeaker may be dissimilar to the directivity of the 
actual acoustic sources on the aircraft (that is the pro- 
pellers). This approximation of a could introduce 
considerable error into the structureborne sound power 
prediction <TT >. (as indicated by the error analysis 

b U 

in Appendix VIII). Hypothetical situations such as these 
encourage the investigator to perform some preliminary 
measurement and analysis on the particular problem of 
interest before attempting to use the proposed separation 
method. 
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Appendix I 

BACKGROUND INFORMATION AND DEFINITIONS 

A^_ FUTURE TRENDS IN AIRCRAFT DESIGN 

Recent research on fuel efficient aircraft designs has 
focused on the new advanced turboprop (ATP) propulsion sys- 
tems (sometimes called propfans) which have been proposed 
for future use in both the commercial and military fleets in 
the 1990 's and beyond. A model of one of these advanced 
turboprop propulsion systems is shown in figure (AI-1). 

(See reference 3.) These propulsion systems were first pro- 
posed for use in high speed aircraft in 1974, and are char- 
acterized by their unusual number of blades (8 or more), 
their use of advanced airfoil design, and their use of high- 
ly swept blades (which reduce the blade tip Mach No. of the 
blades thereby increasing the propulsive efficiency of the 
design). They have since been the subject of numerous air- 
craft studies as a part of the NASA Aircraft Energy Effi- 
ciency Program. The results of these studies suggest that, 
through the use of these advanced turboprops, it may be pos- 
sible to reduce the average fuel consumption of an 0.8 Mach 
No. advanced aircraft by 20 to 40 percent, when compared to 
aircraft equipped with high bypass turbofan engines of equi- 
valent technology (see references 4 and 5). At the same 
time, one of the potentially severe problem areas that was 
identified in these studies was the probable high interior 
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noise levels that could be expected in an aircraft with this 
type of turboprop propulsion system. 

Another technological development which is expected to 
have a dramatic impact on the future design of aircraft is 
the emergence of composites materials. Current composite 
applications to large transport aircraft include secondary 
structure such as the wings and empennage as well as fair- 
ings and nacelle components (see reference 6). Future air- 
craft can be expected to utilize significant quantities of 
composites in the primary structures as well. A recent 
study on an all composite fuselage which meets all the 
design requirements for a 1990 's large transport aircraft 
estimated a 32 percent weight savings for the fuselage when 
compared to a baseline aluminum shell. (See reference 7.) 
Roughly 1/2 of this weight savings was in the design of skin 
panels alone. 

The trends toward all composite fuselage structures is 
considerable cause for concern for the noise control engi- 
neer since relatively few studies have been performed to 
quantify the noise radiative properties of these materials 
(see references 8-14), and because the results of prelimi- 
nary studies suggest that composites (because they are 
lighter) transmit airborne noise more readily than their 
metal counterparts. (See references 13-14). Preliminary 
studies also indicate that the noise transmissive properties 



of composite panels are surprisingly sensitive to parameters 
such as panel lay-up or fiber orientation. (See reference 
14.) 

CHARACTERISTICS OF AIRBORNE PROPELLER NOISE 

The sound pressure field produced by a propeller is 
characterized by its complex periodic nature. The pressure 
pulses produced by the propeller repeat themselves (with 
some unsteadiness) with each rotation of the propeller 
shaft. This periodic nature of the signal produces a sound 
pressure spectrum which consists of discrete tones, begin- 
ning with the fundamental blade passage frequency (propeller 
shaft rpm times the number of blades) which is lowest in 
frequency, followed by tones which have frequencies that are 
integer multiples of the fundmental (referred to as the pro- 
peller harmonics). An example of both the time history and 
the spectrum of the pressure signal which is produced by a 
ATP propulsion system is shown in figure (AI-2). (See ref- 
erence 15.) The predicted and measured data in figure (AI- 
2) are for a two foot diameter model of an ATP tested at 0.8 
Mach No. cruise speed. The propeller was mounted on top of 
an aircraft fuselage and microphones were flush mounted in 
the fuselage. As suggested by the data, the fundamental 
blade passage frequency (BPF) and the first few harmonics 
are the most significant of the tones. Aircraft interior 
noise studies have indicated that the higher harmonics 
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(tones with frequencies above the fifth harmonic of the BPF) 
are overwhelmed by other sources of noise in the cabin such 
as the turbulent boundary layer noise (see references 1-2). 
Thus, the crux of the propeller noise problems lies in the 
low frequency regime. And for nearly all propeller aircraft 
configurations presently under consideration, the fundament- 
al and first few harmonics of the propeller BPF fall in the 
0 - 1000 Hz frequency range. 

Another important aspect of the sound pressure field 
produced by a propeller is its highly directional character. 
The spatial characteristics of the sound field show a strong 
preference for sound radiation in the plane of the propeller 
blades. Figure (AI-3) shows the spatial sound field distri- 
bution for a high speed propeller (ATP) in flight. (See 
reference 16.) Theoretical predictions and measurements 
such as these indicate that the vast majority of the sound 
which impinges on the fuselage of an aircraft is essentially 
confined to an area within one propeller diameter on either 
side of the plane of the propeller. Figure (AI-4) (from 
reference 17) shows that the preferred direction of radi- 
ation of a typical ATP is at an angle of roughly 45 degrees 
in the static configuration. The figure shows that the 
effect of forward flight is to move the region of peak rad- 
iation to a position in space that is close to the plane of 
the propeller. The effect of the boundary layer on the 
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fuselage wall is to further bend the apparent angle of inci- 
dence of the incoming sound waves. Thus, the apparent angle 
of incidence of the noise impinging on the fuselage is gen- 
erally oblique. The wavelengths of the incoming sound waves 
are typically very long (1 to 3 meters) due to the low 
frequency character of the propeller noise. Thus, despite 
the fact that the angle of incidence of the waves is 
oblique, the distribution of pressure over the surface of 
the fuselage (with respect to phase) is slowly varying. 

C^_ CHARACTERISTICS OF STRUCTUREBORNE PROPELLER NOISE 

Since very few studies have been funded in the area of 
structureborne noise in propeller driven aircraft, the 
characteristics of structureborne noise are not as well 
understood and therefore are more difficult to generalize. 
Two separate studies recently completed by NASA and its 
contractors (references 18 and 19-21) and a third on-going 
study (reference 22), however, have identified two major 
contributors to the structureborne noise in propeller driven 
aircraft. Studies completed by Southwest Research Corpor- 
ation (references 19-21) indicated that, at least on some 
single engine propeller aircraft, roughly 1/2 of the noise 
enters the aircraft cabin via vibrational energy that begins 
as engine vibration and propagates along the engine mounting 
system and into the aircraft cabin where it is ultimately 
radiated as noise. A second study performed in house at 
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NASA Langley Research Center (reference 18) indicated that 
on some wing mounted propeller airplanes, the propeller 
downwash effects on the wing can be a significant source of 
interior noise. The study suggests that the trailing vor- 
tices shed from the propeller tips collide with the wing on 
each rotation of the propeller shaft. This sets up vibra- 
tion in the wings which propagates along a structural path 
into the aircraft cabin. 

The frequency character of the structureborne noise is 
generally the same as the airborne noise. Both sources of 
noise manifest discrete frequency tones at the BPF and its 
harmonics. Past research on aircraft with reciprocating 
engines (see references 1-2) has shown that quite often, 
even the piston firing harmonics of the engine coincide in 
frequency with the propeller harmonics. 

D^_ SEPARATION OF AIRBORNE AND STRUCTUREBORNE NOISE 

Three measurement methods which have most recently been 
proposed for separation of airborne and structureborne noise 
in aircraft have included: 

1) Measuring noise levels in the cabin with the engines run- 
ning normally and with the engines detached. (See refer- 
ences 19-21.) 

2) The conventional lead wrapping technique. (See refer- 
ences 18 and 23-24.) 
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3) The partial coherence function method. (See references 
25-27 . ) 

The obvious limitation of the first method is that the 
tests cannot be performed in flight. Detaching the engines 
may also prove to be a prohibitively expensive method for 
determining the relative contributions of the airborne and 
structureborne noises for a diverse fleet of aircraft. 

Attempts to identify noise paths in aircraft using the 
second method (the lead wrapping technique) have been large- 
ly unsuccessful because of the poor transmission loss of 
lead at low frequencies. As discussed earlier, the distin- 
guishing characteristic of the noise generated by a propel- 
ler is its low frequency tonal nature. 

The use of partial and multiple coherence function 
methods to determine the principle noise paths in aircraft 
proves inadequate if the various noise generating mechanisms 
of the aircraft are fully coherent. Unfortunately, for many 
of the aircraft in question, the airborne and structureborne 
components are fully coherent since they both ultimately 
have the same source, viz. the propellers. Thus, because of 
the fully coherent nature of the noise propagated along 
these two paths, the problem of separating and predicting 
the relative contributions of the airborne and structure- 
borne noises in propeller driven aircraft has remained 
largely unsolved. 
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E. ACOUSTIC QUANTITIES 


p(r,t) = 
v(r,t) = 
a(r,t) = 


instantaneous acoustic pressure as a 
function of position in space, 
instantaneous acoustic fluid particle 
velocity vector as a function of position, 
instantaneous acoustic fluid particle 
acceleration vector as a function of position. 


Fj_ FOURIER AND INVERSE FOURIER TRANSFORMS 

The Fourier and inverse Fourier transforms shall be 
defined as 


X(f) = F{x( t) } = 


- j2irft 

x(t) e dt , 


(Al.l) 


x(t) = F -1 {X ( f ) } = 


j2irf t 

X(f) e df 


(Al. 2) 


G. FINITE FOURIER TRANSFORMS 


In general 


X(f) = F{x( t) } = 


- j2irf t 

x(t) e dt . 


(A1.3) 


The finite Fourier transform of the acoustic pressure is 
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given by 


T 


P(r,f) = F (p(r, t) } = 


-j2irft 

p(r,t) e dt . 


The finite Fourier transform of the acoustic fluid 
velocity is given by 


V(r, f ) = F{v(r, t) } = 


+ - j2irf t 

v(r,t) e dt . 


0 


The finite Fourier transform of the acoustic fluid 
acceleration is given by 


A(r , t) = F{a(r, t) } = 


-j2irft 

a(r,t) e dt . 


H. RELATIONSHIP BETWEEN VELOCITY AND ACCELERATION 


a(t) = d/dt{v(t)} 


A ( f ) = - j2irfv( f ) = ju)V(f) . 


V(f) = - j/( 2irf ) £(f) = -j/u) X(f) 


I. AUTO CORRELATION FUNCTIONS 


In general 


R (t) 

XX 


lim 1/T 


x(t)x(t+T) dt 


(A1.4) 

particle 

( Al . 5 ) 

particle 

(A1.6) 

(A1.7) 
(Al.8) 
(Al. 9) 

(Al. 10) 
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An estimate of R (t) with finite data is 

xx 

T-x 

R xx C x ) = l/(T-x) x ( t ) x ( t+ x ) dt . (Al.ll) 

0 

Important properties : 

If the variable x(t) is real, then R (x) is a real, 

X X 

even function of x. 

CROSS CORRELATION FUNCTIONS 
In general 

T 

R (x) = lim x(t)y(t+x) dt . (Al.12) 

T-voo J 
0 

An estimate of R (x) with finite data is 

x y 

T-X 

R X y ( T ) = 1/ (T-x ) x(t)y(t+x) dt . (A1.13) 

0 

The cross correlation between acoustic pressure and fluid 
particle velocity is given by 

T 

R (x) = lim 1/T p(t)v(t+x) dt . C Al .14) 

P rp+CO J 

0 

Evaluating the preceding equation at x=0, one obtains 

T 

R (0) = lim 1/T p(t)v(t) dt . (Al.15) 

T->oo J 
0 
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Important properties: 

If variables x(t) and y(t) are real, R (x) is also 

xy 

real. R (x), in general, is neither an even nor an 
xy 

odd function of x. 


K^_ SPECTRAL DENSITY FUNCTIONS 

The auto spectral density function is defined by 


s * x ( £) = 


-j2irfx 

R xx (x) e dT • 


(Al.16) 


From the theory of Fourier transforms (see reference 52), 


R XX^ - 


j 2 XT f X 

S XX (T) e d£ . 


(A1.17) 


The cross spectral density function is defined by 


S (f ) = 
xy v 


- j 2 tt f x 

R xy (x) e dx . 


(A1.18) 


From the theory of Fourier transforms (see reference 52), 


V T) * 


j2nf T 

S xy (f) e df 


(A1.19) 


Important properties: 

Given that R (x) is a real function, S (f) is a 

XX X X 

real, even function of f. Given that R xy (x) is real, 
S X y(f) ^- n general, be complex with its real part 
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being a even function of f and its imaginary part being an 

odd function of f. Note that S (f) is defined for both 

xy 

positive and negative frequencies. For practical problems 
of interest, negative frequencies have no physical interpre- 
tation . 


ONE-SIDED AUTO SPECTRAL DENSITY FUNCTIONS 
Recall that by definition 

j2irf t 

s xx (f) e df . 

Given that R (t) and S (f) are both real, even 

aX XX 

functions 


R (t) = 
xx ^ ' 


R (t) 

XX v 1 


S (f) COS ( 2 TT f T ) df . 
XX 


Since the integrand in (Al.21) is even, 


(Al.20) 


(Al.21) 


R (x) 
xx v 


= 2 


S (f ) 
xx ^ 


cos ( 2 TT f T ) df 


2S (f ) 

XX 


j2irf t 


df 


(Al.22) 


Thus it is convenient to define the one-sided auto spectral 
density as 
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( 0<f <°°) . 


(Al. 23) 


G (f) = 2S (f) 

XX XX 


Thus , 


R xx< T ’ * 


j2uf x 

G xx (£) e d£ . 


(Al. 24 ) 


And by a similar analysis, 


G C f ) = 2 
xx J 


- j2irf x 

R xx (t) e dt . 


(A1.25) 


ONE-SIDED CROSS SPECTRAL DENSITY FUNCTIONS 
Recall that by definition 


R xy<^ * 


j 2irf x 

s X y( f ) e df . 


(A1.26) 


Given that R^Cx) is real, it follows that 


R (x) = 
xy ^ J 


Re{S }cos(2irfx) - Im{S }sin(2irfx) df 
xy xy 


Re{S ( f ) }cos ( 2irf x ) df 
xy 


Im{Sxy ( f ) }sin(2irfx ) df 


(Al.27) 
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R xy (T) 3 


2Re{S ( f ) }cos (2irft ) df 
xy 


2Im{S (f )}sin(2Ttfx) df . 
xy 


Therefore , 


R (t) = 
xy 


j 2 tt f T 

2S xy (f) e df . 


Thus, the one-sided cross spectral density is defined 


G (f) = 2S (f) 
xy xy v 


(0<f<«) 


Therefore, 


R (t) = 
xy 


j2irf t 

G xy (f) e df • 


And by a similar analysis 


G (f) = 2 
xy 


- j2irf x 

R xy ( t ) e dx . 


Furthermore, the co-spectral density C (f) and the 
quad-spectral density Q (f) are defined as 


(A1.28) 


(Al. 29) 
as 

(Al. 30) 


(Al. 31) 


(Al. 32) 



(Al. 33) 


G ( f ) = C ( f ) - jQ ( f ) , 
xy xy J xy ’ 

where C and Q are both real so that 
xy xy 


oo 

C ( f ) cos ( 2rf x ) + Q (f )sin(2irfx) 

xy xy 

o 


df . 


(Al. 34) 


N^_ EXPECTED VALUE OPERATOR 

The expected value of any real single valued continuous 
function g of the random variable x is defined as 


E{g(x) } 


00 

g(x)p(x) dx , 

— 00 


(Al. 35) 


where p(x) is the probability density function of the random 
variable x. (See reference 53, pp. 28-34 or reference 54, 
pp. 83-97 and pp. 165-173 for a discussion of the probabil- 
ity density function of a random variable.) 


0. ESTIMATION OF THE SPECTRAL DENSITY FUNCTIONS 

The one-sided spectral density functions can be defined 
in terms of the Fourier transform as (see reference 55) 

G (f) = lim 2/T E{X* ( f )X ( f ) } , (Al.36) 

T-*oo 

and 

G (f) = lim 2/T E {X* ( f ) Y ( f ) } . (Al.37) 

x y t -*- 00 
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Estimates of the spectral density functions are given by 
G xx ( f ) = 2/T |x(f)| 2 , (A1.38) 

and 


G xy (f) = 2/T 


X (f)Y(f) 


(Al.39) 


where X(f) and Y(f) are Finite Fourier transforms. 


P^_ ACOUSTIC INTENSITY 

The instantaneous acoustic intensity vector at a point 
in space is defined by 

I(r,t) = p(r,t)v(r,t) . (Al.40) 

The time averaged acoustic intensity vector is given by 


<1 (r ) > . = 


lim 1/T 
T "^ 00 


p(r,t)v(r,t) dt 


(A1.41) 


Hence 

<l(r)> = <p(r )v(r ) > t . (Al.42) 

Letting the spatial dependence be implicit and confining the 
analysis to the magnitudes of the vectors, 


<I> 


t 


lim 1/T 

‘P-VOO 


T 

p(t)v(t) dt . 

0 


(Al. 43 ) 
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Recalling the definition of the cross correlation function 


R ( t ) = lim 1/T 


p ( t ) v( t+x ) dt , 


and substituting, 


(Al.44) 


<I> t = R pv ( 0 ) . (A1.45) 

Recalling the definitions of the co-spectral and the quad- 
spectral density functions and their relation to the cross 
correlation function (equation (A1.34)), 



00 


(t) = 


C pv (f )cos(2irfx) + Q pv (f )sin(2TTfx) 


df , 


0 

and substituting x = 0 into the last equation, 


(A1.46) 


<I> t = R pv (0) = 


C (£) 

pv v 


df = Re{ 


G (f ) 
pv v 


df} 


(Al. 47) 


Thus, the magnitude of the time averaged acoustic intensity 
vector is the real part of the cross spectrum between 
pressure and particle velocity summed over all frequencies. 


Q_:_ SOUND POWER 

The instantaneous sound power, IT, radiated by a sound 
source is defined as 
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TT( t) 


f(r,t)*n dS , 


(Al.48) 


where I is the instantaneous acoustic intensity vector, 
n is a unit vector normal, and the integral is over a 
control surface enclosing the sound source. The time 
averaged sound power is given by 


<TT> = lim 1/T 


TT(t) dt 


= lim 1/T 


I-n dS 


dt . 


(Al.49) 


Interchanging the order of integration, 


<TT> . 


<f*n> t dS . 


(A1.50) 


ACOUSTIC IMPEDANCE 

The specific acoustic impedance, Z, is defined by 
Z(r,f) = P(r, f )/v(r, f ) . (Al.51) 

Important properties: 

In general, Z is complex and is a function of both frequency 
and position in space. At a great distance from the source 
of an acoustic disturbance, acoustic pressure and fluid par- 
ticle velocity are in phase and Z = P 0 c Q . 
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S. MEAN SQUARE SURFACE VELOCITY 


The space- time averaged mean square surface velocity of 
a vibrating surface is given by 


< v 


2 


> 

r 


t 


1/S 


lim 1/T 

T-*ao 


T 

v 2 (r , t) 

0 


dt 


dS , 


(A1.52) 


where S is the surface area of the structure. Recalling the 
definition of the correlation function and substituting, the 
result is 


<v 2 > rit = 1/S 


R vv (0) as = < 


G (f) df > . 

; r 


(A1.53) 


ACOUSTIC RADIATION EFFICIENCY 

The acoustic radiation efficiency of a vibrating 
surface, a, is defined as 


a = <TT> t /(p o c o <v > r#t S) 


(Al. 54) 


Recall from section M. that <TT> t is given by 
<TT> 1 _ = | | <I*n> dS = <I*n> . S , 

t IT j w 


(Al. 55) 


where the chosen control surface is the surface of the 
structure, and the intensity vector is measured at this 
surface. Substituting the previous result into the first 
equation of this section the result is 


= <1 > ./(p c <v 2 > . ) 

n r, t' ^ M o o r , t J 


(Al. 56) 
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where it should be understood that I represents the 
magnitude of the acoustic intensity vector measured normal 
to the vibrating surface of the structure. 


CRITICAL FREQUENCY 

The coincidence frequency is that frequency of vibra- 
tion at which flexural waves in the plate travel (in a given 
direction) at the same speed as acoustic (compressional) 
waves in air. The critical frequency is that special case 
of the coincidence frequencies for which the wavelength of 
the flexural waves in the plate are also equal to the wave- 
length of the acoustic waves (which is the case of grazing 
incidence of the acoustic wave). The critical frequency is 
important because it is lowest in frequency of all the coin- 
cidence frequencies. The critical frequency for an infinite 
flat isotropic plate is given by the equation 


f 

c 


= Cq/(2tt) 


ph 12 ( 1-p 2 ) /Eh 3 



(Al. 57) 


where 

c q = the speed of sound in the acoustic medium, 

p = the density of the plate material, 

h = the plate thickness, 

E = the modulus of elasticity (Young's modulus), 

y = Poisson's ratio for the plate material. 

For an orthotropic or an anisotropic plate, the critical 
frequency depends on the direction of travel of the wave. 
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ACOUSTIC SIGNATURE AND SPECTRUM OF 
AN ADVANCED PROPELLER 


□ EXPERIMENTAL * EXPERIMENTAL, OASPL= 145.0 dB 

o THEORETICAL I THEORETICAL, OASPL = 144.7 dB 



TIME/PERIOD HARMONIC NUMBER 

Figure (AI-2) — Noise signature of an advanced turboprop in a quiet wind tunnel (from 

reference 15). 
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JETSTAR FUSELAGE SURFACE NOISE 

8 - BIADED SR - 3 PROPFAN 



-2.0 -1.6 -1.2 -.8 -.4 0 .4 .8 1.2 1.6 2.0 

X/D 

Figure (AI-3) -- Spatial sound field distribution for an advanced turboprop (from 

reference 16). 
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Visual position Radiation position 



Microphone 


Figure (AI-4) -- Sound field directivity and angle of incidence for an advanced turboprop 

(from reference 17). 


Appendix II 

FREE VIBRATIONAL RESPONSE OF A SIMPLY SUPPORTED 
RECTANGULAR ORTHOTROPIC PLATE 

A. GOVERNING DIFFERENTIAL EQUATION 

The general free vibrational response of a thin shell 
structure is governed by the system of equations (see 
reference 56, eqn. 8.1.2) 

-L i (u 1 ,u 2 ,u 3 ) + X i u i + phu i =0 (i = 1,2,3) , (A2.1) 

where u^ = u^ ( , a 2 , t) . 

L^ = linear stiffness operator for the ith equation. 

= coordinates in 3-D space for the shell or plate. 
a^,a 2 ,a 3 form an orthogonal curvilinear 
coordinate system with and a 2 being the 
surface coordinates located on the neutral surface 
and everywhere tangent to the midsurface of the 
shell. is everywhere outwardly normal to the 

shell . 

u^ = displacement of the structure in the dir- 
ection. 

p = mass per unit area of the shell, 
h = thickness of the shell. 

= equivalent viscous damping factor in the ith 
direction . 
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It is well known in the theory of vibration that the natural 
frequencies associated with the transverse (bending) vibra- 
tional mode shapes occur at much lower frequencies than 
their corresponding in plane (compressional) vibrational 
counterparts. (See reference 56, pg.97). It is also well 
known in the theory of sound radiation from structures that 
the transverse (bending) vibrational modes of the structure 
couple much more strongly to the sound field than their in 
plane (compressional) vibrational counterparts. (See refer- 
ence 57, pg. 116). For these reasons we devote our atten- 
tion solely to the transverse or flexural (bending) vibra- 
tional response of the structure. The system of equations 
then simplifies to a single equation, viz., 

-L 3 (u 3 ) + Xu 3 + phu 3 = 0 . (A2.2) 

B^_ FREE RESPONSE FOR THE UNDAMPED CASE 

Consider a rectangular, orthotropic flat plate of 
dimensions d^ and d 2 in the cartesian coordinate system 
a^,a 3 ,a 3 . The stiffness operator L 3 (for 
flexure) in this case is given by (see reference 56, pp. 
341-354, or reference 58, pp. 140-146) 



D i;L 3 4 /9aJ + 2(D 12 + 2D 66 )3 4 /3aj32 + D 22 3 4 /3a 4 , (A2.3) 
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where the are the plate rigidity constants that 

relate the internal bending and twisting moments of the 
plate to the twists and curvatures they induce. For the 
flat orthotropic plate 


Dn = 

h 3 /12 

( 1 - y 2 ) y 

(A2. 4) 

D 22 = 

h 3 /12 

E 2 / d-yiy 2 ) * 

(A2.5) 

°12 = 

h 3 /12 

E l w 2 // ^ 1_y l lJ 2' 1 = h3 / 12 E 2 p 1^^ 1_,j 1 vi 2 ^ ’ 

(A2.6) 

°66 = 

h 3 /12 

G , 

(A2. 7) 

where 





= the effective modulus of elasticity (Young's 
modulus) in the ith direction, 

= the effective Poisson's ratio in the ith 
direction , 

G = the shear modulus of elasticity (modulus of 
rigidity) , 

h = is the thickness of the plate. 

Thus the governing equation for the transverse vibration of 

a flat orthotropic plate is given by 

D Ha 4 u/3 a J + 2(D 12 + 2D 66 )9 4 u/3a^3a2 + D 22 3 4 u/3a 4 

+Xu + phu = 0 , (A2.8) 
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where it should be understood that u=u^. And if the 
damping term is temporarily neglected, the equation is 

D n 3 4 u/3aJ + 2(D 12 +2D 66 )3 4 u/3c^3oi2 + D 22 3 4 u/3a 4 

+ phu = 0 . (A2.9) 

If the plate is simply supported on all four edges, the 
boundary conditions are zero displacement and zero bending 
moment along the edges. It can be shown that these two 
boundary conditions can be expressed as (see reference 56) 


at 

II 

0 

i — 1 

'V 

O 

II 

0 , 

(A2.10) 

and 

D n3 2 u/3a^ + 

D 12 3 2 u/3a 2 = 0 , 

(A2.ll) 

and at 

a 2 = °’ d 2 

u = 0 , 

(A2.12) 

and 

D 12 d 2 u/3a^ + 

D 22 3 2 u/3a 2 = 0 . 

(A2.13) 


It is well known in the theory of vibration that the general 
solution of the homogeneous partial differential equation is 
given by a series expansion in the natural modes of the 
structure, i.e. 

00 

u(a lf a 2 ,t) = ^ ( t ) U i (a 1 ,a 2 ) , (A2.14) 

i= i 

where n^(t) are the modal participation factors, and 
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U i^ a l ,a 2'* are t * ie e ^-9 en f unctions or modeshapes. 

In the case of the rectangular, simply supported, ortho- 
tropic flat plate, the eigenfunctions are given by 

U i^ a l ,a 2^ = U mn = sin ( m ™]_/ d ]_) sin ( nira 2/ d 2^ * 

m, n = 1 , 2 , 3 , . . . 00 , (A2.15) 

and the modal participation factors are of the form 


T1 ± C t ) = 


j (to t+0 ) 

J mn mn 


'mn 


(A2.16) 


where n are the influence coefficients of the modes 

mn 

in units of displacement (determined by the initial 
conditions), w are the radian natural frequencies of 
the modes, and © mn are the phase angles associated 
with the modes (determined by the initial conditions). 
Thus, the general solution to the homogeneous partial 
differential equation is given by 


u(a 1 ,a 2 ,t) = 


00 00 j(<D t+6 ) 

1 I n mn sin(mTTa^/d^)sin(nTT0!2/d2)e mn mn . (A2.17) 

m = i n= i 


And it is easily shown that this solution satisfies all of 
the simple support boundary conditions. Substituting this 
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solution into the governing partial differential equation of 
motion, one obtains 


00 00 

1 5 

m= i n= i 


D^Cmir/dj.) + 


* D 22 (nir/a 2 ) 4 - v 1 ph 


ifo) t+0 ) 

JV mn mn ; A 
e = u . 


2(D 12 +2D 66 ) (m 1 r/d 1 ) 2 (mr/d 2 ) 2 
n mn sin(mira 1 /d 1 )sin(mra 2 /d 2 ) 


(A2.18) 


If the series is zero for arbitrary values of n and 

mn 

0 , then each coefficient in the series must be 

mn 

identically zero. Thus 


D 11 (mTt/d 1 ) 4 + 2(D 12 + 2D 66 )(mir/d 1 ) 2 (m7/d 2 ) 2 


+ °22 (nil/d 2 )4 ' u mn ph = 0 * 


(A2.19) 


Rearranging this last equation algebraically, 


2 4 / . 

a) = it /ph 
mn 


D^Cm/d ^ 4 


+ 2(D 12 +2D 66 )(m/d 1 ) 2 (n/d 2 ) 2 + D 22 (n/d 2 ) 4 


(A2.20) 


where the a) are the radian natural frequencies of 
mn 


the plate. 


C^_ FREE RESPONSE FOR THE DAMPED CASE 

Now returning to the governing equation and reinstating 
the damping term, 
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(A2.21) 


D 11 3 4u / 3 a J + 2(D 12 + 2D g6 ) 3 4 u/3a 2 3a 2 + D 22 3 4 u/3a 4 
+Xu + phu = 0 . 

Now define the modal damping coefficient as 

5 mn = X/(20h %n ) . (A2.22) 

Substituting into the governing equation, one obtains 

D n3 4 u/3aJ + 2(D 12 +2D 66 )3 4 u/3a^3a 2 + D 22 3 4 u/3o 4 

+ ph2£ u + phu = 0 . (A2.23) 

mn mn 

We now assume a solution of the form 
u(a 1 # a 2 ,t) = 

00 oo yt+je 

1 I Tl mn sin(mna^/d^)sin(mra 2 /d 2 ) e mn , (A2.24) 

m = i n= i * * 


where y may be complex. Substituting into the governing 
equation , 


OO 00 

m= i n=i 


D 11 (mir/d 1 ) 4 + 2 (D 12 + 2D g6 ) (mtr/dj^ ) 2 (mr/d 2 ) 2 


+ D 22 (mr/d 2 ) 4 + Ph 2 C mn 0 ) mn y + phy 2 


ri mn sin (mir/d^ ) sin (nir/d 2 ) = 0 . 


(A2.25) 
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Now recall from section B that 


ph “mn = D l 1 tm l'/d 1 J 4 


+ 2(D 12 +2D 6g ) (mTT/d 1 ) 2 (mr/d 2 ) 2 + D 22 (mr/d 2 ) 4 


(A2. 26 ) 


Substituting, the result is 


00 oo 

1 l 

m= i n= i 


2 2 
pho3_„ + ph2E y + phy 

mn mn mn 1 1 


ri mn sin(m-rra 1 /d 1 )sin(mTa 2 /d 2 ) = 0 . 


(A2.27) 


Therefore, for arbitrary influence coefficients ri , 

mn 

it must be that 


phu3 2 + ph2f tjj y + phy 2 = 0 , 

mn mn mn 1 1 ’ 


(A2. 28) 


or 


y 2 + 21 a) y + 03 2 = 0 

1 ^mn mn r mn 


(A2.29) 


Solving this quadratic, 


y = - £ 03 - 

1 ^mn mn 


_2 2 2 

r 03 - 03 

^mn mn mn 


1/2 


(A2. 30) 


or 


Y = 


- F oj - -io) 
mn mn J mn 


1-5 


mn 


1/2 


(A2. 31) 
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Now define the damped radian natural frequencies as 

1/2 


ft =0) 
mn mn 


1-C 


mn 


(A2. 32) 


The general solution to the damped free vibration problem 
will then be 


u(a 1 ,a 2 ,t) = 


I I n mn sin(mira 1 /d^)sin(mra 2 /d 2 ) 


m= i n= i 


-F a) t t+0 ) 

s mn mn J mn mn J 
e e 


(A2.33) 
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Appendix III 

FORCED VIBRATIONAL RESPONSE OF A SIMPLY SUPPORTED 
RECTANGULAR ORTHOTROPIC PLATE 

A^ GOVERNING DIFFERENTIAL EQUATION 

The general forced vibrational response of a simply 
supported, rectangular orthotropic plate is governed by the 
non- homogeneous partial differential equation 

D^3 4 u/3aJ + 2(D 12 + 2D g6 ) 3 4 u/3a^3a \ + D 22 3 4 u/3a 4 

+ Xu + phu = q(a^,a 2 ,t) . (A3.1) 

B_j_ DEFINITION OF THE FORCING FUNCTION 

For the purposes of this study, the forcing function 
shall be defined as a normally incident (90° angle of 
incidence) acoustic plane wave impinging on the surface of 
the plate (refered to as the airborne component) combined 
simultaneously with a point vibrational input (referred to 
as the structureborne component). For this special case, 
the acoustic (airborne) input will be modeled as a simple 
harmonic forcing function, uniformly distributed over the 
surface of the plate and the vibrational (structureborne) 
input will be modeled as a simple harmonic forcing function, 
located at a discrete point on the plate. Thus, the right 
side of the governing equation is 
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q(a 1 ,a 2 ,t) = 

r -j<L “j<f> I 3 wt 

f g e 5(a 1 -c 1 )6(a 2 -c 2 ) + e e 

where 

f = the magnitude of the structureborne input 
(units of force), 

<j> = the relative phase of the structureborne input 

f = the magnitude of the airborne input (units of 

force per area), 

<j> a = the relative phase of the airborne input, 

c^ = the location of the point vibrational input 

in the coordinate direction, 

c 2 = the location of the point vibrational input 

in the coordinate direction, 

oj = the radian frequency of the simple harmonic 
forcing function, 

and the symbol 5 denotes the dirac delta function. 
Expanding this forcing function into a series 
representation , 

r -3*. _ ^a 1 

f s e 6(a 1 -c 1 )5(a 2 -c 2 ) + e 

oo ao 

1 q kl sin(kTTa 1 /d 1 )sin(liTa 2 /d 2 ) . 
k = i l=i x 

Multiplying the equation by sin(mTra^/d^) 


(A3. 2) 


(A3. 3) 
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• sin ( nTTo^/^ ) and integrating over the plate area 


?1 ?2 


0 0 


f e 
s 




5 C a^- c 1 )6(a 2 -c 2 ) + f a e 


-ik 


sin (miTa^/d^ ) sin (mra 2 /d 2 ) da^da 2 = 


00 00 

k l J 

k= l l=i 


*kl 


h ?2 


0 0 


sin(kTTa^/d 1 )sin(lTra 2 /d 2 ) 


sin (mTra^/d^^ ) sin (nTra 2 /d 2 ) da^da,, . 


(A3. 4) 


Integrating the left side of the equation and using the 
orthogonality principle on the right side of the equation, 

- j<t> 

q mn = 4f e /(d 1 d 2 ) sin (imrc 1 /d 1 ) sin( mrc 2 /d 2 ) 


-j<k 


+ 4f e a /(mmr^) ( 1-cos (mir )) (1-cos (nir) ) 

a. 


(A3. 5) 


Thus, the forcing function becomes 


q(a 1 ,a 2 ,t) = 


00 oo 

1 I 

m= i n= i 


4f e 
s 




/(d^d 2 ) sin (mire., /d., ) sin (nir^/d,, ) 




'2 2 


- j<l> o 

+ 4f e /(mnir ) ( 1-cos (mir ))( 1-cos ( nir ) ) 

3 


• sin(mna^/d^)sin(mra 2 /d 2 ) 


jwt 


(A3. 6) 
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C. STEADY STATE RESPONSE TO THE FORCING FUNCTION 


Now assume a solution of the form 


u( <x 1 , a 2 , t) 


00 oo 


j(u)t-e ) 

2 2 n mn sin(mira^/d^)sin(mra 2 /d 2 ) e mn 

m = i n= i 


(A3. 7) 


where the transient (complementary) solution has been 
neglected, since it exponentially damps out to zero, and 
only the steady state (particular) solution is considered. 
Substituting this solution into the governing equation and 
equating the coefficients of the resultant series, 


- j© 

(phu) +ph 25 a) (ja))-phu) 2 ) n e mn 
mn mn mn J mn 


- j<f) 

4f„ e 3 /(dj^) sin (mire., /d., ) sin (mrc 0 /d,, ) 


V "*1 


'2 2 


-j<|> o 

+ 4f e /(mnir ) ( 1-cos (mir )) (1-cos (rnr ) ) 


(A3. 8) 


Now utilizing the theory of complex variables, redefine the 
first term in brackets in the previous equation as 


(phw mn +ph2£ mn ( ju))-plW 


) = 


phw 


mn 


(i-wWr 

mn 


( 2 £ co /u) ) 

v s mn mn 


11/2 


D Y 


mn 


(A3. 9) 
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(A3. 10) 


where the phase angle v is defined to be 

3 mn 


Y = tan 
' mn 


-1 


( 2 W ,/ “mn )/(1 -“ 2/u mn ) 


Dividing through by this term, the equation for the 
influence coefficients is seen to be 


-je. 


^mn e 


mn 




s ■ JT a 

s e + a e 

mn mn 


1 -3Y 


mn 


where 


s mn = 4f s //( ‘ d l d 2 phu) mn' 1 sin (mirc^/c^ ) sin ( mrc 2 /d 2 ) 


/ 


(l-0) 2 /w 2 ) 2 + (2E to /to ) 2 
k / mn ; 1 s mn ' mn J 


1 1/2 


and 


a = 4f /(mnir phco ) 
mn a K mn 


(1-cos (miT ) ) (1-cos (nit ) ) 


/ 


(l-^ 2 //) 2 * (25 mnM / M]nn ) 2 


1 1/2 


Thus, the steady state solution of the problem is 


00 CO 


u( a l , a 2 » t ) =22 


m = i n= i 


s mn e + a m e 

mn mn 


sin(mira^/d^)sin(mra 2 /d 2 ) e 


-j<L 

+ a 

i 


modal 


(A3. 11) 


(A3. 12) 


(A3. 13) 


(A3. 14) 
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EL TRANSVERSE VIBRATIONAL VELOCITY OF THE PLATE 
The transverse vibrational velocity of the plate is 
by the relation 

v(a 1 ,a 2 ,t) = d/ dt {u(a^ , a , t) } , 

thus 


v(a 1 ,a 2 ,t) = 


00 00 
I I 

m=i n = 1 



-34>s 

-3^ 1 

jw 

s mn e 

+ a mn e 


j (wt-Y mn ) 

• sin(m-rra^/d^)sin(nira 2 /d 2 ) e 


Rearranging this last equation, and recalling that 

U = sin (imra, /d, ) sin(nira 0 /d„ ) , 

mn ll zz 

the equation for the velocity becomes 
v(a 1 ,a 2 ,t) = 


00 oo 

5 I 

m = i n= i 


U) 


-j*s , ~^a 

s mn e a mn e 


U e 
mn 


j( M t-Y n , n * 1 ' /2 ) 


00 oo 


- 5 I 

m = i n= i 


-j0 -j (o)t+ir/2 ) 

J mn „ J 
u) n e U e 

mn mn 


defined 


(A3. 15) 


(A3. 16) 


(A3. 17) 


(A3. 18) 
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Appendix IV 

POWER FLOW FOR COMBINED AIRBORNE AND STRUCTUREBORNE INPUTS 
TO A SIMPLY SUPPORTED RECTANGULAR ORTHOTROPIC PLATE 


Aj_ SPACE-TIME AVERAGED MEAN SQUARE SURFACE VELOCITY 

In Appendix III, it was found that the velocity of the 
plate is given by 


v(a 1 ,a 2 ,t) = 


00 00 


5 5 

m= i n= i 


con 


mn 



U 


mn 


j (cot+ir/2 ) 
e 


(A4.1) 


Now from the theory of complex variables, the time averaged 
mean square value of the surface velocity, (which is simple 
harmonic in character) is given by 


<v 2 ( 


a l» a 2 ^ > t 


1/2 Re{v (a. 


,a 2 ,t) 


v(a. 


> <*2 , t ) } 


= 1/2 


5 J 


+ je 


con 


kl 


L k= i l=i 


kl 


U kl ^ a l ’ 


a 2> 


- j (ajt + TT/2 ) 

e 


j© 


coin 


mn 


l m = i n= i 


mn 


U (a, , O 
mn 1 2 


+ j (oot+Tr/2 ) 
e 


00 oo 


= 1/2 


+ je 


1 5 wn kl e u kl ( ' a l ,a 2^ 


k= i l=i 


00 00 


X V /^mn 


l m = i n= i 


“ n « e lJ mn (a l’"2 ) 


(A4. 2 ) 
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(A4.3) 


00 GO 00 OO 


< V 


w /2 J ,1 I I 

k= i 1=1 m = i n = i 


+ j© k1 -j® mn 

n kl e T1 — e O^Ca^a^J-U^Ca^a,) 


'mn 


kl 1* 2 J ' J mn vu l , "2 


Now finding the space averaged mean square surface velocity 


2 

<v C a 1 , a 2 , t)> t = 


?2 ?1 


0 0 


<v ( , (^2 > t) > t da^dc^ 


00 OO 00 00 


= *\1 5 I 1 

k=i l=i m=i n=i 


+ j 0 ki -j 0 

^kl e n mn e 


mn 


1/S 


d 2 d l 


0 0 


U kl (a l» c, 2 ) * U inn (<l l» a 2 ) da l da 2 


Now recalling the definitions of U, , and U , 
3 kl mn 


d d 


o o 


D kl t “l>"2 ) 0 mn ( V a 2 ) d0, l da 2 


(A4. 4) 


?1 


sin (kira^/d^ ) sin (nura^/d^ ) da^ 


h 


sin (Ittc^A^ ) sinCmro^A^ ) da 2 • 


(A4. 5) 
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Utilizing the principle of orthogonality, it can be seen 
that the integral over the eigenfunctions is identically 
zero in all cases except the single case where k=m and l=n. 
Therefore, all of the terms in the series expansion are zero 
except those terms where k=m and l=n. Thus, the series over 
k, 1, m, and n, is simplified to a double series summing 
over m and n as follows: 


< v (VV^r.t 


00 00 


(JL) 


! /2 l l 


m= i n= i 


+10 -i© 

J mn J mn , /_ 

T1 e „ 6 * 1/S 

mn mn 


2 

sin (imra^/d^) da 1 


2 

sin (mra 2 /d 2 ) da 2 


(A4.6) 


Evaluating the integrals and noting that S = d]_d 2 , 
<v 2 (a 1 ,a 1 ,t)> r>t = 


- 00 00 +10 - 1 © 

2 /0 v v J mn J mn 

(0/8) )ne ne 

L L 'mn mn 

m= i n= l 


(A 4 . 7 ) 


Now substituting from Appendix III for the values of the 
influence coefficients for the case of combined airborne and 
structureborne inputs, 


175 



<v 2 > 


r, t 


00 00 


03 


/® 5 1 


m= i n= i 

-jf 




j<f> ; 


s e s + a e a 
mn mn 


DY 


mn 


•s ~^a 

s e + a e 

mn mn 


- DY 


mn 


(A4. 8) 


Multiplying out terms and simplifying, 


<v 2 > = uj 2 /8 

r, t 


00 00 
m= i n= i 


s 2 + a 2 + 2s a cosfi -4> ) 
mn mn mn mn ^ Y s T a 


(A4. 9) 


Now define the following terms 


00 00 


<v 


>r,t* I 5 

J m= i n= i 


mn 


(A4.10) 


00 00 


<V‘ 


>r,f " 2 * 5 5 

* m= i n= i 


mn 


(A4. 11) 


<v 2 > 

sxa 


00 00 


, = 03/8 y y 2s a cos(<b -d> ) 

r,t L L mn mn vr s T a 

9 m = i n= i 


(A4. 12) 


Thus, it can be easily seen that 


2 2 2 2 

<v> , = <V > . + <V > , + < V > . 

r,t sr,t a r,t sxa r,t 


(A4. 13) 
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B. POWER INPUT 


Recall from Appendix III that the forcing function is 
of the form 


q(a 1 ,a 2 ,t) = 


f s 6 




*5(o(2.'~c^)6(o( 2 — c 2 ) 


+ f a e 




(A4.14) 


The velocity produced by this input was given by 

CD CD 

v(a 1 ,a 2 ,t) = 55 


0) 


m = i n= i 


-j<|> -j<|> 

s e + a e 

mn mn 


U e 
mn 


j(wt-Y mn + ir/ 2 ) 


(A4. 15) 


Now the total time averaged power input to the plate is 
given by 


? 2?1 


<W(t)> t = < 


F • V da^da 2 >£ 


0 0 


d~ d n 


lim 1/T 

T+CO 




r ** 

Re{q(r, t) }Re{v(r , t) } do^do^ 

0 

0 

\ 0 

i 


dt 


?2 , d l 


0 0 


1/2 Re{q (r,t)v(r,t)} da-^do^ . 


(A4.16) 
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Substituting the appropriate expressions 
<w(t)> t = 



2 

i 

+ j<L + j<t> 1 

< 


1/2 Re{ 

f s e 5(a l“ c l )5(a 2- c 2 ) + f a e 

0 

6 

1 



-jut * « 

■® 2 2 “ 

m = 1 n= 1 


-j<t> I 

s e + a e U 

mn mn mn 


J(“ t -W T/2) . . . „ 

e } da 1 da 2 > fc . 


(A4. 17) 


Multiplying this out 


<W(t)>. 


d 2 d l 


0 0 


00 00 


1/2 Re { J j f ( 


■j<t> 


m = i n= i 


^ a l -C l^ a 2 -C 2^ “ 


-j<t>s -3*a 

s mn e + a mn e 


j(»t- Y mn + "/2) 

U mn e }d«, 1 do, 2 > t 


?2 ?1 


+ < 


00 00 


1/2 Re{ J j f. 




0 ) 


0 0 


m = i n= i 




s e 
mn 




+ e 
mn 


„ j (ut-y +it/2 ) 

U mn e > d “l da 2 > t 


(A4. 18) 
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Expanding terms and bringing the integration inside, 


<W(t)>, 


00 00 


• 1/2 Re{ j j £ 


0 ) 


m= i n= i 


s + a e 
mn mn 




mn 


?2 ?1 


stoj-opsto^Cj) u mn d« 1 ao 2 )> t 


0 0 

1/2 Re{ j j f 


00 00 


+ < 


0 ) 


m= i n= i 


s mn e + a mn 




d 2 d l 

r ^ <• x 


U„ da, da 0 }>. . 

mn 1 Z t 


(A4. 19) 


0 0 

Evaluating the integrals 


<w(t)> t = 


00 00 


:l/2 Re{ J ]> f 


(0 


m= i n= i 


s + a e 
mn mn 


-3<W 




CO CD 


sinCmirc^/d^sinCnTTCj/dj) }> fc + <1/2 Re{ ^ ^ f 


03 


m= i n= i 


-jtV't’s' . 

s e + a 

mn mn 




( 1 - cos ( mir ) ) / ( mir /d^ ) 


( 1-cos (mr ) /(nir/d 2 ) 


}>. 


(A4. 20) 
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Now taking the real part and noting that there is no longer 
any time dependent terms left in the equations, 


<w(t)>. 


00 00 


f s w/2 l l 


m= 1 n = i 


s cos ( Tr/2-y mr , ) + a m cos (ir/2-y_ +d> -<b ) 
mn ' mn mn 'mn T a T s 


sin(m'irc^/d^)sin(mrc2/d2 ) 


00 00 


+ f a s " /2 I I 


m = i n= l 


s cos ( ir/2-Y_+(j> -<b ) + a cos (ir/2-y ) 

mn mn T s T a mn mn 


( 1-cos (imr ) ) (1-cos (n-nr ) ) /(mnu ) 


(A4.21) 


Now note that in the special case of zero airborne input 


CO oo 


<i s (t)> t * £ s" /2 I l s mn cos( ’ /2 - Y m n ) 


m= i n= i 


sin(imrc^/d^)sin(nTTC2/d2) • 


(A4. 22) 


And in the special case of zero structureborne input 


00 00 


< w a (t) > t = f a S«/2 l l a mn cos (^/^-y^) 


m= i n = i 


( 1-cos (mir ))( 1-cos (mi ) )/(mmr ) . 


(A4. 23) 
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Note that the power input for the combined airborne and 
structureborne forcing function is not equal to the sum of 
the powers input for the airborne and structureborne cases 
acting alone. Or, in terms of the equations, 

<W(t)> fc / <W s (t)> t + <W a (t)> t . ( A4 . 24 ) 


C. POWER DISSIPATED 

For simple harmonic motion, the time averaged power 
dissipated is given by (reference 56, pg. 331) 


d 2?l 


<W(t)>. = < 


0 0 


Re{F(r , t) } *Re{v(r , t) } da 1 da 2 > fc 


= lim 1/T 
T+oo 


T 


.2 

r 1 


. j 

j 

Re{F v (r,t) }Re{v(r, t) } do^do^ 

0 

0 

i 0 

1 


dt 


= 1/2 Re{F* (r , t) v(r, t) } da^do^ (A4.25) 

0 0 

where F is the resistive force of dissipation per unit 
area due to structural and fluid damping effects, and F v 
is the component of this force in the direction of the 
transverse plate velocity. From the partial differential 
equation governing the plate's dynamics, given in appendix 
III, the component of this dissipative force in the direc- 
tion of the transverse velocity is given by 
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(A4. 26) 


F v (r,t) = Au(r,t) = Av(r,t) . 
Substituting, the result is 
<w(t)> t = 


?2 ?! 


0 0 
thus 


<w(t)> t = 


* ,-*• 


1/2 Re{Av (r,t)v(r,t)} dc^do^, 


? 2?1 


1/2 Re{v (r,t)v(r,t)} do^dc^ 


0 0 


= X 


?2 ?1 


0 0 


(A4. 27) 


2 4. 2 -v 

<v (r,t)> t da^da 2 = XS<v (r,t)> r t . 


(A4. 28) 


Substituting from Appendix IV section A. for <v > . 

r , u 

<W(t)>^ = oj^/8 • 


00 00 

5 2 ^ 


m = 1 n= 1 


2 2 

s + a + 2s a cos(<b -<b ) 
mn mn mn mn T s “a 


(A4. 29) 
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For the special case of zero airborne input 

OO CO 

<w (t)> t = 0) 2 /8 5 5 \S S 2 n . (A4.30) 

m= i n= 1 “ 1 

And for the special case of zero structureborne input 

OO OO 

<w (t)> t = 0) 2 /8 j 5 IS a 2 n . (A4.31) 

m = i n= i 

So it is again clear that the power dissipated for the case 
of the combined airborne and structureborne inputs is not 
equal to the sum of the powers dissipated for the cases of 
the structureborne and airborne inputs acting alone. Or, in 
terms of the equations, 

<W( t) > t / <W s (t)> t + <W a (t)> t . (A4.32) 
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Appendix V 

SOUND RADIATION OF A POINT SOURCE 
MONOPOLE SOUND RADIATION 

The ideal fluid potential function, <j>, for a simple 
harmonic point vibrating source (monopole) is well known in 
acoustics theory (see reference 59, pp. 153-171, reference 
60, pp. 153-183, or reference 61, pp. 111-118) and is given 
by the relation 

j (wt-kr ' ) 

5 = Q/(4irr' ) e (A5.1) 

where 

Q = the volumetric flowrate of the fluid, 
r' = the distance to the point of observation, 

a) = the radian driving frequency of the source, 

c q = the speed of sound in air, 

k = w / c 0 is the wavenumber, 

t = time . 

This potential function satisfies the wave equation 

jut 

V 2 5 - 1/c 2 9 2 $/3t 2 = Q e S(r') . (A5.2) 

The corresponding pressure and particle velocity produced by 
this point source are found from the equations 

p(r',t) = p Q 3<|/3t = j(jjp Q J > (A5.3) 

v(r,t) = -V 5 = -9i/3r' e r , = (jkr'+l)/r' J ® r i • (A5.4) 
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B. MONO POLE NEAR A REFLECTING SURFACE 


Now if the point source is placed a small distance h 
away from a perfectly rigid planar surface which extends to 
infinity in all directions as shown below, 



Figure (AV-1) -- Monopole near a reflecting surface. 


and given that the distance to the point of observation 
r'>>h, so that the distance from the point of observation to 
the point source and its image are both approximately r', it 
can be shown that the potential function for this combina- 
tion of monopole and image monopole is given by 
j (wt-kr ' ) 

5 = Q/(2irr' ) e cos (khcos (©' ) ) • (A5.5) 
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Now in the limit as h-*-0, this potential function becomes 
j (wt-kr ' ) 

5 = Q/(2irr' ) e . (A5.6) 

And the corresponding expressions for the pressure and the 
particle velocity are given by 

j (wt-kr ' ) 

p(r',t) = p Q 3<|/3t = jp Q w Q/(2irr' ) e 

= jP Q w 5 > (A5.7) 

. . j(wt- kr' ) 

v(r',t) = -V I = (jkr'+D/r' Q/(2irr' ) e e r , 

= ( jkr+l)/r' 5 e r , . (A5.8) 

Since, in this case, the monopole is simple harmonic in 
character, the time averaged acoustic intensity of this 
source is given by 

<l(r',t)> t = 1/2 Re{ p*(r',t) v(r',t) } . (A5.9) 

Now note that 
lim v(r ' , t) 

r 1 -voo 

j (wt-kr' ) 

= lim (jkr’+l)/r' Q/(2irr' ) e e , 

r 1 +°° 

j (wt-kr ' ) ^ ^ 

= jk Q/(2irr' ) e e r = jk e r « • (A5.10) 
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Multiplying by P q /p o 

lim v(r',t) = j (p o w)/(p o c o ) J e fl 
r 1 -*00 ^ 

= p(r' , t)/(p Q c o ) e r , . (A5.ll) 

Thus, in the acoustic far field, as r'-*-® 

v(r ' , t) = p(r' , t)/(p Q c o ) e r , , (A5.12) 

and substituting into the expression for intensity 
lim <1 (r ' , t) >. 

r ' ->oo t 

= 1/2 Re{ p* (r ' , t) p(r ' , t ) /(p Q c 0 ) } e r , . (A5.13) 

Thus, the magnitude of the far field time averaged acoustic 
intensity vector in the radially outward direction is given 
by 

<I r ,(r',t)> t = |p(r' ,t) | 2 /(2p o c q ) . (A5.14) 
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Appendix VI 

SOUND RADIATION FROM A SIMPLY SUPPORTED 
RECTANGULAR FLAT ORTHOTROPIC PLATE 

A;_ GEOMETRIC CONSIDERATIONS 

Consider the following geometry: 



Figure (AVI-1) -- Geometry of the simply supported plate. 


Note the coordinate system has its origin at the 
corner of the plate. The coordinate system has its 
origin at the center of the plate. 
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B^_ DERIVATION OF THE EQUATION FOR PRESSURE 

Utilizing Huygens principle of superposition of simple 

sources along with the equations developed in Appendix V, 

the sound generation of each incremental area, dS, of the 

plate is modeled as a monopole lying on a rigid planar 

reflecting surface. Recall that the incremental pressure in 

this case was given by 

^ j(o)t-kr') 

dp(r',t) = ju)P o Q/(2Trr' ) e 

-jkr' ju)t 

= jujp Q / C 2 tt r * ) e Q e . (A6.1) 

Recalling that Q represents the volumetric flowrate, it is 
seen that for simple harmonic vibration of the plate 
jwt 

Q e = v(B 1 ,6 2 >t) dB 1 dB 2 (A6.2) 

where v(B^>8 2 ,t) is the simple harmonic transverse 
plate velocity. Thus 

dp(r',t) = 

- jkr ' 

jcop o /(2nr') e v(8 1 ,6 2 ,t) d8 1 d8 2 . (A6.3) 

In the 8^ system, the coordinates of the point of 
observation are (8^,B 2 ,8 3 ) = ( r i> r 2 ,r 3^» 
and the coordinates of the incremental plate area are 
(8 1 ,8 2 ,8 3 ) = (8 1 ,8 2 ,0). Thus 
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(A6. 4) 


r ~ r l + r 2 + r 3 ’ 


and 


r' 2 = * tr 2- 6 2> 


2 2 
+ r 3 


- ' 2 - 2r l 8 l - 2r 2 B 2 * B l + S 2 ‘ 


(A6. 5) 


And from geometry 


r^ = r sin (e)cos (<j>) , 
r„ = r sin(e)sin(<t>) 


(A6.6) 
(A6. 7) 


Substituting this 


r ' = r 


1-2 sin(0)cos(<f>) C B-j^/r ) - 2 sin(e) sin(<f>) ( 6 2 /r ) 
1/2 


+ (S^r) 2 + (B 2 /r) 2 


(A6. 8) 


Now if we use the far field approximations 


1/ r 1 ^ 1/r << 1 , 


(A6. 9) 


so that 


r >> d^ and r >> d ^ 


(A6. 10) 


then 


I ~ 


1 - sin (0)cos (<}>) C S^/r ) - sin(0)sin((|)) ( & 2 /r ) 


. (A6.ll) 
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Substituting into the expression for the incremental 
pressure caused by the incremental plate area 


dp(r,t) = jwp o /(2iTr) 

- jkr(l-sin(e)cos(<{>) (e 1 /r)-sin(e)sin(<f>) (8 2 /r) ) 

• e 

• v(8 1 ,B 2 ,t) dB 1 dB 2 • 

Integrating over the entire area of the plate 

P(r,t) = 

d 2 /2 d x /2 

jwp Q / ( 2trr ) 

-d 2 /2 -d x /2 

- jkr ( 1-sin (0) cos (<fj) ( B^/r )-sin(0 ) sin(<j)) ( B 2 /r ) ) 

• e 


• vCB 1 ,B 2 ,t) 


dB 1 dB 2 . 


(A6. 12) 


- jkr 

= jwp / ( 2trr ) e 


d 2 /2 d 1 /2 

-d 2 /2 -d x /2 

• v(B 1 ,B 1 ,t) 


r jk sin(0)( B^os (<j>) +B 2 sin(<f>) ) 
e 


dB^dB^ . 


(A6.13) 
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Now make the following change of variables 


B 1 = a l - V 2 ’ B 2 = a 2 ■ d 2 /2 


dB-L = da^ , 


dB 2 = da 2 . 


(A6.14) 

(A6.15) 


Substituting this change of varibles 

^ -j kr 

p(r,t) = jo)p o /(2irr) e 
d, d i 

f z r 1 r jk sin(©) (( d 1 /2 )cos (<))) + ( a 2 -d 2 /2 )sin(<f>) ) 

e 


0 0 

• v(a 1 ,a 2 ,t) 
thus 


da^da 2 , 


(A6.16) 


P(r,t) = ja)p o /(2irr) 


-jk(r + d^/2 sin(e)cos (<j>) + d 2 /2 sin(e)sin(4>) ) 


?2 ?1 


0 0 


v(a 1 , a 2 , t) 


jk sin(e) ( a^cos (<|>) + a 2 sin(<f>) ) 


da^da 2 . 


Recall from Appendix III equation (A3. 16) 

jmt 

v(a 1 ,« 2> t) = jto e 


(A6. 17) 


00 oo 

1 1 n 

m = i n= i 


-je, 


mn 


mn 


sin(mira^/d^)sin(mra 2 /d 2 ) 


(A6. 18) 
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Substituting this into our expression for pressure 
p(r,t) = -0) 2 p o /(2iTr) 

j(wt-k(r + d^/2 sin(9)cos C 4> ) + 


r 2 


. J 

0 0 

i 


jk sin (0) ( c^cos (<t>) + a 2 sin (<)>)) 


OO 00 

Y Y n 

m = 1 n= i 


-je 


ntn 


mn 


sin (mira^/d^ ) sin (niro^A^ ) 


da^df^ • 


(A6.19) 


Interchanging summation and integration 
p(r,t) = -u 2 p o /(2nr) 

j(cot-k(r + d^/2 sin (0)cos (<j)) + ^2^ sin(0)sin(<()) ) ) 


00 00 

I 1 n 

m= i n= i 


-je 


mn 


mn 


d d 


0 0 


sin(mira^/d 1 )sin(mra2/d2 ) 


r jk sin(0) ( a^cos (<f>) +ai 2 sin (<J)) ) 


da 1 da 2 • (A6.20) 


Separating integrals and rewriting 
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p(r,t) = -co p Q /(2Trr) 


j(<ot-k(r + d^/2 sin(6)cos (<J>) + d^/2 sin (0) sin (<j>) ) ) 



oo 

00 

L 

m= 1 

1 

n= i 


h 

jk 

• 


e 

0 

1 



d 2 

jk 

• 


e 


-j© 


n mn e 


mn 


jk sin(e)cos(<(>) a-^ 


sinCimra^/d^^) da 1 


jk sin(©)sin(<|>) a 2 


sin(mra 2 /d 2 ) da 2 


(A6.21) 


C _ EVALUATION OF THE INTEGRALS 

The resulting integrals found in the last section can 
be written as 
d. 


Z = 
m 


l l ^ 
e 


jk sin (0)cos (<}>) 


sin (mira^/d^ ) da^ 




cos(k sin (0)cos (<j)) a-^) sintimra^/d^) da 1 


, d i 


sin(k sin (©)cos C 4> D a^) sin (miTa^/d^ ) da^ . (A6.22) 


and 
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sin(mra 2 /d 2 ) da 2 


?2 


jk sin(©)sin(<|)) a. 


?2 


cos(k sin(©)sin(<j>) a 2 ) sin(nira 2 /d 2 ) da 2 


sin(k sin(©)sin(<f>) a 2 ) sin(mra 2 /d 2 ) da 2 . (A6.23) 


With the integrals represented in this form, they can be 
evaluated by making a trigonometric substitution followed by 
straightforward integration, or more simply by looking them 
up in tables of integrals. (For example see reference 62.) 
From reference 62, pp. 435-437, formulae 296, 316, 318 and 
319, the following formulae were obtained: 

sin (ax) cos (bx) dx 

= -(cos( (a-b)x)/(a-b) + cos ( ( a+b )x)/(a+b) ) /2 for a 2 4 b^ , 


= (l/2a sin^(ax)) sgn(a) for a 2 = b^ . 


(A6.24) 
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sin(ax)sin(bx) dx 


= ( sin ( ( a-b )x) /(a-b ) - sin( (a+b )x) /(a+b ) ) /2 for a 2 / b 2 , 

= (x/2 - sin(2ax)/(4a) ) sgn(ab) for a 2 = b 2 . (A6.25) 


Evaluating the real part of Z m : 


If (k sin(0)cos (<f>) ) 2 / (miT/dj^) 2 


R e{Z m } = 


sin(miTa^/d 1 )cos(k sin(©)cos (<j)) a^) da^ 


= (l-cos(miT-k sin(0)cos(<j))d 1 ) )/(2mi:/d 1 -2k sin (0)cos C <|> ) ) 

+ ( 1-cos (miT+k sin(©)cos(<f>)d 1 ) )/(2mir/d 1 +2k sin(0)cos (<f>) ) . 

(A6.26) 


And with several lines of algebra and trigonometric 
substitutions, the final result is 


Re(Z m } 


1 - cos (mir )cos (k sin(©)cos C<f))d^) 


• (mir/d^)/((mir/d^) 2 -(k sin (0 ) cos (<f>) ) 2 ) . 


(A6.27) 


If (k sin (0 ) cos ( cj>) ) 2 = (miT/d^) 2 


Re( V 


s in (mir a ^/d^) cos (k sin(0)cos (<}>) a^) da.^ 
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(A6.28) 


= d^/(2mir) sin 2 (mir) - sin 2 (0) = 0 . 

Evaluating the imaginary part of Z m : 

If (k sin (e)cos (<J>) / (mir/d^) 2 
d l 

Im{Z m ) = sin(mira^/d^)sin(k sin(6)cos(<j))a^) da^ 

0 

= (sin(mir-k sin(6)cos(<j>)d^)/(2mir/d^-2k sin(0)cos C <|> ) ) 
-sin(mir+k sin(0)cos(<f>)d^)/(2imr/d^-2k sin(0)cos (<j)) ) . 

(A6.29) 

And with several lines of algebra and trigonometric 
substitutions, the final result is 

Im{Z m ) = -cos (mir ) sin (k sin (0)cos (<j>) d^ ) 

• (mn/d^)/((mir/d^)2-(k sin (0)cos (<j>) ) 2 ) . (A6.30) 

If (k sin (0)cos (<j>) ) 2 = (mir/d-^) 2 
d ! 

Im{Z m } = sin(nnra^/d^)sin(k sin(0)cos(<J>) a^) da-^ 

0 

= d^/2 - d 1 /(4mir) sin(2imr) sgn(sin(0)cos C <j> ) ) 

= d-^/2 sgn(sin(0)cos((|>) ) . (A6.31) 
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Summarizing the evaluation of Z m : 

If (k sin (e)cos (<}>)) 2 / (mir/d^) 2 then 

jk sin(0)cos(<j))d^ - 

Z m = 1 - cos(nnr) e 

• (mir/d 1 )/((mTr/d 1 ) 2 -(k sin(0)cos (<}>) ) 2 ) . (A6.32) 

If (k sin (e)cos (<()) ) 2 = (mir/d^) 2 then 

Z m = jd^/2 sgn(sin(©)cos(<f>) ) . (A6.33) 

An exactly analogous analysis can be performed on the 
integral Z . The results are given by: 

If (k sin(0)sin((j>) ) 2 / (niT/d 2 ) 2 then 

jk sin(©)sin(<}>)d 2 ' 

Z n = 1 - cos(nir) e 

• C n tt /d 2 ) / C ( nir/d 2 ) ( k sin (©) sin( <}>) ) 2 ) . (A6.34) 

If (k sin(0)sin(<j>) ) 2 = (niT/d 2 ) 2 then 

Z n = jd 2 /2 sgn ( sin (© ) sin (<jj) ) . (A6.35) 
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D. CALCULATION OF SOUND POWER RADIATED 


Recalling the last equation in section B of this 
appendix, the equation for pressure is now given by 

p(r,t) = -o) 2 p o /(2irr) 

j(wt-k(r + d^/2 sin(©)cos (<}>) + d£/2 sin(0)sin(<j)) ) ) 


00 


• I 

m= 1 


00 


1 

n= i 


-j®, 


ran 


mn 



(A6.36) 


where the values of the complex numbers Z , and Z can 
be evaluated using the equations developed in section C of 
this appendix. From this last equation, the acoustic pres- 
sure at an arbitrary point in space can be calculated where 
it has been assumed that the point of observation is in the 
far acoustic field, i.e. r >> d^ and r >> d 2 . Now 
recall from Appendix V, section B, that in the far acoustic 
field, the time averaged magnitude of the acoustic intensity 
vector in the radially outward direction at a point in space 
for a simple harmonic, monopole source is given by 

<I r ,(r',t)> t = |p(r ' , t) | 2 /(2p Q c o ) . (A6.37) 

Likewise, for simple harmonic motion of a simply supported 
rectangular plate, the far field time averaged magnitude of 
the acoustic intensity vector in the radially outward direc- 
tion at a point in space is given by 
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(A6. 38) 


<I r (r,t)> t = |p(r, t) | 2 /(2p Q c o ) , 

where r >> and r >> d2 . 

Making the appropriate substitutions, the far field time 
averaged magnitude of the acoustic intensity vector in the 
radially outward direction at a point in space for the plate 
vibrating in simple harmonic motion is given by 

<I r (r,t)> t = p*(r, t)p(r, t)/(2p o c o ) 


= Cp o u) 2 /(2 1 rr)) 2 /(2p o c o ) 


00 00 


v v + ^ e kl * * 

) ) n k1 e *-*- z,_ z. 

L k=i 1=1 K1 


J k 1 


00 00 


} l n 


l m= i n= i 


e Z Z 

mn m n 


(A6.39) 


Rearranging the order of summation 

<I r (r,t) > t = (p u 2 /(2irr)) 2 /(2p c ) J JJ 1 

k = i l=i m = i n = i 


00 00 00 oo 


+ j©L 


-j 0 r 


kl J ^mn * * 

n, . e x T1 e Z, Z, z z 

'kl 'mn k 1 m n 


(A6. 40) 
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Now recall from Appendix I, section Q that the time averaged 
sound power is given by 


<TT(t)> t 


<I*n> t dS . 


(A6. 41) 


In the present case under consideration, the control surface 
to be integrated over is a hemispherical surface in the far 
acoustic field, and the normal unit vector is in the 
radially outward direction. Thus 
2ir it/2 


<1T(t)> t = 


< v t > > t 


sin(6) d©d<| > 


(A6.42) 


Substituting the expression for acoustic intensity, the time 
averaged sound power radiated by the plate is 


<TT ( t ) > t = 


2tt it/2 
0 0 


00 00 00 00 

(p o o) 2 /(2Tr)) 2 sin(0)/(2p c ) \ \ \ \ 

° ^ ° k = i 1= i m= i n= i 



e 


•j© 


kl 




dQd<j) . 


(A6. 43) 


201 



E. SOUND POWER FOR THE SPECIAL CASE OF COMBINED INPUTS 


Now recall from Appendix II, section C, that the 
equations for the influence coefficients are given by 

+ j©; 


n kl e 


'kl 


+ j<l>s + ^ a 

s kl e + a kl e 


+ ^ ki 
e , 


(A6. 44) 


and 


-j© 


ri e 
mn 


mn 


-j4> a 

s e + a e 

mn mn 


- DY 


mn 


( A6 . 45 ) 


Thus 


+ ^®kl ’^mn -^^kl~^mn^ 

n kl e n mn e = e 


s, , s + a, , a 
kl mn kl mn 


i(* s -V + KVV 

+ s. . a„ e + a. , s e 

kl mn kl mn 


(A6.46) 


Substituting this equation for the product of the influence 
coefficients into the equation for the sound power 

<TT( t ) > . = 


2tt tt/2 



1 


0 

6 

I 


OO OO OO 00 


(p o 03 2 /( 2ir) ) 2 sin(e)/(2p Q c o ) I I I 5 

k = i l=i m = i n=i 


j(<t> s -4> a ) j(4> a -<t> s ) 

s kl s mn + a kl a mn + s kl a mn e + a kl s mn e 


j ^ki - Y m n^ * * 

e Z, Z, Z Z 

k 1 m n 


d0dc}> . 


(A6. 47) 
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Now since the sound power is a real quantity (note that it 
is calculated from the conjugate product of the pressure 
times itself), it is known that the result of the series in 
the above equation is a real number. Using this knowledge, 
and some simple complex algebra, the equation for sound 
power can be written as: 


<TT(t)>. 


2 it it/2 


0 0 


(p o w 2 /(2ir)) 2 sin(e)/(2p Q c o ) 


00 00 

k S J 

k= l l=i 


K ' 2 ' 2 i ' 2 


00 oo 

+ Re{ l j 


m= i n = i (kl/mn) 


j(Y kl _Y mn ) * * 

s. . s e K1 Z.Z. Z Z } 

kl mn k 1 m n 


00 00 




a kl l* k l 2 l z ll 2 


00 CO 


V V j(Y kl“ Y mn ) * * 

+ Re { ) ) a. . a e K1 n Z. Z, Z Z } 

L i , kl mn klmn 

m= i n= i (kl/mn) 


00 00 


+ I I 2 s kl a kl cos l Z k^ ^ Z l^ 

k= i l=i L 


OO oo 

+ Re{ 5 I 


m= i n = i (kl/mn) 


s, . a e 
kl mn 


3(4> s -<l> a ) ^ Y kl _Y mn^ 


\ h 2 m 2 n > 


d6d<|) 


(A6. 48) 
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Now define the terms in equation (A6.48) as follows: 


<T Vt 


2 it tt/2 


0 0 


(p o 03 2 /(2tt)) 2 sin(6)/(2p o c o ) 


00 OO 


J A , 

k= i 1=1 L 


4i i\i 2 i^i 2 


00 00 




m= i n = i (kl/mn) 


j(Y k l"Y mn ) * * 

S., S m n e Z V Z m Z n > 

kl mn k 1 m n 


d©d<j> . 


(A6.49) 


<TT > . 
a t 


2ir tt/2 

(p o co 2 /(2tt)) 2 sin(0)/(2p o c o ) • 

0 0 


J. X 




2 


Re{ l l 


m= i n = i (kl/mn) 


^kl’W * * 

a. , a e Z,Z, Z Z } 

kl mn k 1 m n 


d0d<|> . 


(A6. 50) 
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<TT >. 

sxa t 


2tt tt/2 


0 0 


(p Q w /(2iOr sin(0)/(2p o c o ) 


OO GO 


+ J. A 2 s ki a ki cosC4> s -4> a ) |z k l 2 |z 1 i 

k= i 1= i L 


OO OO 


„ , \ \ . 3 C<t> s -<t> a 3 j(T kl -Y mn ) 

Re i 2 , . s kl a mn 

m= i n = i (kl/mnj 


* * 

Z, Z Z Z } 
k 1 m n 


d0d<(> . 


(A6. 51) 


Thus we see that the total sound power radiated can be 
written as 


<TT> , = <TT > , + <TT > , + <T T >. . 

t st at sxa t 


(A6.52) 
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Appendix VII 


CALCULATION OF THE BENDING RIGIDITY CONSTANTS 
FOR A SYMMETRIC LAMINATE 
WITH MULTIPLE GENERALLY ORTHOTROPIC LAYERS 


A;_ STRESS-STRAIN RELATIONS 

The stress-strain relations for plane stress in a 
single orthotropic layer (from reference 42, pg.46) is given 
by 


J 

H* 

J 

II 

Q ll 

Q 12 

i 

o 

Q 

to 

II 

°12 

Q 22 

0 

i — 1 

j 

II 

1 

o 

0 

°66 


2 

*12 


(A7.1) 


where the are the reduced stiffnesses and are given 

in terms of the Young's moduli, the shear modulus, and the 
Poisson's ratios of the orthotropic layer by the relations 

Q 11 = E 1 /U-P 12 1J 21^ ’ 

Q 12 = Pi2 E 2 / ^ 1 ^12 li 21 ) = lI 21 E l /(1 - li 12 U 21 ) 


q 22 = E 2 /(i-u 12 y 21 ), 


Q 66 = G 12 * 


(A7.2) 


For an orthotropic layer of arbitrary orientation 
(angle) the stress-strain relations (from reference 42, pg. 
51) are given by 
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(A7.3) 


a l 

II 

Qn 

^12 ^16 


e l 

a 2 


Ql2 

^22 ^26 


e 2 

T 12 

= 

Ql6 

^26 ^66 


Y 12 


where the transformed reduced stiffnesses Q. . are 

1 3 

given by the relations 


'll 


A O O >1 

Q 11 cos 4 © + 2 (Q^ 2 + 2Q 66 )sin z 0cos‘ i e + Q 22 sin 4 0, 


'12 


/ . . 2 2 , . 4 4 . 

(Qll+Q 2 2* 4 Q66 ^ sin 0COS 0 + Q 12 *- Sin 0+cos 0 J> 


'22 


A 0 0/1 

Q^sin^e + 2(Q 12 + 2Qgg)sin^0cos <: e + Q^cos^Q, 


'16 


( ^ 2 — 2Qg 6 ) ® i n0co s 0 + (Q^ 2 ~ Q 22 + 2Qgg ) sin ©cos©, 


'26 


3 3 

(Qll-Ql2-2Q 6 6^ sin 0COS0 + ^l2"^22 +2 ^66 ) sin9cos 9 > 


'66 


(Qn+Q 22 “ 2 Qi 2 - 2 Q 66 ) s i n20cos20 + Qgg (sin 4 0 +cos 4 0 ) . 


(A7.4) 


B. CALCULATION OF THE BENDING RIGIDITIES 

Once the transformed reduced stiffnesses Q. . are 

1 3 

known, the bending rigidities may be calculated 

(from reference 42, pp. 162-166) with the equation 



1/3 


N 



(Q. .), ( h l 

v ij'k v k 


h 


3 

k-1 



(A7.5) 
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where N is the total number of layers, h^ denotes the 
distance from the midsurface to the lower surface of the kth 
layer, and h^ ^ denotes the distance from midsurface 
to the upper surface of the kth layer. 
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Appendix VIII 

MEASUREMENT ERROR ANALYSIS 
A. DESCRIPTION OF THE ANALYSIS 


Recall from Chapter II equations (3.29) and 


<Tr a > t = <TT> t (<J a /o) (a - o s )/(a a - <x s ) , 


< " n "s > t ■ <TT> t (o s /o) (a a - o)/(o a - “s 5 


(3.31): 


(3.29) 

(3.31) 


These two equations represent the estimates of the 
airborne and structureborne sound power components of the 
measured combined sound power radiated by a structure. 

These estimates involve three independent measured quanti- 
ties, viz. the airborne radiation efficiency, a , the 
structureborne radiation efficiency, a g , and the radia- 
tion efficiency of the combined sound, a. Measurement 
errors in any of these three quantities will cause errors in 
the estimates of the airborne and structureborne components. 
The results of an error analysis which investigates the 
sensitivity of the estimates to each type of possible mea- 
surement errors are presented here. In some cases, the 
severity of the error is seen to be heavily dependent on the 
difference in the actual airborne and structureborne radia- 
tion efficiencies or on the actual radiation efficiency of 
the combined noise. Therefore, the results are presented as 
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a with 
s 


different plots for different values of - 
each plot showing a family of curves for different values of 
a. In each figure it is assumed that the actual value of 
a lies somewhere between the actual values of a g and 
a . Thus, each figure has a family of curves with each 
curve representing a different value of the fractional com- 
bined efficiency, a^, defined by the equation: 

a f = (c - a s ) /Ca a - a g ) (A8.1) 

where takes on values between 0 and 1 exclusively. 

It should be noted that while under most circumstances a 
has a value less that a. and greater than a _, this 
is not always the case. Furthermore, it can shown both 
analytically and experimentally that under relatively rare 
circumstances the value of a ^ can fall outside the range 
0 to 1 . 

B. ERRORS IN THE STRUCTUREBORNE ESTIMATE 

Figures (AVIII-1) through (AVIII-3) show the results of 

the error analysis on equation (3.31) for the errors 

expected in the structureborne sound power estimate, 

<TT >. , due to errors in the measurements, 
s t 

Figures (AVIII-la) through (AVIII-lc) show how errors 

in the measurement of a effects the structureborne 

s 

sound power estimate, <TT s > t , for the cases of a 5 
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dB, 15 dB, and 25 dB difference in the airborne and struc- 

tureborne radiation efficiencies (a - c ), respec- 

tively. These three figures show that the relationship 

between measurement errors of o and the estimate 

s 

<TT s > t are essentially linear and independent of the 
radiation efficiency of the combined noise (since all six 
curves in each of these figures fall almost exactly on top 
of one another). 

Figures (AVIII-2a) through (AVIII-2c) show how errors 

in the measurement of a effects the structureborne 

a 

sound power estimate, <TT >. , for the cases of a 5 

b U 

dB, 15 dB, and 25 dB difference in the airborne and struc- 
tureborne radiation efficiencies (a - a ), respec- 
tively. These three figures show that the relationship 

between measurement errors of a and the estimate 

a 

<TT >. is nearly independent of the value of (a 

S t 3 

- cf s ), but is extremely sensitive to the value of the 

radiation efficiency of the combined noise a. If the 

combined noise has a radiation efficiency nearly equal to 

that of the purely airborne input, then small errors in the 

measurement of a can cause large errors in the estimate 

a 

<TT > . 
s t 

Figures (AVIII-3a) through (AVIII-3c) show how errors 
in the measurement of the radiation efficiency of the com- 
bined noise, o, effect the structureborne sound power 
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estimate, <TT >^, for the cases of a 5 dB, 15 dB, and 
25 dB difference in the airborne and structureborne radia- 
tion efficiencies (a - a ), respectively. These 

a S 

three figures show results that are similar to the results 
found in figures (AVIII-2a) through (AVIII-2c) in that 
measurement errors in a produce errors in the estimate 
<TT >. that are nearly independent of the value of 

b U 

(a - a ), but are extremely sensitive to the actual 
value of the radiation efficiency of the combined noise. If 
the combined noise has a radiation efficiency nearly equal 
to that of the purely airborne 1 input, small errors in the 
measurement of a can produce large errors in the estimate 

<Vt- 

ERRORS IN THE AIRBORNE ESTIMATE 

Figures (AVIII-4) through (AVIII-6) show the results of 
the error analysis on equation (3.29) for the errors ex- 
pected in the airborne sound power estimate, <TT >. , 
due to errors in the measurements. 

Figures (AVIII-4a) through (AVIII-4c) show how errors 

in the measurement of a effects the airborne sound 

s 

power estimate, <TT > , for the cases of a 5 dB, 15 

cL U 

dB, and 25 dB difference in airborne and structureborne 

radiation efficiencies (a - a ), respectively. 

These three figures indicate that measurement errors in 

a cause sizable errors in the airborne estimate, 
s 
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<TT >. , only if the difference in the airborne and 

cL t 

structureborne radiation efficiencies (a - a ) is 

a S 

small (5 dB or less) and if the radiation efficiency of the 

combined noise, a is nearly equal to that of a purely 

structureborne input. (See figure (AVIII-4a)). As the 

difference in the airborne and structureborne radiation 

efficiencies (a -a ) increases, the airborne 

estimate, <TT >. , becomes increasingly insensitive to 
a u 

the value of the combined radiation efficiency, a, and 

eventually becomes completely independent of errors in the 

structureborne radiation efficiency a g altogether (as 

indicated in figures (AVIII-4b) and (AVIII-4c)). 

Figures (AVIII-5a) through (AVIII-5c) show how errors 

in the measurement of a effects the airborne estimate, 

a. 

<TT >. , for the cases of a 5 dB, 15 dB, and 25 dB 
a l. 

difference in the airborne and structureborne radiation 

efficiencies (a -a ) , respectively. Figure (AVIII- 
as 

5a) indicate that errors in the measurement of a cause 

a 

small errors in the estimate <TT >, when the differ- 

a t 

ence in the airborne and structureborne radiation efficien- 
cies (a - a ) is small (5 dB or less). The errors 
a s 

produced in the estimate <TT >, become increasingly 

a c 

smaller as (a - a ) becomes larger (as indicated by 
a s 

figures (AVIII-5b) and (AVIII-5c)). Furthermore, the errors 
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produced in the estimate are seen to be virtually independ- 
ent of the value of the combined radiation efficiency a. 

Figures (AVIII-6a) through (AVIII-6c) show how errors 
in the measurement of a effect the airborne sound power 
estimate, <TT >. , for the cases of a 5 dB, 15 dB, and 
25 dB difference in the airborne and structureborne effi- 
ciencies (a - a ), respectively. Figure (AVIII-6a) 
shows that the errors in the estimate can be quite large if 
the difference (a -a ) in small (5 dB or less) and 

3 S 

if the actual radiation efficiency of the combined noise, 
a, is nearly equal to the structureborne radiation effi- 
ciency, a . As the difference (a - a ) becomes 
s as 

larger, however, the errors in the estimate become less 
severe and are almost independent of the actual value of a 
(as indicated by figures (AVIII-6b) and (AVIII-6c)). 

Dj_ CONCLUSIONS 

Based on the results of the error analysis, the most 
severe type of error possible occurs when the combined sound 
power is strongly dominated by the airborne component (a 
is nearly equal to o ) and errors are made in the mea- 
surement of a or in the measurement of a . When this 

3 

happens, the estimate <TT s > t will always be highly in 
error . 

The second most severe type of error occurs when the 
combined sound power is strongly dominated by the 
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structureborne component ( a is nearly equal to a g ) 
and errors are made in the measurement of a or in the 
measurement of a g . When this happens, the estimate of 
the airborne sound power, <7T >, will be largely in 
error only if the difference in the airborne and structure- 
borne radiation efficiencies (a - a ) is small 

3 S 

(5 dB or less). 

All other errors in the estimates due to measurement 
errors are insignificant in comparison with the aforemen- 
tioned two most serious types. Furthermore, from these 
results, it can be seen that when applying these estimates 
to an arbitrary noise path identification problem, measure- 
ment errors are more likely to cause inaccuracies in the 
estimate of the structureborne sound power component than in 
the estimate of the airborne sound power component. 
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Appendix IX 
COMPUTER PROGRAMS 


A. PROGRAM MODAL 

The purpose of program MODAL was to calculate the nat- 
ural frequencies of a plate of either orthotropic or regular 
symmetric angle-ply laminate construction. Data files re- 
quired as input to this program include file DATPCR which 
contains the physical characteristics of the plate, file 
DATPSP which contains the plate stiffness parameters, and 
file DATNFL which contains the uppermost frequency limit for 
the calculation of the natural frequencies. Information 
contained in file DATPCR includes the length, width, thick- 
ness, and density of the plate. Information in file DATPSP 
tells the program whether the plate is orthotropic or a 
laminate and includes information such as the Young's 
moduli, the modulus of rigidity, Poisson's ratios, and the 
number of plies and their orientation if the plate is a 
laminate. File DATNFL tells the program where it can stop 
calculating natural frequencies by specifying an uppermost 
frequency above which no natural frequencies shall be calc- 
ulated . 

The calculation scheme of this program is as follows: 
The main program calls subroutines RPCHR, RSPAR , and RFLIM 
which read and subroutines WPCHR , WSPAR, and WFLIM which 
write files DATPCR, DATPSP, and DATNFL respectively. If the 
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plate is orthotropic, subroutine RSPAR calculates the bend- 
ing rigidity constants D^j using equations (A2.4) 
through (A2.7). If the plate is a laminate, subroutine 
RSPAR calls subroutine CBC which calculates the bending con- 
stants using the equations outlined in Appendix VII. Sub- 
routine CFREQ is then called to calculate the natural freq- 
uencies. The natural frequencies are calculated using equa- 
tion (3.6). The algorithm begins by setting mode numbers m 
and n equal to 1. Natural frequencies are then calculated 
with mode number n being incremented by 1 each time through 
the loop until the natural frequency calculated exceeds the 
uppermost frequency limit set by file DATNFL . Mode number m 
is then incremented by 1 and the entire calculation proc- 
edure is then repeated. Calculations cease when for some 
value of m, n = 1 and the calculated natural frequency ex- 
ceeds the uppermost frequency limit. With this calculation 
scheme, all flexural modeshapes which have natural frequen- 
cies less than or equal to the uppermost frequency limit are 
obtained. After the calulations are complete, program con- 
trol is returned to the main program and subroutine SORT is 
called which sorts out the modes in their order of ascending 
natural frequency. Subroutine WFREQ is then called to write 
this sorted list of mode numbers versus natural frequency in 
file DATMNF which is used in subsequent calculations. 
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B. PROGRAM SOUND 


The purpose of program SOUND was to calculate the sound 
power radiated by a plate and its radiation efficiency due 
to combined airborne and structureborne inputs. Data files 
required as input to this program include files DATPCR, 
DATPSP, DATNFL, and DATMNF (discussed in the previous sec- 
tion) and files DATFF and DATDMP which contain the forcing 
function and the damping model, respectively. Information 
contained in DATFF includes the magnitudes and relative 
phases of the airborne and structureborne inputs and the 
location of the point vibrational (structureborne) input. 
File DATDMP contains the two damping coefficients of equa- 
tion (3.11) or (3.12). 

The calculation scheme of this program is as follows: 
First, the program initializes the value of various cons- 
tants such as the characteristic acoustic impedance of the 
fluid medium, and the number of lines of resolution desired 
in the frequency domain, etc... Second, the program calls 
several subroutines which both read and write the informa- 
tion contained in the input files. The program checks the 
forcing functions and if possible will take advantage of 
symmetry in the numerical integration which is to be per- 
formed later on. The program then begins its outermost loop 
which sets the forcing frequency. Subroutines SBIC and ABIC 
are called to calculate the modal influence coefficients for 
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the structureborne and airborne inputs at the specified 
forcing frequency using equations (3.16) and (3.17). (The 
natural frequency information contained in file DATMNF, the 
damping information in DATDMP, and the forcing function 
information in DATFF is used to perform these calculations.) 
Subroutine STAMSSV is then called to calculate the space- 
time averaged mean square surface velocity of the plate due 
to the combined inputs at the specified frequency using 
equation (3.19). Finally, subroutine CAREDB is called to 
perform the numerical integration and calculate the time 
averaged sound power radiated by the plate due to the com- 
bined inputs at the specified forcing frequency using equa- 
tion (A6.43). At this point, control returns to the main 
program and the final results are calculated. The last step 
in the outer loop was to write the results including sound 
power level, radiation efficiency, mean-square surface velo- 
city, and two error checks on the numerical integration to 
the output file named DATA. The entire outer loop is then 
repeated for different values of the forcing frequency. In 
addition to file DATA, the various input parameters used in 
the program were stored on a separate file named BANNER. 

The numerical integration subroutine CAREDB utilizes an 
adaptive quadrature scheme to integrate over a two-dimen- 
sional hemispherical surface. This routine uses the trape- 
zoidal approximation with interval halfing and utilizes 
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Richardson extrapolation to improve its estimates. To use 
subroutine CAREDB, the user must provide a function subrou- 
tine containing the integrand. Function subroutine FTAI was 
provided to calculate the value of the integrand at an arbi- 
trary polar and azimuthal angle at a point on the hemispher- 
ical surface using equation (A6.43). Subroutine FTAI con- 
tains the innermost loop in the program and, according to 
equation (A6.43), must sum over all the modes of the plate 
to obtain the value of the integrand. Since modes whose 
natural frequencies are much, much greater than the forcing 
frequency are known to have a very small response, and since 
these modes are being driven in their stiffness controlled 
regions, their contribution to both the dynamic and acoustic 
responses of the plate are expected to be insignificant. 
Thus, when calculating the value of the integrand, function 
FTAI considers only those modes which have natural frequen- 
cies in the 0-1000 Hz frequency range plus any additional 
modes whose resonance frequencies are less than or equal to 
twice the forcing frequency. This method of calculating the 
integrand has virtually no effect on the accuracy of the 
calculations and reduces the necessary computation time sig- 
nificantly. 
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C. PROGRAM PREDICT 


The purpose of program PREDICT was to implement equa- 
tions (3.29) and (3.31) for the separation and prediction of 
the airborne and structureborne sound power components of 
some combined inputs case of interest. The code was written 
so that it could be used with the output data files (DATA) 
generated by the analytical computer model (program SOUND), 
and so that it could also be used with data files generated 
by the experimental study. 

The calculation scheme of PREDICT is as follows: The 
user is prompted for the names of the three data files con- 
taining (1) the structureborne radiation efficiency, (2) the 
airborne radiation efficiency, and (3) the combined sound 
power and combined radiation efficiency. The program then 
uses equation (3.29) to predict the airborne component and 
equation (3.31) to predict the structureborne component of 
the combined sound power. No attempt at predicting the com- 
ponents is made if the difference in airborne and structure- 
borne radiation efficiencies, (a -a ), is less than 

a s 

2 dB, or if either of the predicted sound powers is nega- 
tive. The indication that one of these two events has 
occured is signaled when the program returns a value of 0.0 
for either or both of the airborne and structureborne sound 
power predictions. The requirement that the difference in 
efficiencies (a -a ) be greater that 2 dB was 

3. S 


227 



selected on the basis that the experimental measurements of 
radiation efficiency are at best only accurate to plus or 
minus 1 dB. The requirement that the estimates be non- 
negative was dictated by the fact that one can not take the 
logrithm of a negative number (which is necessary on a dB 
scale), not because the estimates can not or should not ever 
be negative. In fact, both the analytical and experimental 
evidence suggests that the estimates will quite routinely 
have negative values. 
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The following pages are a listing of program MODAL 


PROGRAM MODAL ( INPUT, OUTPUT, TAPE5=INP0T, TAPE6=0UTPUT ) 
DIMENSION OMEGA ( 1000 ) , NMAX ( 50 ) 

DIMENSION MNSORT (1000 , 2 ) ,FREQ ( 1000 ) 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NU1 , NU2 
COMMON/LIMIT/FMAX 
COMMON/NF/MODES ,MMAX , NMAX , OMEGA 
COMMON/PCHR/D1 ,D2,H, RHO 

COMMON/ SPARA/Dll , D22 , D66 , D12 , D16 , D26 , DK , 

*E1 ,E2 ,G, NU1 , NU2 
COMMON/SORT/MNSORT , NOM , FREQ 
COMMON/TYPE/IPTYPE , NPLYS , IPO 
OPEN ( UNIT=1, FILE= ' DATPCR' , STATUS= ' OLD' ) 

OPEN ( UNIT=2 ,FILE= ' DATPSP ' , STATUS= ' OLD ' ) 

OPEN ( UNIT=3 ,FILE= ' DATNFL ' , STATUS= 1 OLD ' ) 

OPEN ( UNIT=7 , FILE= ' DATMNF ' , STATUS= ' NEW ' ) 

CALL RPC HR 
CALL WPCHR 
CALL RSPAR 
CALL WSPAR 

IF ( IPTYPE . NE . 0 . AND . IPTYPE . NE . 1 ) THEN 

GOTO 10 

END IF 

CALL RFLIM 

CALL WFLIM 

CALL CFREQ 

CALL SORT 

CALL WFREQ 

CLOSE ( 1 ) 

CLOSE ( 2 ) 

CLOSE ( 3 ) 

CLOSE ( 7 ) 

10 STOP 
END 


SUBROUTINE RPCHR 
COMMON/PCHR/D1 , D2 , H , RHO 
READ (1,10) Dl , D2 

10 FORMAT (//4X,E10.4/4X,E10. 4) 
READ (1, 20 )H, RHO 
20 FORMAT (3X,E10.4/5X,E10.4) 
RETURN 
END 
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o o on 


SUBROUTINE WPCHR 
C0MM0N/PCHR/D1 , D2 , H , RHO 
WRITE (6 ,10 ) 

10 FORMAT (IX, 14 HPANEL GEOMETRY/) 

WRITE ( 6 , 20 ) Dl , D2 

20 FORMATdX, 3HD1=,E10 . 4/1X, 3HD2=,E10 . 4 ) 
WRITE (6, 30 )H, RHO 

30 FORMAT ( IX , 2HH= , E10 . 4/1X , 4HRHO= , E10 . 4// ) 
RETURN 
END 


SUBROUTINE RSPAR 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NU1 , NU2 
COMMON/PCHR/D1 , D2 , H , RHO 
COMMON/SPARA/D11, D22 ,D66 ,Dl2,Dl6,D26 ,DK, 
*E1 ,E2 , G ,NUl , NU2 
COMMON/TYPE /IPTYPE , NPLYS , IPO 
READ (2,10) IPTYPE , NPLYS , IPO 
10 FORMAT (//7X, 11/8 X, 1 2/13 X, 814 ) 

READ ( 2 , 20 ) El , E2 , G 

20 FORMAT ( 4X, E10 . 4/4X , E10 . 4/4X , E10 . 4) 

READ (2,30) NUl , NU2 
30 FORMAT (5X,F7.3/5X,F7. 3) 

Dll=0 . 0 
D22=0.0 
D66=0 . 0 
D12=0 . 0 
D16=0 . 0 
D26=0.0 

IF ( IPTYPE. EQ.0) THEN 
CONST=H**3/12 .0 
Dll=CONST*El/ ( 1 . 0-NUl*NU2 ) 
D22=CONST*E2/(1.0-NU1*NU2 ) 

D66=CONST*G 

Dl2=CONST*El*NU2/( 1 . 0-NUl*NU2 ) 

DK=2 . 0* (D12+2 . 0*D66 ) 

END IF 

IF ( IPTYPE. EQ.l) THEN 

CALL CBC 

END IF 

RETURN 

END 


SUBROUTINE CBC 
INTEGER IPTYPE, NPLYS,IPO(8) 
REAL NUl , NU2 
COMMON/PCHR/D1 ,D2 ,H, RHO 
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COMMON/ SPARA/Dll , D22 , D66 , D12 , D16 , D26 , DK , 
*E1,E2,G,NU1,NU2 
COMMON/TYPE/IPTYPE , NPLYS , IPO 
PARAMETER (PI=3. 14159265359) 

Q11=E1/ ( 1 . 0-NU1 *NU2 ) 

Q22=E2/(1.0-NU1*NU2) 

Q12=NUl*E2/ ( 1 . 0-NUl*NU2 ) 

Q66=G 
N2=NPLYS/2 
DZ=H/FLOAT ( NPLYS ) 

DO 10 K=1,N2 

ZK=H/2 . O-FLOAT (K ) *DZ 

ZKM1=ZK+DZ 

A=FLOAT ( IPO ( K ) ) *PI/180 .0 
DUM=2 . 0/3 . 0* ( ZKM1**3-ZK**3 ) 

QBll=Qll*COS ( A) **4+ 

*2.0*(Q12+2.0*Q66 ) *SIN ( A) **2*COS(A) **2+ 
*Q22*SIN(A)**4 

QB12= (Q11+Q22-4 .0*Q66)*SIN(A) **2*COS ( A) **2+ 
*Q12*(SIN(A) ** 4+COS (A) **4) 

QB22=Q11*SIN( A) **4+ 

*2.0*(Q12+2.0*Q66)*SIN(A) **2*COS(A) **2+ 

*Q22*COS(A) **4 

QB16= (Q11-Q12-2 .0*Q66)*SIN(A) *COS (A) **3+ 

* (Q12-Q22+2.0*Q66)*SIN(A) **3*COS(A) 

QB26= (Qll-Q12-2.0*Q66)*SIN(A) **3*COS( A)+ 

* (Q12-Q22+2.0*Q66 )*SIN(A) *COS(A) **3 
QB66=(Q11+Q22-2.0*(Q12+Q66 ) )*SIN( A) **2*COS (A) **2+ 
*Q66* (SIN (A)** 4+COS (A) **4) 

D11=D11+QB11*DUM 
D22=D22+QB22*DUM 
D66=D66+QB66*DUM 
D12=D12+QB12*DUM 
D16=D16+QB16 *DUM 
D26=D26+QB26 *DUM 
10 CONTINUE 

DK=2.0*(D12+2.0*D66 ) 

RETURN 

END 


SUBROUTINE WSPAR 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NUl , NU2 
COMMON/PCHR/D1 , D2 , H , RHO 

COMMON/SPARA/D11 ,D22 ,D66 ,D12 ,D16 ,D26 ,DK, 
*E1,E2,G,NU1,NU2 
COMMON/TYPE/IPTYPE , NPLYS , IPO 
WRITE (6, 10 ) 

10 FORMAT ( IX, 9HSTIFFNESS/) 
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WRITE (6,20) IPTYPE , NPLYS , IPO 
20 FORMAT (IX, 5HTYPE=, IX, I 1/1 X, 6HNPLYS= , IX, I2/1X , 
*12HORIENTATION=, 814) 

WRITE (6, 30 )E1,E2,G 

30 FORMAT (IX, 3HE1=, E10 . 4/1X, 3HE2=, E10 . 4/1X, 2HG=,E10. 4 ) 
WRITE (6,40) NUl , NU2 

40 FORMAT ( IX, 4HNUl= , F6 . 3/lX , 4HNU2 = , F6 . 3 ) 

WRITE ( 6 , 50 ) Dll , D22 , D66 , D12 , D16 , D26 , DK 
50 FORMAT (/IX, 4HDll=,E10 . 4/1X, 4HD22=,E10 . 4/1X, 4HD66=, 
*E10 . 4/lX, 4HD12=, E10 . 4/1X, 4HD16=, E10 . 4/1X, 4HD26=, 
*E10 . 4/1X, 4HDK= ,E10.4//) 

RETURN 

END 


SUBROUTINE RFLIM 
COMMON/LIMI T/FMAX 
READ (3, 10 ) FMAX 
10 FORMAT (//6X,F7.0) 
RETURN 
END 


SUBROUTINE WFLIM 
COMMON/LIMI T/FMAX 
WRITE (6, 10 ) 

10 FORMAT ( IX, 11HFREQ. LIMIT/) 

WRITE (6, 20 ) FMAX 
20 FORMAT ( IX, 5HFMAX=, F7 . 0// ) 

RETURN 
END 


SUBROUTINE CFREQ 

DIMENSION OMEGA ( 1000 ) , NMAX ( 50 ) 

REAL NUl , NU2 
COMMON/LIMIT/FMAX 
COMMON/NF/MODES ,MMAX , NMAX , OMEGA 
COMMON/PCHR/D1 , D2 , H , RHO 

COMMO N/ S PARA/ Dll , D2 2 , D6 6 , Dl 2 , Dl 6 , D2 6 , DK , 
*E1 , E2 , G , NUl , NU2 
PARAMETER (PI=3 . 14159265359 ) 

MMAX=1000 

MODES=0 

DO 30 M=1 , 50 

NMAX (M)=0 

DO 10 N=1 , 50 

I=MODES+N 

RHOH=RHO*H 

CONST=PI ** 2/SQRT ( RHOH ) 
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DUM=Dll* (M/Dl ) **4+DK* (M/Dl ) **2* (N/D2 ) **2 
*+D22* (N/D2 ) **4-2 . 0*D16* (M/Dl ) **3* (N/D2 ) - 
*2 . 0*D26* (M/Dl ) * (N/D2 ) **3 
RNFREQ=CONST*SQRT (DUM) 

TEST=RNFREQ/ ( 2 . 0 *PI ) 

IF ( TEST . GT . FMAX ) THEN 

GOTO 20 

ELSE 

NMAX ( M )=N 
OMEGA ( I ) =RNFREQ 
END IF 

IF(I.EQ.IOOO) THEN 

MODES=1000 

MMAX=M 

GOTO 40 

END IF 

10 CONTINUE 
20 IF(N.EQ.l) THEN 
MMAX=M-1 
GOTO 40 
END IF 

MODES=MODES+NMAX ( M ) 

30 CONTINUE 
40 RETURN 
END 


SUBROUTINE SORT 

DIMENSION OMEGA (1000) , NMAX (50) ,NMIN(50) 
DIMENSION MNSORT (1000,2), FREQ ( 1000 ) 
COMMON/NF/MODES , MMAX , NMAX , OMEGA 
COMMON/SORT/MNSORT , NOM , FREQ 
PARAMETER (PI=3. 14159265359 ) 

NOM=MODES 
DO 10 M=1 , MMAX 
NMIN(M)=1 
10 CONTINUE 

DO 30 1=1, MODES 
J=0 

FMIN=1 . 0E+10 
DO 20 M=1 , MMAX 
K=J+NMIN( M) 

IF ( NMI N ( M ) . GT . NMAX ( M ) ) THEN 

TEST=1 . 0E+10 

ELSE 

TEST=OMEGA(K) 

END IF 

IF (TEST .LT .FMIN) THEN 
FMI N=OMEGA ( K ) 

MNSORT ( I, 1)=M 
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MNSORT (1,2 )=NMIN (M) 

FREQ ( I )=OMEGA (K ) / ( 2 . 0*PI ) 

END IF 

J=J+NMAX(M) 

20 CONTINUE 

M=MNSORT(I,l) 

N=MNSORT (1,2) 

NMIN(M)=N+1 
30 CONTINUE 
RETURN 
END 


SUBROUTINE WFREQ 

DIMENSION MNSORT (1000,2), FREQ ( 1000 ) 
COMMON/SORT/MNSORT , NOM , FREQ 
WRITE (7, 10) 

10 FORMAT ( IX, 1 9 HN AT URAL FREQUENCIES/) 

WRITE (7,20) NOM 

20 FORMAT (1X,16HNUMBER OF MODES=,l4/) 

WRITE (7, 30) 

30 FORMAT ( 4X, 1HM, 13X, 1HN, 10X, 13HFREQUENCY , HZ/) 
DO 50 1=1, NOM 
M=MNSORT (1,1) 

N=MNSORT (1,2) 

WRITE (7,40 )M, N,FREQ ( I ) 

40 FORMAT ( IX, I4,10X,I4,12X,F9.4) 

50 CONTINUE 
RETURN 
END 
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PROGRAM SOUND ( INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT) 
DIMENSION MNSORT (300,2) , FREQ (300 ) 

DIMENSION AMN( 300 ) ,PSIMN(300 ) 

DIMENSION SMN(300),PHIMN(300) 

DIMENSION CLIM ( 5 ) , EPS ( 3 ) 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL MU,NU,NU1,NU2 
EXTERNAL FTAI , DUMMY 
COMMON/AIC/AMN , PSIMN 
COMMON/DAMP/Z ETAl 1 , Z ETA 9 9 
COMMON/DIST/R 
COMMON/FFREQ/W 

COMMON/FFUNC/Al , A 2 , F0 , MU , NU , P0 
COMMON/LIMI T/FLIM 
COMMON/MED I UM/RHO 0 , CO 
COMMON/NI/NPOINTS 
COMMON/PCHR/D1 , D2 , H , RHO 
COMMON/ SIC/ SMN , PHIMN 
COMMON/SORT/MNSORT , NOM , FREQ 
COMMON/SPARA/Dll ,D22 ,D66 ,D12 ,D16 ,D26 ,DK , 

*El,E2,G,NUl, NU2 
COMMON/TYPE/IPTYPE , NPLYS , IPO 
PARAMETER (PI=3. 14159265359) 

OPEN ( UNIT=1 , FILE= * DATPCR 1 , STATUS= ' OLD ' ) 

OPEN ( UNIT=2 ,FILE= ' DATPSP 1 , STATUS= ' OLD ' ) 

OPEN ( UNIT=3 ,FILE= 1 DATNFL ' , STATUS= 1 OLD ' ) 

OPEN ( UNIT=4 , FILE= ' DATDMP ' , STATUS= ' OLD ' ) 

OPEN ( UNIT=7 , FILE= ' DATMNF ' , STATUS= ' OLD ' ) 

OPEN ( UNIT=8 , FILE= ' DATFF ' , STATUS= ' OLD ' ) 

OPEN ( UNIT=9 , FILE= ' DATDIR ' , STATUS= ' NEW ' ) 

OPEN ( UNIT=10 , FILE= ' DATA ' , STATUS= ' NEW ' ) 

R=20 .0 
RHO0=l . 21 
C0 = 343 .0 
LINES=500 
FMAX=1000.0 
DF=FMAX/LI NES 
CLIM ( 1 )=0 . 0 
CLIM ( 2 )=0 . 0 
CLIM ( 3 )=PI/2 . 0 
CLIM ( 4 ) =0 . 0 
EPS ( 1 ) = 0 . 1 
EPS ( 2 ) = 0 . 0 
EPS ( 3 ) = 0 . 0 
ITEXT=0 
CALL RPCHR 
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CALL RSPAR 
CALL RFLIM 
CALL RDAMP 
CALL RFREQ 
CALL RFFUNC 
CALL WPCHR 
CALL WSPAR 
CALL WFLIM 
CALL WDAMP 
CALL WFFUNC 
TESTl=Al/Dl 
TEST2=A2/D2 

IF ( TEST1 . EQ. .5 .AND .TEST 2 . EQ. . 5 ) THEN 
SYM=4 . 0 

CLIM ( 5 )=PI/2 . 0 

ELSE 

SYM=1 . 0 

CLIM( 5 )=PI*2 . 0 
END IF 
WRITE ( 6 / 10 ) 

10 FORMAT (IX, 10 HRESOLUTION) 

WRITE (6, 20) LINES, FMAX,DF 
20 FORMAT (/IX, 6HLINES= , I4/1X , 5HFMAX=, 

♦F6.0/1X, 3HDF=,F10 .4) 

WRITE (10, 30) 

3 0 FORMAT ( //I X , 6H " FREQ " , 4X , 5H " PWL " , 3X , 

* 10 H" SIGMA, DB" , 4X , 6H" <V2> " , 

*6X, 5H"PWR" , 6X , 8H"EPS ( 3 ) " , 2X , 6H"IERR" , IX, 6H"PNTS"/ ) 

DO 50 1=1, LINES 

FF=I *DF 

W=2.0*PI*FF 

V2=0 . 0 

PWR=0 . 0 

IERR=0 

PWL=0 . 0 

SIGLGT=0 . 0 

NPOINTS=0 

CALL SBIC 

CALL ABIC 

CALL STAMSS V ( NOM , MU , NU , V2 ) 

CALL CAREDB ( CLIM , DUMMY , FTAI , EPS , ITEXT , PWR , IERR ) 

PWR=S YM* PWR 

SIGMA=PWR/(RHOO*CO*V2*Dl*D2 ) 

PWL=10 . 0*ALOG10 ( PWR ) +120 . 0 
SIGLGT=10 .0*ALOG10 (SIGMA) 

WRITE ( 10 , 40 ) FF , PWL , SIGLGT , V2 , PWR , EPS ( 3 ) , IERR , NPOINTS 
40 FORMAT (1X,F8.3,3X,F5.1,3X,F6.1,3X,E10.4,3X,E10.4,3X, 
*E10 . 4 , 4X, II , 3X,I5 ) 

50 CONTINUE 
CLOSE ( 1 ) 
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CLOSE ( 2 ) 
CLOSE ( 3 ) 
CLOSE ( 4 ) 
CLOSE ( 7 ) 
CLOSE ( 8 ) 
CLOSE ( 9 ) 
CLOSE (10 ) 
STOP 
END 


SUBROUTINE RPCHR 
COMMON/PCHR/Dl , D2 , H , RHO 
READ (1,10 )D1,D2 

10 FORMAT (//4X,E10.4/4X,E10. 4) 
READ ( 1 , 20 )H,RHO 
20 FORMAT (3X,E10.4/5X,E10.4) 
RETURN 
END 


SUBROUTINE RSPAR 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NUl , NU2 
COMMON/PCHR/Dl , D2 , H , RHO 

COMMON/ SPARA/Dl 1 , D2 2 , D6 6 , Dl 2 , Dl 6 , D2 6 , DK , 
*El , E2 , G , NUl , NU2 
COMMON/TYPE/ IPTYPE , NPLYS , IPO 
READ( 2, 10 ) IPTYPE, NPLYS, IPO 
10 FORMAT (//7X, 11/8 X, I 2/13X, 814 ) 

READ (2, 20 )E1,E2,G 

20 FORMAT ( 4X, E10 . 4/4X , E10 . 4/4X , E10 . 4) 

READ ( 2 , 30 ) NUl , NU2 
30 FORMAT (5X,F7.3/5X,F7. 3) 

D11=0 . 0 
D22=0 . 0 
D66=0 . 0 
D12=0 . 0 
D16=0 . 0 
D26=0 . 0 

IF ( IPTYPE. EQ.0) THEN 
CONST=H** 3/12 . 0 
D11=CONST*E1/(1.0-NU1*NU2 ) 

D22=CONST*E2/ ( 1 . 0-NUl *NU2 ) 

D66=CONST*G 

D12=CONST*E1*NU2/(1.0-NU1*NU2) 

DK=2.0* (D12+2.0*D66 ) 

END IF 

IF ( IPTYPE. EQ.l) THEN 
CALL CBC 
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END IF 
RETURN 
END 


SUBROUTINE CBC 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NUl , NU2 
COMMON/PCHR/D1 , D2 , H , RHO 

COMMON/SPARA/D11 , D2 2 , D6 6 , Dl 2 , D16 , D26 , DK , 
*E1,E2,G,NU1,NU2 
COMMON/TYPE/ IPTYPE , NPLYS , IPO 
PARAMETER ( PI =3 .14159265359) 

Q11=E1/ ( 1 . 0-NUl*NU2 ) 

Q22=E2/(1.0-NU1*NU2) 

Q12=NU1*E2/ ( 1 . 0-NUl*NU2 ) 

Q66=G 
N2=NPLYS/2 
DZ=H/FLOAT (NPLYS ) 

DO 10 K=1,N2 

ZK=H/2 . O-FLOAT (K ) *DZ 

ZKM1=ZK+DZ 

A=FLOAT ( IPO (K ) )*PI/180.0 
DUM=2 . 0/3 . 0* ( ZKM1**3-ZK**3 ) 

QB11=Q11 *COS (A) **4+ 

*2.0* (Q12+2 . 0*Q66 ) *SIN( A) **2*COS (A) **2+ 

*Q22*SIN( A) **4 

QB12= (Q11+Q22-4 .0*Q66)*SIN(A) **2*COS ( A) **2+ 

*Q12 * ( SIN ( A) **4+COS ( A) **4 ) 

QB22=Q11*SIN (A) **4+ 

*2.0*(Q12+2.0*Q66)*SIN(A) **2*COS(A) **2+ 

*Q22*COS (A) ** 4 

QB16= (Q11-Q12-2 ,0*Q66)*SIN(A) *COS ( A) **3+ 

* (Q12-Q22+2.0*Q66 )*SIN(A) **3*COS ( A) 
QB26=(Q11-Q12-2.0*Q66 )*SIN( A) **3*COS( A)+ 
*(Q12-Q22+2.0*Q66)*SIN(A) *COS ( A) **3 
QB66= (Q11+Q22-2 . 0* (Q12+Q66 ) )*SIN(A) **2*COS ( A) **2+ 
*Q66*(SIN(A) ** 4+COS ( A) ** 4 ) 

D11=D11+QB11*DUM 
D22=D22+QB22*DUM 
D66=D66+QB66*DUM 
D12=D12+QB12*DUM 
D16=D16+QB16 *DUM 
D26=D26+QB26*DUM 
10 CONTINUE 

DK=2 . 0* (D12+2 . 0*D66 ) 

RETURN 

END 
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SUBROUTINE RFLIM 
COMMON/LIMIT/FLIM 
READ (3,10) FLIM 
10 FORMAT (//6X,F 7.0) 

RETURN 
END 


SUBROUTINE RDAMP 
COMMON/DAMP/Z ETAl 1 , Z ETA 9 9 
READ (4,10) ZETAll 
10 FORMAT (//8X,E10. 4) 

READ (4,20) ZETA99 
20 FORMAT (8X,E10. 4) 

RETURN 

END 


SUBROUTINE RFREQ 

DIMENSION MNSORT (300,2) ,FREQ(300) 
COMMON/SORT/MNSORT , NOM , FREQ 
READ( 7,10) NOM 
10 FORMAT ( //I 7 X, 14///) 

DO 30 1=1, NOM 

READ (7,20) MNSORT (1,1), MNSORT (1,2), FREQ ( I ) 
20 FORMAT (1X,I4,10X,I4,12X,F9.4) 

30 CONTINUE ' 

RETURN 

END 


SUBROUTINE RFFUNC 
REAL MU, NU 

COMMON/FFUNC/A1 , A2 , F0 ,MU , NU , P0 
READ (8, 10) F0 , P0 

10 FORMAT (//4X,E10.4/4X,E10. 4) 
READ (8, 20) MU , NU 

20 FORMAT (////4X,E10.4/4X,E10.4) 
READ (8,30) A1,A2 

30 FORMAT (// //4X, Ell .5/4X,E11.5) 

RETURN 
END 


SUBROUTINE WPCHR 
COMMON/PCHR/D1 , D2 , H , RHO 
WRITE (6, 10 ) 

10 FORMAT ( IX, 14HPANEL GEOMETRY/) 

WRITE (6, 20 )D1,D2 

20 FORMAT (IX, 3HD1=,E10 . 4/1X, 3HD2=,E10 . 4 ) 
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WRITE (6,30) H,RH0 

30 FORMAT ( IX , 2HH= , E10 . 4/1X , 4HRHO= , E10 . 4// ) 

RETURN 
END 


SUBROUTINE WSPAR 

INTEGER IPTYPE , NPLYS , IPO ( 8 ) 

REAL NUl , NU2 
COMMON/PCHR/D1 , D2 ,H, RHO 

COMMON/SPARA/D11 , D22 , D66 , D12 , D16 , D26 , DK , 
*E1,E2,G,NU1,NU2 
COMMON/TYPE/IPTYPE , NPLYS , IPO 
WRITE (6, 10 ) 

10 FORMAT ( IX, 9HSTIFFNESS/) 

WRITE (6,20) IPTYPE , NPLYS , IPO 
20 FORMAT ( IX , 5HTYPE= ,1X,I1/1X, 6HNPLYS= , 

*1X,I2/1X, 12HORIENTATION=, 814 ) 

WRITE (6, 30 )E1,E2,G 

30 FORMAT ( IX, 3HEl= , E10 . 4/lX , 3HE2= , E10 . 4/1X , 2HG= , E10 . 4 ) 

WRITE (6, 40) NUl, NU2 

40 FORMAT (IX, 4HNUl=,F6 . 3/lX, 4HNU2=,F6 . 3 ) 

WRITE ( 6 , 50 ) Dll , D22 , D66 , D12 , D16 , D26 , DK 
50 FORMAT ( /IX, 4HDll=, E10 . 4/lX, 4HD22=, 

*E10 . 4/lX, 4HD66=,E10 . 4/lX, 4HD12=,E10 . 4/lX, 4HD16=, 
*E10 . 4/lX, 4HD26=,E10 . 4/lX , 4HDK= ,E10 . 4// ) 

RETURN 

END 


SUBROUTINE WFLIM 
COMMON/LIMIT/FLIM 
WRITE (6, 10) 

10 FORMAT (IX, 11 HFREQ. LIMIT/) 

WRITE (6,20) FLIM 
20 FORMAT (IX, 5HFLIM=, F7 . 0// ) 
RETURN 
END 


SUBROUTINE WDAMP 
COMMON/DAMP/Z ETAll , Z ETA 9 9 
WRITE (6, 10 ) 

10 FORMAT ( IX, 7HDAMPING/) 

WRITE (6, 20 ) ZETA11 
20 FORMAT ( IX , 7HZETAll=, E10 . 4 ) 
WRITE (6,30) Z ETA 9 9 
30 FORMAT (IX, 7HZETA99= , E10 . 4// ) 
RETURN 
END 
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SUBROUTINE WFFUNC 
REAL MU, NU 

COMMON/FF UNC/ Al , A 2 , F 0 , MU , NU , P 0 
WRITE (6/10 ) 

10 FORMAT ( IX, 9HMAGNITUDE/) 

WRITE (6, 20) F0,P0 

20 FORMAT ( IX, 3HF0=, E10 . 4/1X, 3HP0=,E10 . 4 ) 

WRITE (6, 30 ) 

30 FORMAT (//IX, 5HPHASE/) 

WRITE (6,40) MU , NU 

40 FORMAT ( IX, 3HMU= , E10 . 4/1X, 3HNU= , E10 . 4 ) 

WRITE (6,50) 

50 FORMAT ( //IX, 15HSHAKER LOCATION/) 

WRITE (6, 60) A1,A2 

6 0 FORMAT ( IX , 3HAl= , Ell . 5/1X , 3HA2= , Ell . 5// ) 

RETURN 
END 


SUBROUTINE SBIC 

DIMENSION MNSORT (300,2) ,FREQ( 300 ) 

DIMENSION PHIMN (300), SMN (300) 

REAL MU , NU , NUM 
COMMON/DAMP/ZETA11 , Z ETA 99 
COMMON/FFREQ/W 

COMMON/FF UNC/ Al ,A2,F0,MU,NU,P0 
COMMON/PCHR/D1 , D2 , H , RHO 
COMMON/SIC/SMN, PHIMN 
COMMON/SORT/MNSORT , NOM , FREQ 
PARAMETER (PI=3. 1415926 535 9) 

DO 10 1=1, NOM 
M=MNSORT (1,1) 

N=MNSORT (1,2) 

WMN=2 . 0 *PI *FREQ ( I ) 

ZMN=ZETAll*FREQ ( 1 ) /FREQ ( I ) +ZETA99 

DEN=SQRT ( ( 1 . 0-W** 2/WMN* * 2 ) ** 2+ ( 2 . 0 *ZMN*W/WMN ) ** 2 ) 
NUM=SIN(M*PI*A1/D1)*SIN(N*PI*A2/D2 ) 

DUM=NUM/DEN 

SMN( I )=4 . 0*FO/(D1*D2*RHO*H*WMN**2 ) *DUM 
NUM=2 . 0*ZMN*W/WMN 
DEN=1 . 0-W** 2/WMN* *2 

PHIMN( I)=ATAN2 (NUM,DEN)+MU*PI/180 .0 
10 CONTINUE 
RETURN 
END 


SUBROUTINE ABIC 
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DIMENSION MNSORT (300,2), FREQ ( 300 ) 

DIMENSION AMN( 300 ) , PSIMN ( 300 ) 

REAL MU , NU , NUM 
COMMON/AIC/AMN, PSIMN 
COMMON/DAMP/ZETAll , ZETA99 
COMMON/FFREQ/W 

COMMON/FFUNC/A1 , A2 , F0 ,MU , NU , P0 
COMMON/PCHR/D1 ,D2 ,H, RHO 
COMMON/SORT/MNSORT , NOM , FREQ 
PARAMETER (PI=3. 1415926535 9) 

DO 10 1=1, NOM 
M=MNSORT (1,1) 

N=MNSORT (1,2) 

WMN=2 . 0 *PI *FREQ ( I ) 

ZMN=ZETAll*FREQ ( 1 ) /FREQ ( I ) +ZETA99 

DEN=SQRT ( ( 1 . 0-W**2/WMN**2 ) **2+ ( 2 . 0*ZMN*W/WMN) **2 ) 
NUM= ( 1 . 0-COS ( M*PI ) ) * (1 . 0-COS (N*PI ) ) 

DUM=NUM/DEN 

AMN( I)=4.0*P0/(M*N*PI**2*RHO*H*WMN**2)*DUM 
NUM=2 . 0*ZMN*W/WMN 
DEN=1 . 0-W** 2/WMN** 2 

PSIMN(I)=ATAN2 ( NUM , DEN ) +NU*PI/180 .0 
10 CONTINUE 
RETURN 
END 


SUBROUTINE STAMSSV ( NOM, MU, NU , V2 ) 

DIMENSION SMN( 300 ) , AMN( 300 ) 

REAL MU , NU 

COMMON/AIC/AMN 

COMMON/FFREQ/W 

COMMON/SIC/ SMN 

PARAMETER ( PI=3 . 14159265359 ) 

PCOEFF=COS ( ( MU-NU ) *PI/180 . 0 ) 

DO 10 1=1, NOM 

V2INC=SMN( I ) **2+AMN( I ) **2+2 . 0*SMN( I ) *AMN( I ) * 
*PCOEFF 
V2=V2+V2INC 
10 CONTINUE 

V2=V2*W**2/8.0 

RETURN 

END 


FUNCTION FTAI (THETA, PHI ) 

DIMENSION AMN( 300 ) , PSIMN ( 300 ) 
DIMENSION SMN( 300 ) , PHIMN( 300 ) 
DIMENSION MNSORT (300, 2 ) ,FREQ( 300) 
REAL K 
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COMPLEX PCMPLX, PSUM,P2 , JM, JN 
COMMON/AIC/AMN, PSIMN 
COMMON/DIST/R 
COMMON/FFREQ/W 
COMMON/MED I UM/RHO 0 , CO 
COMMON/NI/NPOINTS 
COMMON/PCHR/D1 , D2 
COMMON/SIC/ SMN , PHIMN 
COMMON/SORT/MNSORT , NOM , FREQ 
PARAMETER (PI=3. 1415 9265359 ) 

K=W/C0 

FF=W/(2.0*PI) 

TEST=2 . 0 *FF 

CONST=-RHOO*W**2/( 2 . 0*PI*R) 

PSUM= ( 0 . 0 , 0 . 0 ) 

TAI=0 . 0 

NPOINTS=NPOINTS+l 
DO 10 1=1, NOM 

IF ( ( I . GT . 30 ) . AND . ( FREQ ( I ) . GT . TEST ) ) GOTO 20 
IF(SMN(I) . NE . 0 . 0 . OR . AMN ( I ) . NE . 0 . 0 ) THEN 
M=MNSORT (1,1) 

N=MNSORT (1,2) 

ARG1=K*SIN ( THETA ) *COS ( PHI ) 

ARG2=ARG1 ** 2 

ARG3=M*PI/D1 

ARG4=ARG3**2 

IF ( ARG2 . EQ. ARG4 ) THEN 

VR=0 . 0 

VI=Dl/2 . 0*SIGN (1.0 , ARG1 ) 

ELSE 

DUM=ARG3/ ( ARG4-ARG2 ) 

ARG=ARGl*Dl 

VR=DUM* ( 1 . O-COS ( M*PI ) *COS ( ARG ) ) 

VI=DUM* ( -COS (M*PI ) *SIN( ARG) ) 

END IF 

JM=CMPLX ( VR , VI ) 

ARG1=K*SIN ( THETA ) *SIN ( PHI ) 

ARG2=ARG1**2 

ARG3=N*PI/D2 

ARG4=ARG3**2 

IF ( ARG2 ,EQ. ARG4 ) THEN 

VR=0 . 0 

VI=D2/2*SIGN( 1 . 0 , ARGl ) 

ELSE 

DQM=ARG3/( ARG4-ARG2 ) 

ARG=ARG1*D2 

VR=DUM* ( 1 . O-COS ( N*PI ) *COS ( ARG ) ) 

VI=DUM* ( -COS ( N*PI ) *SIN( ARG) ) 

END IF 

JN=CMPLX(VR, VI) 
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VR=SMN ( I ) *COS ( PHIMN ( I ) ) +AMN ( I ) *COS ( PSIMN ( I ) ) 
VI=-SMN(I)*SIN(PHIMN(I) ) -AMN( I ) *SIN( PSIMN( I ) ) 
PSUM=PSUM+CMPLX ( VR , VI ) * JM* J N 
END IF 

10 CONTINUE 
20 PCMPLX=CONST*PSUM 

P 2=PCMPLX*CON JG ( PCMPLX ) 

TAI=REAL(P2)/(2.0*RHO0*C0 ) 

FTAI=TAI *R**2*SIN( THETA ) 

RETURN 

END 

C 

C 

SUBROUTI NE CAREDB ( CLIM , FLIM , FX , EPS , ITEXT , SUM , IERR ) 
c******************************************************** 
C* * 

C* PURPOSE: * 

C* TO EVALUATE THE DOUBLE INTEGRAL OF F(X,Y)DYDX USING * 
C* CAUTIOUS ROMBERG INTEGRATION. (THE INTEGRAL IS COM- * 
C* PUTED BY CALCULATING THE INTEGRAL OVER SUITABLY * 

C* SMALL SUBINTERVALS OF THE INTERVALS OF INTEGRATION, * 
C* SIMULTANEOUSLY ALONG THE INNER INTEGRAL AT GIVEN * 
C* OUTER LIMIT VALUES AND ALONG THE ENTIRE OUTER INTE- * 
C* GRAL, AND SUMMING THE RESULTS. THE INTEGRAL OVER * 
C* EACH SUBINTERVAL IS COMPUTED THROUGH ROMBERG * 

C* EXTRAPOLATION WHERE THE VALUE OF THE INTEGRAL IS * 
C* ACCEPTED WHEN THE COMPUTED ERROR IS LESS THAN: * 

C* MAXIMUM (ABSOLUTE ERROR, RELATIVE ERROR * CURRENT * 
C* ESTIMATE OF THE INTEGRAL). ) * 

C* * 

C* * 

C* USE: * 

C* CALL CAREDB (CLIM, FLIM, FX, EPS, ITEXT, SUM, IERR) * 

C* * 

c* * 

C* PARAMETERS: * 

C* * 

C* CLIM AN INPUT ONE-DIMENSIONAL ARRAY DIMENSIONED* 

C* FIVE. SPECIFIES THE INTEGRATION LIMITS. * 

C* * 

C* CLIM ( 1 ) INTEGRATION LIMIT INDICATOR. * 

C* =0. F1(X) AND F2 (X) ARE CONSTANT * 

C* FUNCTIONS. * 

C* =1. F1(X) AND F2 (X) ARE NON-CONSTANT * 

C* FUNCTIONS. * 

C* * 

C* CLIM ( 2 ) LOWER LIMIT OF INTEGRATION FOR X, A. * 

C* * 

C* CLIM ( 3 ) UPPER LIMIT OF INTEGRATION FOR X, B. * 

C* * 


244 



C* CLIM ( 4 ) LOWER LIMIT OF INTEGRATION FOR Y, FI(X), * 
C* WHEN CLIM(1 )=0 . * 
C* * 
C* CLIM ( 5 ) UPPER LIMIT OF INTEGRATION FOR Y, F2(X), * 
C* WHEN CLIM(1)=0. * 
C* * 
C* IF CLIM ( 1 )=1 . , THE PARAMETERS CLIM (4) AND CLIM(5) * 
C* ARE NOT USED BY SUBROUTINE CAREDB. * 
C* * 
C* FLIM AN INPUT PARAMETER SPECIFYING THE NAME OF * 
C* A USER-PROVIDED SUBROUTINE WITH ARGUMENTS X, P, AND * 
C* Q USED TO EVALUATE THE INNER LIMITS OF INTEGRATION * 
C* WHEN CLIM ( 1 )=1 . THE ARGUMENTS P AND Q ARE RESPECT- * 
C* IVELY THE LOWER AND UPPER LIMITS F1(X) AND F2(X), * 
C* EVALUATED AT X. THE NAME OF THE SUBROUTINE MUST * 
C* BE DECLARED IN AN EXTERNAL STATEMENT IN THE CALLING * 
C* PROGRAM WHEN CLIM(1)=1. * 
C* * 
C* IF CLIM(1 )=0 . , FLIM IS A DUMMY ARGUMENT. THE USER * 
C* NEED NOT PROVIDE A SUBROUTINE AND THE EXTERNAL * 
C* STATEMENT IS NOT NEEDED. * 
C* * 
C* FX AN INPUT PARAMETER SPECIFYING THE NAME OF * 
C* A USER-PROVIDED FUNCTION SUBPROGRAM WITH ARGUMENTS X* 
C* AND Y USED TO EVALUATE F(X,Y). THE NAME OF THE SUB-* 
C* PROGRAM MUST BE DECLARED IN AN EXTERNAL STATEMENT * 
C* IN THE CALLING PROGRAM. * 
C* * 
C* EPS AN INPUT/OUTPUT ONE-DIMENSIONAL ARRAY * 
C* DIMENSIONED THREE SPECIFYING MAXIMUM ALLOWABLE * 
C* ERROR CRITERIA. * 
C* * 
C* EPS ( 1 ) INPUT VARIABLE CONTAINING THE MAXIMUM * 
C* ALLOWABLE RELATIVE ERROR. * 
C* * 
C* EPS ( 2 ) INPUT VARIABLE CONTAINING THE MAXIMUM * 
C* ALLOWABLE ABSOLUTE ERROR. * 
C* * 
C* EPS ( 3 ) OUTPUT VARIABLE CONTAINING AN ESTIMATE OF * 
C* THE ERROR IN COMPUTING THE INTEGRAL OVER * 
C* THE ENTIRE DOUBLE INTERVAL OF INTEGRATION. * 
C* * 
C* ITEXT AN INPUT INTEGER CODE PROVIDING THE USER * 
C* WITH THE OPTION OF HAVING CAREDB PRINT OUT IRREGULAR* 
C* FUNCTION BEHAVIOR (INCLUDING END POINT SINGULARI- * 
C* TIES, JUMP DISCONTINUITIES AND REGIONS INDICATING * 
C* FUNCTIONAL NOISE), AND INTERVALS OF UNSUCCESSFUL * 
C* INTEGRATIONS WITH CORRESPONDING ESTIMATES. THIS * 
C* PRINTOUT SHOULD ONLY BE REQUESTED FOR PROBLEMS * 
C* WHICH MUST BE RERUN DUE TO UNSATISFACTORY RESULTS * 


245 



THE FIRST TIME 


C* 

C* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 

c* 


* 

ITEXT = 0 NO PRINTOUT REQUESTED. * 

ITEXT = 1 PRINTOUT REQUESTED WITH RESPECT 

TO INNER INTEGRAL ONLY. * 

ITEXT = 2 PRINTOUT REQUESTED WITH RESPECT * 

TO OUTER INTEGRAL ONLY. * 

ITEXT = 3 PRINTOUT ON BOTH INTEGRALS * 

REQUESTED * 

* 

SUM AN OUTPUT VARIABLE CONTAINING THE COMPUTED * 

VALUE OF THE DOUBLE INTEGRAL. (THE VALUES OF SUM * 
AND EPS ( 3 ) SHOULD BE CONSIDERED FOR ALL POSSIBLE * 
RETURNS - SEE IERR. ) * 

* 


IERR AN OUTPUT ERROR PARAMETER * 

=0 INTEGRATION SUCCESSFUL - NORMAL RETURN * 

(RELATIVELY SMOOTH INTEGRAND) * 

=1 INTEGRATION SUCCESSFUL (WITH RESPECT TO AT * 

LEAST ONE VALUE OF X, FUNCTION SINGULARITIES * 

WITHIN THE INNER INTEGRAL WERE SUCCESSFULLY * 
HANDLED DURING THE INTEGRATION PROCESS.) * 

=2 INTEGRATION SUCCESSFUL (IN ONE OR MORE SUB- * 

INTERVALS OF THE INTERVALS OF INTEGRATION, * 

FUNCTION SINGULARITIES IN X, OR BOTH X AND Y, * 
WERE SUCCESSFULLY HANDLED DURING THE INTE- * 

GRATION PROCESS . ) * 

=3 INTEGRATION SUCCESSFUL (WITH RESPECT TO AT * 

LEAST ONE VALUE OF X, NO REGULAR FUNCTION * 

BEHAVIOR WAS DETECTED WITHIN THE INNER INTE- * 
GRAL BUT THE CONVERGENCE CRITERION WAS SATIS- * 
FIED. ) (THE INTEGRAL VALUE SHOULD BE ACCEPTED* 
WITH CAUTION.) * 

=4 INTEGRATION SUCCESSFUL (IN ONE OR MORE SUB- * 

INTERVALS OF BOTH INTERVALS OF INTEGRATION NO * 
REGULAR FUNCTION BEHAVIOR WAS DETECTED IN X, * 
OR BOTH X AND Y, BUT THE CONVERGENCE CRI- * 

TERION WAS SATISFIED.) (THE INTEGRAL VALUE * 
SHOULD BE ACCEPTED WITH CAUTION.) * 

=5 UNSUCCESSFUL INTEGRATION (IN SUBROUTINE * 

QXZ056 , *QXZNA508 * 

THE STORAGE ARRAY, *TS*, CONTAINING THE FUNC- * 
TION VALUES STILL TO BE USED DURING THE INTE- * 
GRATION PROCESS ALONG Y AND AT SOME VALUE OF * 
X IS EXHAUSTED RELATIVE TO SUPPLIED ERROR * 

CRITERION.) * 

=6 UNSUCCESSFUL INTEGRATION (WHILE EVALUATING THE * 

FUNCTION ALONG Y AT A GIVEN VALUE OF X, A * 

SUBINTERVAL SMALLER THAN THE MINIMUM ALLOW- * 

ABLE STEP IS REQUESTED, OR THE MAXIMUM NUMBER * 
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OF SUBINTERVALS STILL TO BE CONSIDERED EX- * 
CEEDS THE ALLOWABLE MAXIMUM SIZE OF THE SUB- * 
INTERVAL STACK.) * 

=7 UNSUCCESSFUL INTEGRATION (THE STORAGE ARRAY, * 
*TS* , IS EXHAUSTED RELATIVE TO SUPPLIED ERROR * 
CRITERION.) * 

=8 UNSUCCESSFUL INTEGRATION (A SUBINTERVAL SMALL- * 
ER THAN THE MINIMUM ALLOWABLE STEP IS REQUES- * 
TED, OR THE MAXIMUM NUMBER OF SUBINTERVALS * 

STILL TO BE CONSIDERED EXCEEDS THE ALLOWABLE * 
MAXIMUM SIZE OF THE SUBINTERVAL STACK.) * 

=9 INTEGRATION SUCCESSFUL (THE INNER INTEGRAL LIMIT* 
SMALLER THAN MINIMUN ALLOWABLE STEP ONE OR * 
MORE TIMES DURING THE INTEGRATION PROCESS. A * 
SINGLE TRAPEZOIDAL SUM IS USED TO ESTIMATE THE* 

INTEGRAL AND INTEGRATION CONTINUES. ) * 

* 

UPON RETURN FROM HIS PARAMETER SHOULD BE * 

TESTED IN THE CALLING PROGRAM. * 


REQUIRED ROUTINES 


QXZ056 


QXZNA509 


LANGUAGE 


FORTRAN 


AUTHOR / IMPLEMENTER 

COMPUTER SCIENCES / ACD PROGRAMMER SUPPORT 
CORPORATION GROUP, EXT. 3548 


DATE RELEASED 


JUNE 1975 


LATEST REVISION 


APRIL 1983 
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The following pages are a listing of program PREDICT: 


PROGRAM PREDICT 

REAL PWL ( 500 ) , SIGMA ( 500 ) , V2 ( 500 ) , PWR ( 500 ) 

REAL SDATA (500,2) , ADATA (500,2), PRED (500,2) 

REAL FF, DEN, DNUM, DUM, DIFF , TEST 
COMMON/INPUT/PWL , SIGMA , V2 , PWR 
CONST=40 .695 
WRITE ( * , 10 ) 

10 FORMAT (/IX, ' PLACE THE INPUT DATA DISK IN SLOT B:'/) 
WRITE (*,20 ) 

20 FORMAT ( /IX , 1 PLACE THE OUTPUT DATA DISK IN SLOT A:'/) 
PAUSE 

WRITE ( * , 30 ) 

30 FORMAT (/IX, 1 ENTER THE NAME OF THE STRUCTUREBORNE DATA 
*FILE ' / ) 

OPEN ( 1 , FILE= ' ' , STATUS= ' OLD ' ) 

WRITE ( * , 40 ) 

40 FORMAT (/IX,' ENTER THE NAME OF THE AIRBORNE DATA 
*FILE ' / ) 

OPEN ( 2 , FILE= ' ' ,STATUS='OLD' ) 

WRITE (*,50) 

50 FORMAT ( /IX, 1 ENTER THE NAME OF THE COMBINED INPUTS DATA 
*FILE ' / ) 

OPEN ( 3 ,FILE= ' ' ,STATUS='OLD' ) 

60 OPEN ( 4, FILE=' PREDICT. PRN' , STATUS= ' NEW ' ) 

WRITE ( * , 70 ) 

70 FORMAT (/IX,' TANDY 2000 READING THE STRUCTUREBORNE DATA 
*FILE ' / ) 

CALL RFILE(l) 

CLOSE ( 1 ) 

DO 80 J=1 , 500 
SDATA ( J , 1 ) =PWL ( J ) 

SDATA (J , 2 )=10 . 0** (SIGMA ( J ) /10 . 0 ) 

80 CONTINUE 

WRITE (*,90) 

90 FORMAT (/IX, ' TANDY 2000 READING THE AIRBORNE DATA 
*FILE 1 / ) 

CALL RFILE ( 2 ) 

CLOSE ( 2 ) 

DO 100 J=1 , 500 
ADATA (J,1)=PWL(J) 

ADATA (J,2)=10.0** (SIGMA (J)/10.0) 

100 CONTINUE 

WRITE (*,110 ) 

110 FORMAT ( /IX , 1 TANDY 2000 READING THE COMBINED INPUTS 
*DATA FILE'/) 

CALL RFILE ( 3 ) 
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CLOSE ( 3 ) 

WRI TE (*,120) 

120 FORMAT ( /IX , ' TANDY 2000 NOW BUSY COMPUTING THE 
♦RESULTS'/) 

DO 130 J=1 , 500 

DNUM=SD ATA ( J , 2 ) * ( PWR ( J ) -ADATA ( J , 2 ) *CONST * V2 ( J ) ) 

DEN=SD ATA ( J , 2 ) -ADATA ( J , 2 ) 

DIFF=10 . 0*ALOG10 ( ADATA ( J , 2 ) /SDATA ( J , 2 ) ) 

TEST=ABS(DIFF) 

IF ( TEST . LT .2.0) THEN 

PRED ( J , 1 ) = 0 . 0 

ELSE 

DUM=DNUM/DEN 
IF(DUM.LE.O.O) THEN 
PRED ( J , 1 ) = 0 . 0 
ELSE 

PRED( J,l)=10 .0*ALOG10 (DUM) +120.0 
END IF 
END IF 

130 CONTINUE 

DO 140 J=1 , 500 

DNUM=ADATA ( J , 2 ) * ( PWR ( J ) -SDATA ( J , 2 ) *CONST*V2 ( J ) ) 
DEN=ADATA ( J, 2 ) -SDATA ( J , 2 ) 

DIFF=10 .0*ALOG10 ( ADATA ( J , 2 ) /SDATA ( J , 2 ) ) 

TEST=ABS ( DIFF ) 

IF ( TEST. LT. 2.0) THEN 

PRED ( J , 2 )=0 . 0 

ELSE 

DUM=DNUM/DEN 
IF (DUM. LE. 0.0) THEN 
PRED ( J, 2 )=0 . 0 
ELSE 

PRED( J, 2) =10 .0*ALOG10 (DUM) +120 .0 
END IF 
END IF 

140 CONTINUE 

WRITE (*,150 ) 

150 FORMAT (/IX, ' TANDY 2000 NOW WRITING THE RESULTS TO 
♦DISK' /) 

DO 170 J=1 , 500 
FF=2 . 0*FLOAT ( J ) 

WRITE (4, 160 ) FF , PWL ( J ) , SDATA ( J , 1 ) , PRED ( J , 1 ) , ADATA ( J , 1 ) , 
♦PRED (J, 2) 

160 FORMAT ( IX, F8. 3 , 1X,F6 . 1 , 1X,F6 . 1 , IX ,F6 . 1 , IX ,F6 . 1 , IX, F6 . 1 ) 
170 CONTINUE 
CLOSE ( 4 ) 

WRITE ( ♦ , 180 ) 

180 FORMAT ( /IX , ' PROGRAM COMPLETE'/) 

STOP 

END 
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o o 


SUBROUTINE RFILE ( IUNIT ) 

REAL PWL ( 500 ) , SIGMA ( 500 ) , V2 ( 500 ) , PWR ( 500 ) 
INTEGER IUNIT 

COMMON/ I NPUT/PWL , SIGMA , V2 , PWR 
READ ( IUNIT, 10 ) 

10 FORMAT (////) 

DO 30 J=1 , 500 

READ ( IUNIT , 20 ) PWL ( J ) , SIGMA ( J ) , V2 ( J ) , PWR ( J ) 
20 FORMAT (13X,F5.1,3X,F6.1,3X f E10.4,3X,E10.4) 
30 CONTINUE 
RETURN 
END 
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Appendix X 

VALIDATION OP THE COMPUTER PROGRAMS 
A_j_ GRAPHICAL APPROACH 

It was decided that several simple tests of programs 
MODAL and SOUND should be performed to insure that the these 
programs were operating correctly (free of bugs). The tests 
were designed so that the results produced by the programs 
were expected to exhibit simple physical principles known in 
advance. The tests were also designed so that the results 
could be presented in simple graphical form. This graphical 
approach provided a reliable method for debugging the pro- 
grams almost at a glance. 

VALIDATION OF PROGRAM MODAL 

It was decided that an acceptable test of program MODAL 
(the program which performs the eigenvalue analysis) would 
be to plot the natural frequencies computed by the program 
versus mode numbers m and n. This should be a good test of 
the program since an important characteristic of this type 
of plot (for plates and shells) is known in advance. In 
reference 56, it is shown that one of the characteristics of 
the flexural modes of all plates and shells is that as mode 
numbers m and n increase without bound, the various curves 
of natural frequency versus mode number m, obtained for 
different values of mode number n, tend to collapse on top 
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of one another. (For example, see reference 56, figure 
5.5.2. ) 

A small computer program was written so that the natur- 
al frequencies computed by program MODAL and stored in file 
DATMNF could be called up and plotted versus mode numbers m 
and n. An example of the output produced by this program is 
given in figure (AX-1). The data shown in figure (AX-1) are 
for the AA 2024 aluminum plate (plate no. 1). From the data 
in this figure, it is seen that the natural frequencies com- 
puted by program MODAL exhibit the expected physical charac- 
teristic. Similar results were obtained for plates 2 
through 6. 

VALIDATION OF PROGRAM SOUND 

Since program SOUND was the workhorse of the analytical 
study, the tests of this program were more extensive. The 
tests hinged on checks of two simple quantities, however. 
First, directivity of the sound radiation was checked by 
plotting acoustic intensity versus azimuthal angle (for 
different values of the polar angle), and secondly, the 
input parameters to program SOUND were adjusted so that the 
coincidence phenomenon could be observed and verified in 
terms of the radiation efficiency of the structure. 
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Checking Sound Field Directivity 


Plots of sound field directivity provided a graphical 
means for checking at least three different physical princi- 
ples. The directivity plots also give the programmer an 
idea of what the numerical integration routine is using to 
compute the overall sound power radiated. (It should be 
noted that the numerical integration routine, CAREDB, is by 
far the largest and most complicated portion of program 
SOUND. ) 

First of all, directivity plots provide a means of 
checking symmetry. For example, the sound field produced by 
two of the three forcing functions used in this study, viz. 
the normally incident acoustic input and the point vibra- 
tional input located at the center of the plate, should be 
symmetric about the C <t> = 0) and the B 2 (<l> = 90) 

axes. Conversely, the sound field produced by the third 
forcing function, viz. the point vibrational input located 
near the corner of the plate, should have no axis of sym- 
metry. 

Secondly, directivity plots give an indication of the 
degree of evanescence in the sound field. Since structure- 
borne noise tends to be an inefficient form of noise radia- 
tion, the sound field produced by a point vibrational input 
to a plate should be very reactive and highly directional in 
character. In contrast, a normally incident, spatially 
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uniform acoustic input produces a vibrational response that 
generates noise much more efficiently. Therefore, the sound 
field associated with the airborne noise is expected to be 
much less directional in character by comparison (more like 
the sound field of a plane wave or a simple monopole). 

Thirdly, directivity plots permit examination of the 
frequency dependent characteristics of the sound field. For 
example, as the forcing frequency becomes larger, more and 
more higher order vibrational modes participate in the noise 
generation. As long as the forcing frequency is below coin- 
cidence, the higher order modes will be inefficient, highly 
directional noise radiators. This means that as the forcing 
frequency is increased, the sound field (and therefore the 
directivity plots) should become more directional. 

As part of the validation procedure, program SOUND was 
modified so that these properties of frequency dependence, 
evanescence, and symmetry (or asymmetry) could be observed 
and verified. The program was altered so that each time the 
numerical integration routine calls function FTAI (the rou- 
tine which calculates the intensity at a point in space), 
the azimuthal angle <j>, the polar angle 0, and the out- 
ward normal component of the acoustic intensity I are re- 
corded in the data file DATDIR. (The portions of the pro- 
gram which take advantage of symmetry were removed. Also 
note that the original listing of program SOUND in Appendix 
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IX defines and opens file DATDIR but does not use it.) 

Since the storage of <f> , 0, and I for each call of func- 
tion FTAI for a 500 line spectra would produce an excessive- 
ly large data base, program SOUND was further modified so 
that data was computed for only one forcing frequency at a 
time (instead of 500). 

For the purposes of validating program SOUND, results 
are presented here for the cases of the forcing frequency 
equal to 500 Hz and the forcing frequency equal to 1000 Hz. 
From the discussion of the calculation scheme given in 
Appendix IX, it can be seen that if the forcing frequency is 
500 Hz, all of the modes with natural frequencies between 0 
and 1000 Hz participate in the analysis. Similarly, if the 
forcing frequency is 1000 Hz, all of the modes with natural 
frequencies between 0 and 2000 Hz participate in the anal- 
ysis. 

In addition to the restrictions on the forcing freq- 
uency, the results given here are limited to the cases of 
polar angle 9 equal to 22.5, 45, 67.5 and 90 degrees. 
Although program SOUND generates data for many, many other 
values of the angle 9 (in fact, it must do so to perform 
the numerical integration), it was decided that the data 
from the four angles stated were fairly representative of 
the entire sound field and certainly adequate for validation 
of the program. 
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The results of test runs of the modified program SOUND 
for the case of the AA 2024 aluminum plate (plate no. 1) 
subjected to a 1 Pa peak, normally incident, spatially uni- 
form acoustic load at forcing frequencies of 500 and 1000 Hz 
are presented in figures (AX-2) and (AX-3), respectively. 
Note that for the case of the forcing frequency equal to 500 
Hz that the sound field is nearly uniform with only small 
variations in I (versus c}>) even for large values of © 

(e.g. 0 = 90). In contrast, the plots of I versus i- n 
figure (AX-3) (forcing frequency equal to 1000 Hz) show more 
directional preference for the noise radiation with the 
larger and more rapid variations occuring at the larger 
values of 0 (e.g. 0 = 67.5, 90). Thus, the two figures 
verify that the sound field becomes more directional as the 
forcing frequency is increased. Also note that the direc- 
tivity plots in the two figures are symmetric about the 8^ 

(<j> = 0) and the 8 ^ ( = 90) axes as expected. 

The results of test runs of program SOUND for the case 
of the AA 2024 aluminum plate subjected to a . 01 N peak, 
point vibrational load located at the center of the plate at 
forcing frequencies of 500 and 1000 Hz are presented in 
figures (AX-4) and (AX-5), respectively. Note that in 
figure (AX-4), the sound field possesses much larger varia- 
tions in the intensity level I versus <b for all values of 
0 when compared to the corresponding cases of the acoustic 
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input (figure (AX-2)). Figures (AX-4) and (AX-5) also 
corroborate the earlier findings for the acoustic input by 
once again indicating that the sound field becomes more 
directional as the forcing frequency is increased. Again 
note that the directivity plots in figures (AX-4) and (AX-5) 
are symmetric about the 6 ^ (<j> = 0 ) and the 82 (<f> = 

90) axes as expected. 

Finally, the results of test runs of program SOUND for 
the case of the AA 2024 aluminum plate subjected to . 01 N 
peak, point vibrational load located near the corner of the 
plate (ot^ = . 06033 m and a 2 = . 13547 m) at forcing freq- 
uencies of 500 and 1000 Hz are presented in figures (AX- 6 ) 
and (AX-7), respectively. The first and most obvious char- 
acteristic of the plots in these two figures is that the 
intensity level I is no longer symmetric (as expected) 
about any axis in <{> or 0. Also, recall that the plots 
in figures (AX-2) through (AX-5) were, for the most part, 
limited to roughly a 15 dB variation in the intensity level 
over the hemispherical surface. (Only one case in figure 
(AX-4) shows substantially more than a 15 dB variation.) 
Figures (AX- 6 ) and (AX-7), however, tend to exhibit larger 
variations than 15 dB in the intensity level. Figure (AX-7) 
also exhibits larger differences in the shape of the curves 
(as the angle 0 is varied) than was observed in the pre- 
vious cases. 
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All of these observations tend to confirm that program 
SOUND is free of errors and producing the correct (and 
expected) results. 

Checking for the Coincidence Effect 

Thus far, no specific evidence has been offered to show 
that the absolute values of the sound power radiated by the 
plates or the vibrational response of the plates predicted 
by program SOUND are correct. Since the radiation effi- 
ciency calculations involve both the acoustic and dynamic 
portions of the analysis, the calculations of a are the 
best candidate for checking to be sure that no shift (up or 
down) has occured in the predicted values of either the 
sound power or the space-time averaged mean-square surface 
velocity. 

It might be argued that the fact that a never exceed- 
ed the value 1 ( 0 on the dB scale) in any of the more than 
10000 separate calculations of radiation efficiency present- 
ed in this paper is very strong evidence in itself that the 
absolute values predicted by program SOUND are correct. The 
close agreement in the radiation efficiency levels predicted 
analytically and those measured for plates 1 through 6 is 
offered as further evidence that the predicted levels are 
correct. In order to further establish the credibility of 
the absolute levels predicted by program SOUND, it was 


258 



decided that a case should be run in which the radiation 
efficiency is expected to exceed the value 1. 

From the discussions contained in reference 48, it is 
clear that the only frequency region in which the radiation 
efficiency a is expected to exceed the value 1 for a 
simply supported flat plate is near the critical frequency. 
Plates 1 through 6 all have critical frequencies which are 
much greater than 1000 Hz (outside the range of analysis). 
Therefore, in order to observe the coincidence phenomenon, 
it was necessary to find (or invent) a plate which has its 
critical frequency in the 0 to 1000 Hz frequency range. 

Since the author was unaware of any material which (when 
made into a plate of comparable size and without the aid of 
stiffeners) would meet this requirement, it was decided to 
invent an imaginary plate with the desired qualities. 
Therefore, program SOUND was run for the case of a plate 
that is the same size as plate no. 1 (.79 mm thick AA 2024 
aluminum) except that it is 5 times thicker and is con- 
structed of an imaginary material. The imaginary material 
has a density that is exactly the same as the aluminum (2.22 
kg/m ) but has moduli of elasticity and rigidity (E and G) 
that are exactly 10 times greater. Simple calculations 
using equation (A2.20) show that this plate has only two 
modes with natural frequencies in the 0 - 2000 Hz range and 
only one of these two occurs within the 0 - 1000 Hz range. 
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Thus, the fundamental frequency of this plate is found to be 
roughly 717 Hz. Using equation (Al.57), the calculated 
critical frequency of this plate is found to be roughly 954 
Hz. Since, in this case, the response of the plate is con- 
trolled almost exclusively by the fundamental mode, there is 
expected to be little or no difference in the airborne and 
structureborne radiation efficiencies of the plate. 

Figures (AX-8) and (AX-9) show the results of the com- 
putations for a 1 Pa peak, normally incident, spatially 
uniform acoustic load, and a . 01 N peak, point vibrational 
load located near the corner of the plate (a^ = .06033 m 
and <*2 = . 13547 m) . Figure (AX-8) shows that, as ex- 
pected, there is virtually no difference in the radiation 
efficiencies resulting from the two forcing functions. 
Futhermore, the radiation efficiency a exceeds the value 1 
(0 dB) slightly in advance of the critical frequency of 954 
Hz as predicted by the theory (see reference 48). Figure 
(AX- 9) shows that the peak noise radiation occurs near the 
fundamental mode at 717 Hz. This result was expected since 
the plate's level of vibration should peak in this range. 
Figure (AX-9) also shows that the acoustic input produces 
more noise radiation than the vibrational input. This 
result was also anticipated since the net force input to the 
plate by the airborne forcing function is roughly ten times 
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the net force that is input by the structureborne forcing 
function . 

The author contends that these results, taken together 
with the results of the directivity studies, and the results 
presented for plates 1 through 6 serve as adequate proof 
that program SOUND is free of bugs. 




Figure (AX-1) — Natural frequencies versus mode number m 

for different values of mode number n for 
plate no. 1. 
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Figure (AX-2) — Plate no. 1 directivity plots at 20 m distance for a 1 Pa peak, normally 

incident, spatially uniform acoustic input with a forcing frequency of 
500 Hz. 
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Figure (AX-3) — Plate no. 1 directivity plots at 20 m distance for a 1 Pa peak, normally 

incident, spatially uniform acoustic input with a forcing frequency of 
1000 Hz. 
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- Plate no. 1 directivity plots at 20 m distance for a . 01 N peak, point 
vibrational input at the center of the plate with a forcing frequency of 
500 Hz. 
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Figure (AX-5) -- Plate no. 1 directivity plots at 20 m distance for a . 01 N peak, point 

vibrational input at the center of the plate with a forcing frequency of 
1000 Hz. 
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Figure (AX-6) -- Plate no. 1 directivity plots at 20 m distance for a . 01 N peak, point 

vibrational input near the corner of the plate with a forcing frequency 
of 500 Hz. 
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Figure (AX-7) — Plate no. 1 directivity plots at 20 m distance for a .01 N peak, point 

vibrational input near the corner of the plate with a forcing frequency 
of 1000 Hz. 
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Figure (AX-8) — Airborne and structureborne radiation efficiencies of an extremely stiff 

plate constructed of an imaginary material. 
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Figure (AX-9) — Airborne and structureborne sound powers radiated by an extremely stiff 

plate constructed of an imaginary material. 


Appendix XI 

TWO MICROPHONE ACOUSTIC INTENSITY MEASUREMENT METHOD 


A_j_ DERIVATION OF THE FUNDAMENTAL EQUATION FOR INTENSITY 
Recall from Appendix I, equation (A1.47) that the 
acoustic intensity is given by 


00 


00 


<J> t = v (0) = 


C (£) 
pv v 


df = Re{ 


G pv (f) 


df} 


(All.l) 


0 0 

and recall that an estimate of G (f) can be obtained 

ir v 

from equation (Al.39) as 


G pv (f) = 2/T 


P*(f)V(f) 


(All. 2) 


Substituting the estimate given by (All. 2) into (All.l) 


<I>. = Re{ 


2/T 

P*(f)V(f) 

0 



df } 


(All. 3) 


Thus, the acoustic intensity as a function of frequency is 
given by 


<l(f)> t = Re{2Af/T 


P*(f)V(f) 


} . 


(All. 4) 


And from the theory of Finite Fourier transforms, it is 
known that Af = 1/T , thus the equation for intensity as 
function of frequency is given by 
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(All. 5) 


<l(f)> t = 2/T 2 Re{P*(f)V(f) } , 

where the Finite Fourier transforms are defined by the equa- 
tions (Al.4) and (A1.5) given in Appendix I, viz. 

T 

+ r — j 2 tt f t 

P(r,f) = F (p(r, t) } = p(r, t) e dt , (All. 6) 

0 

+ -y f — j 2 tt f t 

V(r,f) = F{v(r,t)} = v(r,t) e dt . (All. 7) 

0 

Now recall that the Navier-Stokes equation of momentum 
conservation for a compressible flow is given by 

P Q Dv/Dt = p Q g - Vp + V*T, (All. 8) 

where p Q is the density of the fluid medium, v is the 
particle velocity, g is the acceleration due to gravity, p 
is the acoustic pressure, t is the stress 

tensor due to viscous effects, "V is the gradient operator, 
and D/Dt is the substantial or total derivative given by 

D/Dt = 3/3t + v^3/3x^ + V23/3X2 + v^S/Sx^ . (All. 9) 

It is well known in acoustic theory that the effects of the 
viscous stress tensor on acoustic wave propagation is that 
of a damping factor which is proportional to the acoustic 
driving frequency squared (see reference 59, pg . 518). 

Thus, the damping effects are known to be insignificant at 


272 



low frequency. Neglecting the effects of viscosity, the 
momentum equation simplifies to the familiar Euler's equa- 
tion, viz. 

p Q Dv/Dt = p Q g - Vp • (All. 10) 

Furthermore, gravity has a negligible influence on acoustic 
disturbances except at extremely low driving frequencies 
(see reference 59, pg. 9). Making a small perturbation 
assumption and neglecting the effects of gravity, the momen- 
tum equation becomes 


9v/9t = -l/p Q Vp • 


(All. 11) 


If the analysis is confined to a single dimension, the 
momentum equation becomes 


9v x /9t = -1/p o 9p/9x L . 


(All. 12) 


Making a finite difference approximation for the pressure 
gradient yields 


3v x /9t = -1/ ( p Q Ax) ^ - p 2 ) 


(All. 13) 


where Ax is the spacing between the two microphones. 
Applying the Fourier transform as defined in Appendix I, 
equation (A1.3), the momentum equation becomes 


00 

- jmt 

9v (m)/9t e dt = -l/(p Ax) 

X o 

— oo 


P 1 (to) - P 2 (oj) 


. (All. 14) 
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Integrating the left side of equation (All. 14) by parts, 


DW 


-jut 

v x (w) e dt = -l/(p Q Ax) 


P 1 (u) - P 2 (w) 


. (All. 15) 


Rearranging equation (All. 15) algebraically 


-jut 

v x (u) e dt = j/( P q uAx) 


P x (u) - P 2 (u) 


(All. 16) 


The term on the left side of equation (All. 16) is the 
Fourier transform of the particle velocity. Therefore, 


V(u) = j / ( p Q u Ax ) 


P 1 (u) - P 2 (u) 


(All. 17 ) 


Equation (All. 17) is an approximation of the particle 
velocity at a point midway between two microphones. The 
pressure midway between two microphones can be estimated by 


P(u) = 


P 1 (u) + P 2 (u) 


/2 


(All. 18) 


Substituting equations (All. 17) and (All. 18) into equation 
(All. 5), and letting the frequency dependence be implicit 


<I> t = 2/T 2 Re{P*V} 


= 2/T z Re{ j/(2p Q uAx) 


* * * * 
P 1 P 1 - P 2 P 2 + ? 2 P 1 - P 1 P 2 


} 


(All. 19) 
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Recalling the equations for the estimation of auto and cross 
spectral density functions (Al.38) and (Al.39) and 
substituting into equation (All. 19), the result is 


<I>, 


2/r Re{ ( j/2p Q wAx) • (T/2 ) 


G 11 " G 22 + G 12 ' G 12 


= 1/T Re{ (j/2p Q ajAx) 


5 11 “ G 22 + 2 ^ Q 12 


} . 


} 

(All. 20) 


Taking the real part of the right side of equation (All. 20) 
as indicated, the result is 


<I> t = -1/T Q 12 /(p o <jjAx) = -Q 12 Af/(p o wAx) 

= Im{G 12 }Af/(p o wAx) (All. 21) 


where -Q^ 2 the imaginary part of the one sided cross 
spectral density between microphones 1 and 2 (quad-spectral 
density) . 
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B. ERRORS ASSOCIATED WITH ACOUSTIC INTENSITY MEASUREMENTS 
There are five principle sources of error associated 
with the two microphone cross spectral method of intensity 
measurement in the 0-1000 Hz range. They are as follows: 

1. Instrumentation phase mismatch. 

2. Finite difference error. 

3. Directional effects and errors of interpretation. 

4. Near field effects. 

5. Random errors. 

For the convenience of the reader, a brief discussion of 
each type of error is presented here. More detailed 
discussions are contained in references 30, and 63-65. 

Recall equation (All. 21) which states 

<I>£ = Im{G^2 }Af/(p o uAx) . (All. 21) 

Utilizing equation (Al.39) from Appendix I, the equation for 
time averaged intensity can be written as 

<I> t = Im{ 2/T P*P 2 } Af /( p q oj A x) . (All. 22) 

The complex Fourier transforms in equation (All. 22) can be 
represented in complex polar form and the equation can be 
rewritten as 

j(e,-0 9 ) 

<I> t = Im{ 2/T IpJ |P 2 I e } Af/(p o wAx) 
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PjJ 1 1*2 1 sinCe^e^ /(p q wAx) . 


(All. 23) 


= 2/T 2 


Thus, it is obvious from equation (All. 23) that the magni- 
tude of the intensity vector is proportional to the sine of 
the relative phase difference between the two microphones. 
This measured relative phase difference (0^-02 or 
A0) has a component due to instrumentation phase mis- 
match. The instrumentation phase mismatch causes the larg- 
est errors when it is of the same order of magnitude or 
larger than the physical phase difference. (The physical 
phase difference is the actual phase difference between the 
microphone signals that one hopes to measure.) This leads 
to the conclusion that the instrumentation phase mismatch 
causes its most serious errors at low frequencies where the 
wavelengths are long and the intensity probe is measuring a 
very small phase difference. Hence, this type of error is 
largest primarily in the low frequency regime. Elimination 
of the instrumentation phase mismatch may be approached in 
one of two ways. One method proposed by Chung et. al. (see 
reference 30) uses a microphone interchange technique to 
eliminate this type of error. The more common method is to 
carefully measure the instrumentation phase mismatch and 
compensate for it in subsequent computer calculations. The 
microphone interchange technique was used for the measure- 
ments in this study. This method utilizes two measurements 
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of the cross spectra and uses the following equation to 
calculate the acoustic intensity: 


<I>^ = Im{ 


G 12 * G 21 


1/2 


} Af/(p o u)Ax) 


(All. 24) 


The use of this equation for the intensity measurments 
cancels all instrumentation phase mismatch exactly. 

The second type of error introduced by the two micro- 
phone method is the error associated with the finite differ- 
ence approximation of equation (All. 13). This error occurs 
primarily in the high frequency regime. To assure that this 
error is small, it should be required that 

kAx = w/c Q Ax = 2ir Ax/X << ir/2 , (All. 25) 

so that Ax/X << 1/4 , where Ax is the spacing between 
microphones and X is the wavelength of interest. 

The third type of measurement error stems from misin- 
terpretation of the results. Directional effects and flank- 
ing from multiple sources can result in the measurement of 
components of intensity vectors unintended by the measurer. 
Careful planning and execution of the measurements can help 
to prevent the acquisition of data contaminated with acous- 
tic intensity vector components from unwanted sound sources. 

The fourth type of error is near field measurement 
error. In theory, the large pressure gradients in the near 
field of higher order acoustic sources such as dipoles and 
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quadrapoles can cause considerable error in the measurement 
accuracy of the two microphone method. For a detailed dis- 
cussion of this type of error see reference 63. 

The final category of errors is random or statistical 
errors. Errors in the space averaging and time averaging 
are included in the random error. An extensive analysis of 
the statistical errors encountered in the time averaging of 
acoustic intensity measurements is contained in reference 
64. Random errors due to the selected method of space aver- 
aging are not as well understood. Recent studies (see ref- 
erence 65) indicate that, for the same number of measure- 
ments, the continuous sweep method of space averaging is 
more accurate than the fixed point method of averaging even 
under conditions where the sweep rate is non-constant. 
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Appendix XII 

INSTRUMENT SPECIFICATIONS 
Altec 405-8H Loudspeakers 

Power - 10 Watts music power 

Frequency response - 60-15000 Hz 

Sensitivity - 92 dB SPL at 4 ft. from 1 Watt 

Impedance - 8 Ohms 

Size - 4 3/8 in. X 4 3/8 in. X 2 1/8 in. 

BBN 501 Piezoelectric Accelerometers 
Sensitivity - 10 mV/g 

Amplitude linearity - Within 1% up to 212 g 
Noise floor 

(1 Hz to 25 kHz) - .000020 V 

Frequency response - Flat from 7 to 20000 Hz 
Temp, range - -54 to 121 degrees Celsius 

Peak operating g's - 212 

Weight - 1.8 grams 

Size - 0.31 in. Dia. X 0.42 in. H 
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Bruel Kjaer 2610 Measuring Amplifiers 


Measuring range 
Frequency range 

Phase deviation 

Overall gain 

Preamp input 

Filter mode 

Power required 
Temp . range 
Humidity 

Bruel & Kjaer 2706 
Power output cap. 
Limiting 
Frequency range 
Noise and Hum 
Attenuator 
Gain control 
Power required 
Temp, range 


.00001 V to 30 V FSD 

Flat response 2 Hz to 200 kHz 
(within 0.5 dB) 

Within 5 degrees between any two 
instruments from 20 Hz to 20kHz 

From -30 to +100 dB selectable in 
accurate 10 dB steps plus extra 0 
to 10 dB gain for continuous 
adjustment between steps 

Standard B & K 7 pin microphone 
socket with +200 V polarization 

22.4 Hz high pass filter with 
18 dB/octave fall off 

115 V at 60 Hz 

+5 to +40 degrees Celsius 

0 to 90 % relative humidity 


Power Amplifier 

75 VA into 3 Ohm resistor load 

- 5 A for Vibration Exciter Type 4809 

- 10 Hz to 20 kHz within 0.5 dB 

At least 70 dB below full output 

- 0 to 40 dB in 10 dB steps 

0 to minus infinity logrithmic 
115 V at 60 Hz approx. 140 Watts 
5 to 40 degrees Celsius 
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Bruel & Kjaer 4809 
Rated force 
Frequency range 
Axial resonance 
Max. acceleration 
Max. displacement 
Maximum velocity 
Max. input current 

Bruel & Kjaer 4166 
Measuring range 

Frequency range 
Size 

Fluke 2205A Switch 
Channel capacity 

Isolation 

Cross talk 

Temp, range 
Humidity 
Power required 


Vibration Exciter 

44.5 Newton, 10 lbf sine peak 
10 Hz to 20 kHz bare table 
20 kHz bare table 
736 m/s2 (75 g) 

8 mm (0.315 in) peak-to-peak 
1.65 m/s (65 in/s) peak 
5 A rms 

Condenser Microphone with 2619 preamp. 

25 to 145 dB SPL 

(lin. noise level 20 Hz to 20 kHz) 
Flat response 20 to 20 kHz 
1/2 in. diameter 

Controller (scanner) 

10 channels used in this study 
(expandable to 100 channels) 

- Signal lines to power ground 
greater than 10 billion Ohms 

Less than 30 dB below applied 
signal from DC to 1 MHz when 
terminated with 1 MegaOhm 

0 to 50 degrees Celsius 

0% to 80% relative humidity 

120 V at 60 Hz, 15 VA maximum 


282 



Genrad 1382 Random Noise Generator 


Spectrum 


Output voltage 


Amplitude control 


Power required 


White noise from 20 Hz to 20 kHz 
within 1 dB with 1/2 power points 
at 10 Hz and at 50 kHz 

Maximum open circuit output voltage 
is at least 3 V rms for any 
bandwidth 

Output amplitude is continuously 
adjustable from full output to 
approx. 60 dB below that level 

6W at 100 to 125 V at 60 Hz 


Genrad 2515 Computer-Aided Test System 


Number of channels 
Frequency range 

Dynamic range 
Amplitude flatness 
Channel deviation 

Sensitivity 

Max. Voltage 
Bandwidths 

Storage 

Averaging 


8 (expandable to 16) 

DC to 25.6 kHz with alias 

protection on all channels 

(DC to 1024 Hz used in this study) 

> 70 dB with 12 bit A/D conversion 

Within 0.25 dB over entire range 

Within 0.2 dB amplitude 
< 2.0 degrees phase up to 10 kHz 

8 ranges 

(0.0625 V full scale lowest range) 
(8.0 V full scale highest range) 

Protected to 25 V 

15 selectable ranges 
from 1.28 Hz to 20.5 kHz 
plus 25.6 kHz 

10 Mbytes on Winchester drive and 
0.5 Mbytes on mini-floppy 

Summation averaging with Hanning 
and 50% redundancy used in this 
study 
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Genrad 2515 Computer-Aided Test System (cont'd) 


Processing 

Power required 
Temp . range 
Humidity 

Hewlett Packard 

Attenuation 

Accuracy 

Power capacity 

Hewlett Packard 
Ranges 

Voltage 

measurement 

accuracy 

Response time 


Measurement and calculation of 
transfer functions and auto power 
spectra were used in this study 

120 V at 60 Hz 500 Watts 

10 to 45 degrees Celsius 

20% to 80% relative humidity 


350D Attenuator Set 

110 dB in 10 dB and 1 dB steps 

For the 10 dB section from 0 to 
100 kHz error is less than 0.125 dB 

For the 100 dB section from 0 to 
100 kHz error is less than 0.25 dB 
at any step up to 70 dB and less 
than 0.5 dB for steps above 70 dB 

In 600 Ohms; 5 Watts (55 V rms) 


3403C True RMS Voltmeter 

.01 V 
. 1 V 

1.0 V 

10.0 V 

100.0 V 

1000.0 V 


Error within 0.2% of full scale 
plus 0.2% of reading for the 25 Hz 
to 100 kHz range for the 1 V, 10 V, 
and 100 V settings 

Fast: 1 second 
Slow: 10 seconds 
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Kenwood KA-52 Stereo Integrated Amplifier 


Power output 

Total Harmonic 

Distortion 

(20 Hz to 20 kHz) 

Frequency response 

SNR 

Power required 


55 Watts when both channels driven 
at 8 Ohms 1 kHz 

.09% at rated power into 8 Ohms 
Flat 10 Hz to 100 kHz 
101 dB for 150 mV input 
120 V at 60 Hz 185 Watts 


PCB 221 Force Transducer 

Dynamic range 

Sensitivity 

Resonant frequency 

Full scale ouput 

Temp . range 

Max. vibration 

Rockland 1022F Analog Fi 


.004 to 100 lbs. 

50 mV/lb 
70,000 Hz 
5 volts 

-73 to 121 degrees Celsius 
2000 g 


Filter type 

Frequency range 
Roll-off 
Responses 
Passband gains 


2 analog filters which can be 
selected as either high pass or 
low pass 

10 Hz to 111 kHz 
24 db/octa ve/sec tion 
Butterworth, RC 
Selectable 0, 20, 40 dB 
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Tektronix T922R Oscilloscope 

Type - Dual channel 

Frequency range - DC to 5 MHz 

Max. input voltage - 30 V 

Power required - 120 V at 60 Hz 

Temp, range - 0 to 45 degrees Celsius 
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End of Document 



